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ABSTRACT 

The formation of the first stars and quasars marks the transformation of the universe 
from its smooth initial state to its clumpy current state. In popular cosmological models, 
the first sources of light began to form at a redshift z = 30 and reionized most of the 
hydrogen in the universe by z = 7. Current observations are at the threshold of probing 
the hydrogen reionization epoch. The study of high-redshift sources is likely to attract 
major attention in observational and theoretical cosmology over the next decade. 
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1. Preface: The Frontier of Small-Scale Structure 

The detection of cosmic microwave background (CMB) anisotropies (Bennett et al. 1996; de 
Bernardis et al. 2000; Hanany et al. 2000) confirmed the notion that the present large-scale structure 
in the universe originated from small- amplitude density fluctuations at early times. Due to the 
natural instability of gravity, regions that were denser than average collapsed and formed bound 
objects, first on small spatial scales and later on larger and larger scales. The present-day abundance 
of bound objects, such as galaxies and X-ray clusters, can be explained based on an appropriate 
extrapolation of the detected anisotropies to smaller scales. Existing observations with the Hubble 
Space Telescope (e.g., Steidel et al. 1996; Madau et al. 1996; Chen et al. 1999; Clements et al. 1999) 
and ground-based telescopes (Lowenthal et al. 1997; Dey et al. 1999; Hu et al. 1998, 1999; Spinrad 
et al. 1999; Steidel et al. 1999), have constrained the evolution of galaxies and their stellar content at 
z <^ 6. However, in the bottom-up hierarchy of the popular Cold Dark Matter (CDM) cosmologies, 
galaxies were assembled out of building blocks of smaller mass. The elementary building blocks, i.e., 
the first gaseous objects to form, acquired a total mass of order the Jeans mass (~ 10 4 M Q ), below 
which gas pressure opposed gravity and prevented collapse (Couchman & Rees 1986; Haiman & 
Loeb 1997; Ostriker & Gnedin 1996). In variants of the standard CDM model, these basic building 
blocks first formed at z ~ 15-30. 

An important qualitative outcome of the microwave anisotropy data is the confirmation that 
the universe started out simple. It was by and large homogeneous and isotropic with small fluctu- 
ations that can be described by linear perturbation analysis. The current universe is clumpy and 
complicated. Hence, the arrow of time in cosmic history also describes the progression from simplic- 
ity to complexity (see Figure |l]). While the conditions in the early universe can be summarized on 
a single sheet of paper, the mere description of the physical and biological structures found in the 
present-day universe cannot be captured by thousands of books in our libraries. The formation of 
the first bound objects marks the central milestone in the transition from simplicity to complexity. 
Pedagogically, it would seem only natural to attempt to understand this epoch before we try to 
explain the present-day universe. Historically, however, most of the astronomical literature focused 
on the local universe and has only been shifting recently to the early universe. This violation of the 
pedagogical rule was forced upon us by the limited state of our technology; observation of earlier 
cosmic times requires detection of distant sources, which is feasible only with large telescopes and 
highly-sensitive instrumentation. 

For these reasons, advances in technology are likely to make the high redshift universe an 
important frontier of cosmology over the coming decade. This effort will involve large (30 meter) 
ground-based telescopes and will culminate in the launch of the successor to the Hubble Space 
Telescope, called Next Generation Space Telescope (NGST). Figure || shows an artist's illustration 
of this telescope which is currently planned for launch in 2009. NGST will image the first sources of 
light that formed in the universe. With its exceptional sub-nJy (1 nJy = 10 _32 erg cm -2 s~ x Hz^ 1 ) 
sensitivity in the 1-3.5/xm infrared regime, NGST is ideally suited for probing optical-UV emission 
from sources at redshifts ^ 10, just when popular Cold Dark Matter models for structure formation 
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Fig. 1. — Milestones in the evolution of the universe from simplicity to complexity. The "end of the 
dark ages" bridges between the recombination epoch probed by microwave anisotropy experiments 
(z ~ 10 3 ) and the horizon of current observations (z ~ 5-6). 



- 6 - 



predict the first baryonic objects to have collapsed. 




Fig. 2. — Artist's illustration of one of the current designs (GSFC) of the Next Generation Space 
Telescope. More details about the telescope can be found at |http:/ /ngst.gsfc.nasa.gov/ 

The study of the formation of the first generation of sources at early cosmic times (high 
redshifts) holds the key to constraining the power-spectrum of density fluctuations on small scales. 
Previous research in cosmology has been dominated by studies of Large Scale Structure (LSS); 
future studies are likely to focus on Small Scale Structure (SSS). 

The first sources are a direct consequence of the growth of linear density fluctuations. As 
such, they emerge from a well-defined set of initial conditions and the physics of their formation 
can be followed precisely by computer simulation. The cosmic initial conditions for the formation 
of the first generation of stars are much simpler than those responsible for star formation in the 
Galactic interstellar medium at present. The cosmic conditions are fully specified by the primordial 
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power spectrum of Gaussian density fluctuations, the mean density of dark matter, the initial 
temperature and density of the cosmic gas, and the primordial composition according to Big-Bang 
nucleosynthesis. The chemistry is much simpler in the absence of metals and the gas dynamics 
is much simpler in the absence of both dynamically-significant magnetic fields and feedback from 
luminous objects. 

The initial mass function of the first stars and black holes is therefore determined by a simple 
set of initial conditions (although subsequent generations of stars are affected by feedback from 
photoionization heating and metal enrichment). While the early evolution of the seed density fluc- 
tuations can be fully described analytically, the collapse and fragmentation of nonlinear structure 
must be simulated numerically. The first baryonic objects connect the simple initial state of the 
universe to its complex current state, and their study with hydrodynamic simulations (e.g., Abel et 
al. 1998a, Abel, Bryan, & Norman 2000; Bromm, Coppi, & Larson 1999) and with future telescopes 
such as NGST offers the key to advancing our knowledge on the formation physics of stars and 
massive black holes. 

The first light from stars and quasars ended the "dark ages" ^ of the universe and initiated 
a "renaissance of enlightenment" in the otherwise fading glow of the microwave background (see 
Figure |l]) . It is easy to see why the mere conversion of trace amounts of gas into stars or black holes 
at this early epoch could have had a dramatic effect on the ionization state and temperature of the 
rest of the gas in the universe. Nuclear fusion releases ~ 7x 10 6 eV per hydrogen atom, and thin-disk 
accretion onto a Schwarzschild black hole releases ten times more energy; however, the ionization 
of hydrogen requires only 13.6 eV. It is therefore sufficient to convert a small fraction, ~ 10~ 5 of 
the total baryonic mass into stars or black holes in order to ionize the rest of the universe. (The 
actual required fraction is higher by at least an order of magnitude [Bromm, Kudritzky, & Loeb 
2000] because only some of the emitted photons are above the ionization threshold of 13.6 eV and 
because each hydrogen atom recombines more than once at redshifts z ^ 7). Recent calculations of 
structure formation in popular CDM cosmologies imply that the universe was ionized at z ~ 7-12 
(Haiman & Loeb 1998, 1999b,c; Gnedin & Ostriker 1997; Chiu & Ostriker 2000; Gnedin 2000a), 
and has remained ionized ever since. Current observations are at the threshold of probing this 
epoch of reionization, given the fact that galaxies and quasars at redshifts ~ 6 are being discovered 
(Fan et al. 2000; Stern et al. 2000). One of these sources is a bright quasar at z = 5.8 whose 
spectrum is shown in Figure ||. The plot indicates that there is transmitted flux shortward of the 
Lya wavelength at the quasar redshift. The optical depth at these wavelengths of the uniform 
cosmic gas in the intergalactic medium is however (Gunn & Peterson 1965), 



m e cH(z s ) ' V0.03;V0.3y V 10 

where H ~ lOO/i km s _1 Mpc -1 ^^ (1 + z s ) 3 ^ 2 is the Hubble parameter at the source redshift 



2 The use of this term in the cosmological context was coined by Sir Martin Rees. 
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f a = 0.4162 and X a = 1216A are the oscillator strength and the wavelength of the Lya transition; 
nHi(z s ) is the neutral hydrogen density at the source redshift (assuming primordial abundances); 
Q m and Qf, are the present-day density parameters of all matter and of baryons, respectively; and 
xm is the average fraction of neutral hydrogen. In the second equality we have implicitly considered 
high redshifts (see equations (|9|) and (|i~0|) in § |2.1| ). Modeling of the transmitted flux (Fan et al. 
2000) implies r s < 0.5 or xui <^ 10~ 6 , i.e., the low-density gas throughout the universe is fully 
ionized at z = 5.8! One of the important challenges for future observations will be to identify when 
and how the intergalactic medium was ionized. Theoretical calculations (see § 6.3.1[) imply that such 
observations are just around the corner. 

Figure ||] shows schematically the various stages in a theoretical scenario for the history of 
hydrogen reionization in the intergalactic medium. The first gaseous clouds collapse at redshifts 
~ 20-30 and fragment into stars due to molecular hydrogen (H2) cooling. However, H2 is fragile 
and can be easily dissociated by a small flux of UV radiation. Hence the bulk of the radiation 
that ionized the universe is emitted from galaxies with a virial temperature ;> 10 4 K, where atomic 



cooling is effective and allows the gas to fragment (see the end of §3.3 for an alternative scenario 



Since recent observations confine the standard set of cosmological parameters to a relatively 
narrow range, we assume a ACDM cosmology with a particular standard set of parameters in 
the quantitative results in this review. For the contributions to the energy density, we assume 
ratios relative to the critical density of £l m = 0.3, S1a = 0.7, and Qf, = 0.045, for matter, vacuum 
(cosmological constant), and baryons, respectively. We also assume a Hubble constant Hq = 
100 h km s _1 Mpc _1 with h = 0.7, and a primordial scale invariant (n = 1) power spectrum with 
as = 0.9, where is the root-mean-square amplitude of mass fluctuations in spheres of radius 
8 h^ 1 Mpc. These parameter values are based primarily on the following observational results: 
CMB temperature anisotropy measurements on large scales (Bennett et al. 1996) and on the scale 
of ~ 1° (Lange et al. 2000; Balbi et al. 2000); the abundance of galaxy clusters locally (Viana & 
Liddle 1999; Pen 1998; Eke, Cole, & Frenk 1996) and as a function of redshift (Bahcall & Fan 1998; 
Eke, Cole, Frenk, & Henry 1998); the baryon density inferred from big bang nucleosynthesis (see 
the review by Tytler et al. 2000) ; distance measurements used to derive the Hubble constant (Mould 
et al. 2000; Jha et al. 1999; Tonry et al. 1997; but see Theureau et al. 1997; Parodi et al. 2000); 
and indications of cosmic acceleration from distances based on type la supernovae (Perlmutter et 
al. 1999; Riess et al. 1998). 

This review summarizes recent theoretical advances in understanding the physics of the first 
generation of cosmic structures. Although the literature on this subject extends all the way back 
to the sixties (Saslaw &: Zipoy 1967, Peebles & Dicke 1968, Hirasawa 1969, Matsuda et al. 1969, 
Hutchins 1976, Silk 1983, Palla et al. 1983, Lepp & Shull 1984, Couchman 1985, Couchman & Rees 
1986, Lahav 1986), this review focuses on the progress made over the past decade in the modern 
context of CDM cosmologies. 
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Fig. 3. — Optical spectrum of the highest-redshift known quasar at z = 5.8, discovered by the 
Sloan Digital Sky Survey (Fan et al. 2000). 
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Fig. 4. — Stages in the reionization of hydrogen in the intergalactic medium. 
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2. Hierarchical Formation of Cold Dark Matter Halos 

2.1. The Expanding Universe 

The modern physical description of the universe as a whole can be traced back to Einstein, who 
argued theoretically for the so-called "cosmological principle" : that the distribution of matter and 
energy must be homogeneous and isotropic on the largest scales. Today isotropy is well established 
(see the review by Wu, Lahav, & Rees 1999) for the distribution of faint radio sources, optically- 
selected galaxies, the X-ray background, and most importantly the cosmic microwave background 
(henceforth, CMB; see, e.g., Bennett et al. 1996). The constraints on homogeneity are less strict, 
but a cosmological model in which the universe is isotropic but significantly inhomogeneous in 
spherical shells around our special location is also excluded (Goodman 1995). 

In General Relativity, the metric for a space which is spatially homogeneous and isotropic is 
the Robertson- Walker metric, which can be written in the form 



ds 2 = dt 2 - a\t) 



^ -o2 , -2 a JJ.2 X 



1-kR 2 



R 2 (dO 



+ R" dti* + sin" 



(2) 



where a(t) is the cosmic scale factor which describes expansion in time, and (R, 0, <p) are spherical 
comoving coordinates. The constant k determines the geometry of the metric; it is positive in a 
closed universe, zero in a flat universe, and negative in an open universe. Observers at rest remain 
at rest, at fixed (R, 9, (ft), with their physical separation increasing with time in proportion to a(t). 
A given observer sees a nearby observer at physical distance D receding at the Hubble velocity 
H(t)D, where the Hubble constant at time t is H(t) = dlna(t)/dt. Light emitted by a source at 
time t is observed at t = with a redshift z = l/a(t) — 1, where we set a(t = 0) = 1. 

The Einstein field equations of General Relativity yield the Friedmann equation (e.g., Weinberg 
1972; Kolb & Turner 1990) 

H 2 {t) = —p-^, (3) 



which relates the expansion of the universe to its matter-energy content. For each component of 
the energy density p, with an equation of state p = p(p), the density p varies with a(t) according 
to the equation of energy conservation 

d( P R 3 ) = -pd(R 3 ) . (4) 

With the critical density 

fdt) S (5) 

defined as the density needed for k = 0, we define the ratio of the total density to the critical 
density as 

n = — . (6) 

PC 
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With Q m , Oa, and Q r denoting the present contributions to f2 from matter (including cold dark 
matter as well as a contribution from baryons), vacuum density (cosmological constant), and 
radiation, respectively, the Friedmann equation becomes 



H(t) 



H, 



o 



O O O I 1 / 2 

—3- + "A H t H ; 

a 13 a 4 a 



(7) 



where we define Hq and f2o = £l m + JIa + to be the present values of H and O, respectively, and 
we let 

n k = -A = i - o . (8) 

In the particularly simple Einstein-de Sitter model (O m = 1, ^a = = ^fc = 0), the scale factor 
varies as a(t) oc i 2 / 3 . Even models with non-zero J7a or fl k approach the Einstein-de Sitter behavior 
at high redshifts, i.e., when 



'1 + z) 3> max 



(i - n m - n A )/n m , (^ A /^ m ) 1/3 l (9) 



(as long as fi r can be neglected). The Friedmann equation implies that models with f2/% = 
converge to the Einstein-de Sitter limit faster than do open models. E.g., for f2 m = 0.3 and 
SI A = 0.7 equation (S) corresponds to the condition z 3> 1.3, which is easily satisfied by the 
reionization redshift. In this high-z regime, H(t) ~ 2/(3t), and the age of the universe is 

t ^ 2 -= (1 + z)^ 2 = 5.38 x 10 8 ( ±±±) ^ yr , (10) 



3 Hq \j V 10 

where in the last expression we assumed our standard cosmological parameters (see the end of §|l]). 

In the standard hot Big Bang model, the universe is initially hot and the energy density is 
dominated by radiation. The transition to matter domination occurs at z ~ 10 4 , but the universe 
remains hot enough that the gas is ionized, and electron-photon scattering effectively couples the 
matter and radiation. At z ~ 1200 the temperature drops below ~ 3300 K and protons and 
electrons recombine to form neutral hydrogen. The photons then decouple and travel freely until 
the present, when they are observed as the CMB. 



2.2. Linear Gravitational Growth 

Observations of the CMB (e.g., Bennett et al. 1996) show that the universe at recombination 
was extremely uniform, but with spatial fluctuations in the energy density and gravitational po- 
tential of roughly one part in 10 5 . Such small fluctuations, generated in the early universe, grow 
over time due to gravitational instability, and eventually lead to the formation of galaxies and the 
large-scale structure observed in the present universe. 

As in the previous section, we distinguish between fixed and comoving coordinates. Using vec- 
tor notation, the fixed coordinate r corresponds to a comoving position x = r/a. In a homogeneous 
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universe with density p, we describe the cosmological expansion in terms of an ideal pressure-less 
fluid of particles each of which is at fixed x, expanding with the Hubble flow v = H(t)r where 
v = dr/dt. Onto this uniform expansion we impose small perturbations, given by a relative density 
perturbation 

«M = P -f - i , (ID 

where the mean fluid density is p, with a corresponding peculiar velocity u = v — Hr. Then the 
fluid is described by the continuity and Euler equations in comoving coordinates (Peebles 1980, 
1993): 

^ + Iv-[(l + *)u] = (12) 

^+#u + i(u-V)u = —V<f>. (13) 
at a a 

The potential <p is given by the Poisson equation, in terms of the density perturbation: 

V 2 (j) = 4TrGpa 2 5 . (14) 

This fluid description is valid for describing the evolution of collisionless cold dark matter particles 
until different particle streams cross. This "shell-crossing" typically occurs only after perturbations 
have grown to become non-linear, and at that point the individual particle trajectories must in 
general be followed. Similarly, baryons can be described as a pressure-less fluid as long as their 
temperature is negligibly small, but non-linear collapse leads to the formation of shocks in the gas. 

For small perturbations 5 <C 1, the fluid equations can be linearized and combined to yield 

+ 2H*=«JAJ. (15) 

This linear equation has in general two independent solutions, only one of which grows with time. 
Starting with random initial conditions, this "growing mode" comes to dominate the density evolu- 
tion. Thus, until it becomes non-linear, the density perturbation maintains its shape in comoving 
coordinates and grows in proportion to a growth factor D(t). The growth factor is in general given 
by (Peebles 1980) 

(n A a 3 + n k a + n m ) 1/2 /■« a 3 / 2 da 
D(t) oc ^ / tjk , (16) 

« 3/2 J (n A a3 + n k a + n m f 2 

where we neglect O r when considering halos forming at z <C 10 4 . In the Einstein-de Sitter model 
(or, at high redshift, in other models as well) the growth factor is simply proportional to a(t). 

The spatial form of the initial density fluctuations can be described in Fourier space, in terms 
of Fourier components 

<5 k = / d 3 x5{x)e- ikx . (17) 



-14- 



Here we use the comoving wavevector k, whose magnitude k is the comoving wavenumber which is 
equal to 2tt divided by the wavelength. The Fourier description is particularly simple for fluctuations 
generated by inflation (e.g., Kolb & Turner 1990). Inflation generates perturbations given by 
a Gaussian random field, in which different k-modes are statistically independent, each with a 
random phase. The statistical properties of the fluctuations are determined by the variance of the 
different k-modes, and the variance is described in terms of the power spectrum P(k) as follows: 



In standard models, inflation produces a primordial power-law spectrum P(k) oc k n with n ~ 1. 
Perturbation growth in the radiation-dominated and then matter-dominated universe results in a 
modified final power spectrum, characterized by a turnover at a scale of order the horizon cH^ 1 at 
matter-radiation equality, and a small-scale asymptotic shape of P{k) oc k n ~ 4 . On large scales the 
power spectrum evolves in proportion to the square of the growth factor, and this simple evolution 
is termed linear evolution. On small scales, the power spectrum changes shape due to the additional 
non-linear gravitational growth of perturbations, yielding the full, non-linear power spectrum. The 
overall amplitude of the power spectrum is not specified by current models of inflation, and it is 
usually set observationally using the CMB temperature fluctuations or local measures of large-scale 
structure. 

Since density fluctuations may exist on all scales, in order to determine the formation of 
objects of a given size or mass it is useful to consider the statistical distribution of the smoothed 
density field. Using a window function W(y) normalized so that / d 3 yW(y) = 1, the smoothed 
density perturbation field, / d 3 y<5(x + y)W(y), itself follows a Gaussian distribution with zero 
mean. For the particular choice of a spherical top-hat, in which W = 1 in a sphere of radius R and 
is zero outside, the smoothed perturbation field measures the fluctuations in the mass in spheres 
of radius R. The normalization of the present power spectrum is often specified by the value of 
eg = a(R = 8h~ 1 Mpc). For the top-hat, the smoothed perturbation field is denoted Sr or 8m, 
where the mass M is related to the comoving radius R by M = 4irp m R 3 /3, in terms of the current 
mean density of matter p m . The variance {5m) 2 is 




(18) 



where 5^ is the three-dimensional Dirac delta function. 




(19) 



where ji(x) = (sinx — xcosx)/x 2 . The function a(M) plays a crucial role in estimates of the 
abundance of collapsed objects, as described below. 



2.3. Formation of Nonlinear Objects 



The small density fluctuations evidenced in the CMB grow over time as described in the previ- 
ous subsection, until the perturbation 5 becomes of order unity, and the full non-linear gravitational 
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problem must be considered. The dynamical collapse of a dark matter halo can be solved analyti- 
cally only in cases of particular symmetry. If we consider a region which is much smaller than the 
horizon cH^ 1 , then the formation of a halo can be formulated as a problem in Newtonian gravity, 
in some cases with minor corrections coming from General Relativity. The simplest case is that of 
spherical symmetry, with an initial (t = U <C to) top-hat of uniform overdensity Si inside a sphere 
of radius R. Although this model is restricted in its direct applicability, the results of spherical 
collapse have turned out to be surprisingly useful in understanding the properties and distribution 
of halos in models based on cold dark matter. 

The collapse of a spherical top-hat is described by the Newtonian equation (with a correction 
for the cosmological constant) 

d 2 r TT n^ GM 

^2=^0^-—, (20) 

where r is the radius in a fixed (not comoving) coordinate frame, Hq is the present Hubble constant, 
M is the total mass enclosed within radius r, and the initial velocity field is given by the Hubble 
flow dr/dt = H(t)r. The enclosed 5 grows initially as 5l = 5iD(t)/D(ti), in accordance with linear 
theory, but eventually S grows above 5l- If the mass shell at radius r is bound (i.e., if its total 
Newtonian energy is negative) then it reaches a radius of maximum expansion and subsequently 
collapses. At the moment when the top-hat collapses to a point, the overdensity predicted by linear 
theory is (Peebles 1980) 5l = 1.686 in the Einstein-de Sitter model, with only a weak dependence 
on Q m and Thus a top-hat collapses at redshift z if its linear overdensity extrapolated to the 
present day (also termed the critical density of collapse) is 

A ( \ L686 (01\ 

«*) = , (21) 

where we set D(z = 0) = 1. 

Even a slight violation of the exact symmetry of the initial perturbation can prevent the top- 
hat from collapsing to a point. Instead, the halo reaches a state of virial equilibrium by violent 
relaxation (phase mixing). Using the virial theorem U = —2K to relate the potential energy U to 
the kinetic energy K in the final state, the final overdensity relative to the critical density at the 
collapse redshift is A c = 187T 2 ~ 178 in the Einstein-de Sitter model, modified in a universe with 
Q m + f^A = 1 to the fitting formula (Bryan & Norman 1998) 

A c = 18vr 2 + 82d - 39d 2 , (22) 

where d = Q.^ — 1 is evaluated at the collapse redshift, so that 

= n m (i + z f 

m n m (i + z) 3 + n A + n k (i + z) 2 ' 1 ; 

A halo of mass M collapsing at redshift z thus has a (physical) virial radius 

M 







- 1 / 3 (l + z\ 




18vr 2 _ 


V io J 



h- 1 kpc , (24) 
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and a corresponding circular velocity, 
(GM\ l l 2 



Ve. 



23.4 
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10 s h- 1 Mr. 
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V 10 J 



1/2 



km s 



(25) 



In these expressions we have assumed a present Hubble constant written in the form Hq 
100 h km s^Mpc -1 . We may also define a virial temperature 



2k K 



1.98 x 10 4 
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(26) 



where \i is the mean molecular weight and m p is the proton mass. Note that the value of \x depends 
on the ionization fraction of the gas; [i = 0.59 for a fully ionized primordial gas, /i = 0.61 for a 
gas with ionized hydrogen but only singly-ionized helium, and fi = 1.22 for neutral primordial gas. 
The binding energy of the halo is approximately!] 
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Note that the binding energy of the baryons is smaller by a factor equal to the baryon fraction 

Although spherical collapse captures some of the physics governing the formation of halos, 
structure formation in cold dark matter models proceeds hierarchically. At early times, most of the 
dark matter is in low-mass halos, and these halos continuously accrete and merge to form high-mass 
halos. Numerical simulations of hierarchical halo formation indicate a roughly universal spherically- 
averaged density profile for the resulting halos (Navarro, Frenk, White 1997, hereafter NFW), 
though with considerable scatter among different halos (e.g., Bullock et al. 2000). The NFW profile 
has the form 

Kr) = S(l + ^-^— ^, (28) 



8ttG cnx(1 + c^x) 2 

where x = r/r v i r , and the characteristic density 5 C is related to the concentration parameter cn by 



A, 



3 ln(l + cn) - c N / (1 + On) 



(29) 



The concentration parameter itself depends on the halo mass M , at a given redshift z. We note that 
the dense, cuspy halo profile predicted by CDM models is not apparent in the mass distribution 
derived from measurements of the rotation curves of dwarf galaxies (e.g., de Blok & McGaugh 1997; 
Salucci &: Burkert 2000), although observational and modeling uncertainties may preclude a firm 
conclusion at present (van den Bosch et al. 2000; Swaters, Madore, & Trewhella 2000). 



3 The coefficient of 1/2 in equation ( pTj ) would be exact for a singular isothermal sphere, p(r) oc 1/r 2 
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2.4. The Abundance of Dark Matter Halos 



In addition to characterizing the properties of individual halos, a critical prediction of any the- 
ory of structure formation is the abundance of halos, i.e., the number density of halos as a function 
of mass, at any redshift. This prediction is an important step toward inferring the abundances of 
galaxies and galaxy clusters. While the number density of halos can be measured for particular 
cosmologies in numerical simulations, an analytic model helps us gain physical understanding and 
can be used to explore the dependence of abundances on all the cosmological parameters. 

A simple analytic model which successfully matches most of the numerical simulations was 
developed by Press & Schechter (1974). The model is based on the ideas of a Gaussian random 
field of density perturbations, linear gravitational growth, and spherical collapse. To determine 
the abundance of halos at a redshift z, we use 5m, the density field smoothed on a mass scale M, 



as defined in §2.2. Although the model is based on the initial conditions, it is usually expressed 
in terms of redshift-zero quantities. Thus, we use the linearly-extrapolated density field, i.e., the 
initial density field at high redshift extrapolated to the present by simple multiplication by the 
relative growth factor (see § |2.2f) . Similarly, in this section the 'present power spectrum' refers 
to the initial power spectrum, linearly-extrapolated to the present without including non-linear 
evolution. Since 5m is distributed as a Gaussian variable with zero mean and standard deviation 
a(M) [which depends only on the present power spectrum, see equation fll9|)1, the probability that 
5m is greater than some 5 equals 



OO [ 

c15m j — exp 

8 V2na(M) 



2a 2 (M) 



Hvfc) ■ <30) 



The fundamental ansatz is to identify this probability with the fraction of dark matter particles 
which are part of collapsed halos of mass greater than M, at redshift z. There are two additional 
ingredients: First, the value used for 5 is 5 cr it(z) given in equation (|2l|), which is the critical density 
of collapse found for a spherical top-hat (extrapolated to the present since a(M) is calculated using 
the present power spectrum); and second, the fraction of dark matter in halos above M is multiplied 
by an additional factor of 2 in order to ensure that every particle ends up as part of some halo with 
M > 0. Thus, the final formula for the mass fraction in halos above M at redshift z is 



F(> M\z) =erfcf%^ 



(31) 



This ad-hoc factor of 2 is necessary, since otherwise only positive fluctuations of 5m would 
be included. Bond et al. (1991) found an alternate derivation of this correction factor, using a 
different ansatz. In their derivation, the factor of 2 has a more satisfactory origin, namely the 
so-called "cloud-in-cloud" problem: For a given mass M, even if 5m is smaller than 5 cr i t (z), it 
is possible that the corresponding region lies inside a region of some larger mass Ml > M, with 
5 Ml > 5 CT it(z). In this case the original region should be counted as belonging to a halo of mass Ml- 
Thus, the fraction of particles which are part of collapsed halos of mass greater than M is larger 



than the expression given in equation (|30|). Bond et al. showed that, under certain assumptions, 
the additional contribution results precisely in a factor of 2 correction. 



Differentiating the fraction of dark matter in halos above M yields the mass distribution. 
Letting dn be the comoving number density of halos of mass between M and M + dM, we have 



where v c = S CT n(z)/a(M) is the number of standard deviations which the critical collapse over- 
density represents on mass scale M. Thus, the abundance of halos depends on the two functions 
a(M) and S cr it( z ), each of which depends on the energy content of the universe and the values of 
the other cosmological parameters. We illustrate the abundance of halos for our standard choice of 
the ACDM model with Q m = 0.3 (see the end of §|). 

Figure || shows a(M) and S CI it(z), with the input power spectrum computed from Eisenstein 
& Hu (1999). The solid line is a{M) for the cold dark matter model with the parameters specified 
above. The horizontal dotted lines show the value of <5 C rit(-z) at z = 0, 2, 5, 10, 20 and 30, as indicated 
in the figure. From the intersection of these horizontal lines with the solid line we infer, e.g., that 
at z = 5 a 1 - a fluctuation on a mass scale of 2 x 10 7 Mq will collapse. On the other hand, at 
z = 5 collapsing halos require a 2 — a fluctuation on a mass scale of 3x 10 10 Mq, since a{M) on this 
mass scale equals about half of 5 cr it{z = 5). Since at each redshift a fixed fraction (31.7%) of the 
total dark matter mass lies in halos above the 1 — a mass, Figure |5| shows that most of the mass is 
in small halos at high redshift, but it continuously shifts toward higher characteristic halo masses 
at lower redshift. Note also that a(M) flattens at low masses because of the changing shape of the 
power spectrum. Since a — > oo as M — » 0, in the cold dark matter model all the dark matter is 
tied up in halos at all redshifts, if sufficiently low-mass halos are considered. 

Also shown in Figure [5] is the effect of cutting off the power spectrum on small scales. The short- 
dashed curve corresponds to the case where the power spectrum is set to zero above a comoving 
wavenumber k = lOMpc^ 1 , which corresponds to a mass M = 1.7 x 1O 8 M0. The long-dashed 
curve corresponds to a more radical cutoff above k = IMpc -1 , or below M = 1.7 x 10 11 Mq. A 
cutoff severely reduces the abundance of low-mass halos, and the finite value of a{M = 0) implies 
that at all redshifts some fraction of the dark matter does not fall into halos. At high redshifts 
where S CT it(z) ^ cr(M = 0), all halos are rare and only a small fraction of the dark matter lies in 
halos. In particular, this can affect the abundance of halos at the time of reionization, and thus the 
observed limits on reionization constrain scenarios which include a small-scale cutoff in the power 
spectrum (Barkana, Haiman, & Ostriker 2000). 

In Figures |6] - |9| we show explicitly the properties of collapsing halos which represent 1 — a, 
2 — a, and 3 — a fluctuations (corresponding in all cases to the curves in order from bottom to top), 
as a function of redshift. No cutoff is applied to the power spectrum. Figure |6| shows the halo mass, 
Figure [7| the virial radius, Figure ^ the virial temperature (with [x in equation ( p6| ) set equal to 
0.6, although low temperature halos contain neutral gas) as well as circular velocity, and Figure ^ 




(32) 
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shows the total binding energy of these halos. In Figures |6] and |8[ the dashed curves indicate 
the minimum virial temperature required for efficient cooling (see § |3.3| ) with primordial atomic 
species only (upper curve) or with the addition of molecular hydrogen (lower curve). Figure || 
shows the binding energy of dark matter halos. The binding energy of the baryons is a factor 
~ Qb/£l m ~ 15% smaller, if they follow the dark matter. Except for this constant factor, the figure 
shows the minimum amount of energy that needs to be deposited into the gas in order to unbind 
it from the potential well of the dark matter. For example, the hydrodynamic energy released by 
a single super novae, ~ 10 51 erg, is sufficient to unbind the gas in all 1 — a halos at z ;> 5 and in all 
2 — cj halos at z ^ 12. 




10 a 10 4 10 6 10 8 10 10 10 12 10 14 
M[Solar Mass] 



Fig. 5. — Mass fluctuations and collapse thresholds in cold dark matter models. The horizontal dot- 
ted lines show the value of the extrapolated collapse overdensity <5 C rit(#) at the indicated redshifts. 
Also shown is the value of a(M) for the cosmological parameters given in the text (solid curve), as 
well as a(M) for a power spectrum with a cutoff below a mass M = 1.7 X 1O 8 M (short-dashed 
curve), or M = 1.7 x 10 11 Mq (long-dashed curve). The intersection of the horizontal lines with 
the other curves indicate, at each redshift z, the mass scale (for each model) at which a 1 — a 
fluctuation is just collapsing at z (see the discussion in the text). 
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Redshift 

Fig. 6. — Characteristic properties of collapsing halos: Halo mass. The solid curves show the mass 
of collapsing halos which correspond to 1 — a, 2 — a, and 3 — a fluctuations (in order from bottom 
to top). The dashed curves show the mass corresponding to the minimum temperature required for 
efficient cooling with primordial atomic species only (upper curve) or with the addition of molecular 
hydrogen (lower curve). 
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Fig. 7. — Characteristic properties of collapsing halos: Halo virial radius. The curves show the 
virial radius of collapsing halos which correspond to 1 — a, 2 — a, and 3 — a fluctuations (in order 
from bottom to top). 
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Fig. 8. — Characteristic properties of collapsing halos: Halo virial temperature and circular velocity. 
The solid curves show the virial temperature (or, equivalently, the circular velocity) of collapsing 
halos which correspond to 1 — a, 2 — a, and 3 — a fluctuations (in order from bottom to top). The 
dashed curves show the minimum temperature required for efficient cooling with primordial atomic 
species only (upper curve) or with the addition of molecular hydrogen (lower curve). 
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Fig. 9. — Characteristic properties of collapsing halos: Halo binding energy. The curves show the 
total binding energy of collapsing halos which correspond to 1 — a, 2 — a, and 3 — a fluctuations 
(in order from bottom to top). 
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Fig. 10. — Halo mass function at several redshifts: z = (solid curve), z = 5 (dotted curve), z = 10 
(short-dashed curve), z = 20 (long-dashed curve), and z = 30 (dot-dashed curve). 
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At z = 5, the halo masses which correspond to 1 — a, 2 — a, and 3 — a fluctuations are 
1.8 x 10 Mq, 3.0 x 1O 1O M0, and 7.0 x 10 11 Mq, respectively. The corresponding virial temperatures 
are 2.0 x 10 3 K, 2.8 x 10 5 K, and 2.3 x 10 6 K. The equivalent circular velocities are 7.5 km s _1 , 88 
km s" 1 , and 250 km s _1 . At z = 10, the 1 — a, 2 — a, and 3 — a fluctuations correspond to halo 
masses of 1.3 x 10 3 Mq, 5.7 x 10 7 Mq, and 4.8 x 10 9 Mq, respectively. The corresponding virial 
temperatures are 6.2 K, 7.9 x 10 3 K, and 1.5 x 10 5 K. The equivalent circular velocities are 0.41 
km s~ , 15 km s , and 65 km s" 1 . Atomic cooling is efficient at T v i r ;> 10 4 K, or a circular velocity 
;> 17 km s~ 1 . This corresponds to a 1.2 — a fluctuation and a halo mass of 2.1 x 10 8 Mq at z = 5, 
and a 2.1 — a fluctuation and a halo mass of 8.3 x 10 7 Mq at z = 10. Molecular hydrogen provides 
efficient cooling down to T V - 1T ~ 300 K, or a circular velocity V c ~ 2.0 km s -1 . This corresponds to 
a 0.76 — a fluctuation and a halo mass of 3.5 x 10 5 Mq at z = 5, and a 1.3 — a fluctuation and a 
halo mass of 1.4 x 1O 5 M at z = 10. 

In Figure [ITJ we show the halo mass function dn/dln(M) at several different redshifts: z = 
(solid curve), z = 5 (dotted curve), z = 10 (short-dashed curve), z = 20 (long-dashed curve), and 
z = 30 (dot-dashed curve). Note that the mass function does not decrease monotonically with 
redshift at all masses. At the lowest masses, the abundance of halos is higher at z > than at 
z = 0. 

3. Gas Infall and Cooling in Dark Matter Halos 

3.1. Cosmological Jeans Mass 

The Jeans length Aj was originally defined (Jeans 1928) in Newtonian gravity as the critical 
wavelength that separates oscillatory and exponentially-growing density perturbations in an infinite, 
uniform, and stationary distribution of gas. On scales t smaller than Aj, the sound crossing 
time, i/c s is shorter than the gravitational free-fall time, (Gp) _1//2 , allowing the build-up of a 
pressure force that counteracts gravity. On larger scales, the pressure gradient force is too slow 
to react to a build-up of the attractive gravitational force. The Jeans mass is defined as the mass 
within a sphere of radius Aj/2, Mj = (4-7r/3) / o(Aj/2) 3 . In a perturbation with a mass greater than 
Mj, the self-gravity cannot be supported by the pressure gradient, and so the gas is unstable to 
gravitational collapse. The Newtonian derivation of the Jeans instability suffers from a conceptual 
inconsistency, as the unperturbed gravitational force of the uniform background must induce bulk 
motions (compare Binney & Tremaine 1987). However, this inconsistency is remedied when the 
analysis is done in an expanding universe. 

The perturbative derivation of the Jeans instability criterion can be carried out in a cos- 
mological setting by considering a sinusoidal perturbation superposed on a uniformly expanding 
background. Here, as in the Newtonian limit, there is a critical wavelength Aj that separates oscil- 
latory and growing modes. Although the expansion of the background slows down the exponential 
growth of the amplitude to a power-law growth, the fundamental concept of a minimum mass that 
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can collapse at any given time remains the same (see, e.g. Kolb & Turner 1990; Peebles 1993). 

We consider a mixture of dark matter and baryons with density parameters r2| m = p~dm/Pc 
and = Pb/Pc where pdm is the average dark matter density, pb is the average baryonic density, 
p c is the critical density, and + $7^ = fi^ is given by equation (|2^). We also assume spatial 
fluctuations in the gas and dark matter densities with the form of a single spherical Fourier mode 
on a scale much smaller than the horizon, 

Pdm(r,t) - Pdmjt) 
Pdm(t) 

Pb(r,t) - p h (t) 

Pb(t) 

where p~dm(t) and pb(t) are the background densities of the dark matter and baryons, <5dm(i) and 
5b(t) are the dark matter and baryon overdensity amplitudes, r is the comoving radial coordinate, 
and k is the comoving perturbation wavenumber. We adopt an ideal gas equation-of-state for the 
baryons with a specific heat ratio 7=5/3. Initially, at time t = the gas temperature is uniform 
Tb(r,ti)=Ti, and the perturbation amplitudes are small <5dm,i><5b,i ^ 1- We define the region inside 
the first zero of s'm(kr)/(kr), namely < kr < ir, as the collapsing "object". 

The evolution of the temperature of the baryons Tb(r, t) in the linear regime is determined by 
the coupling of their free electrons to the Cosmic Microwave Background (CMB) through Compton 
scattering, and by the adiabatic expansion of the gas. Hence, Tb(r, t) is generally somewhere 
between the CMB temperature, T 7 oc (1 + and the adiabatically-scaled temperature T a< j oc 
(1 + z)~ 2 . In the limit of tight coupling to T 7 , the gas temperature remains uniform. On the other 
hand, in the adiabatic limit, the temperature develops a gradient according to the relation 

T b oc^- 1} . (35) 



8dm(t) 

Sb(t) 



sm(kr) 
kr 

sin(fcr) 
kr 



(33) 
(34) 



The evolution of dark matter overdensity, <5d m (^)) m the linear regime is described by the 
equation (see §9.3.2 of Kolb & Turner 1990), 

3 

<W + 2i/(5 dm = -H (filA + ^dm^dm) (36) 

whereas the evolution of the overdensity of the baryons, 5b(t), is described by 

5 b + 2H5b = ^H 2 (n h 5 h + n dm 5 Am )-^ - U U b + -(3[8 h - S h 4 ] . (37) 

Here, H(t) = a/a is the Hubble parameter at a cosmological time t, and p = 1.22 is the mean 
molecular weight of the neutral primordial gas in atomic units. The parameter (3 distinguishes 
between the two limits for the evolution of the gas temperature. In the adiabatic limit (3 = 1, 
and when the baryon temperature is uniform and locked to the background radiation, /3 = 0. 
The last term on the right hand side (in square brackets) takes into account the extra pressure 
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gradient force in V{p\ i T) = (TV px> + p^VT), arising from the temperature gradient which develops 
in the adiabatic limit. The Jeans wavelength Aj = 27r/fcj is obtained by setting the right-hand 
side of equation (37) to zero, and solving for the critical wavenumber kj. As can be seen from 



equation (|37|), the critical wavelength Aj (and therefore the mass Mj) is in general time-dependent. 
We infer from equation ( |37| ) that as time proceeds, perturbations with increasingly smaller initial 
wavelengths stop oscillating and start to grow. 



To estimate the Jeans wavelength, we equate the right-hand-side of equation ( [371) to zero. 
We further approximate 5^ ~ 5dmi an d consider sufficiently high redshifts at which the universe 
is matter-dominated and flat (equations (||) and (10) in §2T). We also assume fib <C fi m > where 



fi m = fidm + fife is the total matter density parameter. Following cosmological recombination 
at z ~ 10 3 , the residual ionization of the cosmic gas keeps its temperature locked to the CMB 
temperature (via Compton scattering) down to a redshift of (p. 179 of Peebles 1993) 

l + 137(fi 6 /i 2 /0.022) 2/5 . (38) 

In the redshift range between recombination and Zt, (3 = and 

kj = (2tt/Aj) = [2fc B T 7 (0)/3 / um p ]- 1 / 2 v / fi^^o , (39) 

so that the Jeans mass is therefore redshift independent and obtains the value (for the total mass 
of baryons and dark matter) 



/U, = ( — )%(0) = 1.35 x 10 5 ( ' M & . (40) 



Based on the similarity of Mj to the mass of a globular cluster, Peebles & Dicke (1968) sug- 
gested that globular clusters form as the first generation of baryonic objects shortly after cosmo- 
logical recombination. Peebles & Dicke assumed a baryonic universe, with a nonlinear fluctuation 
amplitude on small scales at z ~ 10 3 , a model which has by now been ruled out. The lack of a 
dominant mass of dark matter inside globular clusters (Moore 1996; Heggie & Hut 1995) makes it 
unlikely that they formed through direct cosmological collapse, and more likely that they resulted 
from fragmentation during the process of galaxy formation. Furthermore, globular clusters have 
been observed to form in galaxy mergers (e.g., Miller et al. 1997). 

At z <; zt, the gas temperature declines adiabatically as [(1 + z)/(l + zt)] 2 (i.e., (3 = 1) and 
the total Jeans mass obtains the value, 

Mj = 5.73 x 10 3 — — M . (41) 

J ^ 0.15 ) \o.o22 y v io y v ; 

Note that we have neglected Compton drag, i.e., the radiation force which suppresses gravitational 
growth of structure in the baryon fluid as long as the electron abundance is sufficiently high to 
keep the baryons dynamically coupled to the photons. After cosmological recombination, the net 
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friction force on the predominantly neutral fluid decreases dramatically, allowing the baryons to fall 
into dark matter potential wells, and essentially erasing the memory of Compton drag by z ~ 100 
(e.g., §5.3.1. of Hu 1995). 

It is not clear how the value of the Jeans mass derived above relates to the mass of collapsed, 
bound objects. The above analysis is perturbative (Eqs. |36| and |37j] are valid only as long as <5b and 
5dm are much smaller than unity), and thus can only describe the initial phase of the collapse. As 
5^ and £d m grow and become larger than unity, the density profiles start to evolve and dark matter 
shells may cross baryonic shells (Haiman, Thoul, & Loeb 1996) due to their different dynamics. 
Hence the amount of mass enclosed within a given baryonic shell may increase with time, until 
eventually the dark matter pulls the baryons with it and causes their collapse even for objects 
below the Jeans mass. 

Even within linear theory, the Jeans mass is related only to the evolution of perturbations at 
a given time. When the Jeans mass itself varies with time, the overall suppression of the growth 
of perturbations depends on a time-averaged Jeans mass. Gnedin & Hui (1998) showed that the 
correct time-averaged mass is the filtering mass Mp = (4-7r/3) p (2-7ra/kp) 3 , in terms of the comoving 
wavenumber kp associated with the "filtering scale". The wavenumber kp is related to the Jeans 
wavenumber fcj by 

J_ J_ /■* ,/ 2(ff Mt') + 2H(t')D(t') ft dt" 
k F (t) D(t)Jo K) kj(t>) J t > a?(t") 

where D(t) is the linear growth factor (§ |2.2|) . At high redshift (where 
simplifies to (Gnedin 2000b) 

1 _ 3 f a da' 
k 2 F (t) ~ a Jo k 2 j{a') 

Then the relationship between the linear overdensity of the dark matter 5^ m and the linear over- 
density of the baryons 5b, in the limit of small k, can be written as (Gnedin & Hui 1998) 

t- l - k i + °^- (44) 

Linear theory specifies whether an initial perturbation, characterized by the parameters k, 
^dm,i) <5b,i and t{, begins to grow. To determine the minimum mass of nonlinear baryonic objects 
resulting from the shell-crossing and virialization of the dark matter, we must use a different model 
which examines the response of the gas to the gravitational potential of a virialized dark matter 
halo. 



(42) 

1), this relation 
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3.2. Response of Baryons to Nonlinear Dark Matter Potentials 

The dark matter is assumed to be cold and to dominate gravity, and so its collapse and 
virialization proceeds unimpeded by pressure effects. In order to estimate the minimum mass of 
baryonic objects, we must go beyond linear perturbation theory and examine the baryonic mass 
that can accrete into the final gravitational potential well of the dark matter. 

For this purpose, we assume that the dark matter had already virialized and produced a gravi- 
tational potential (f)(r) at a redshift z N -„ (with eft — > at large distances, and 4> < inside the object) 
and calculate the resulting overdensity in the gas distribution, ignoring cooling (an assumption jus- 
tified by spherical collapse simulations which indicate that cooling becomes important only after 
virialization; see Haiman, Thoul, & Loeb 1996). 

After the gas settles into the dark matter potential well, it satisfies the hydrostatic equilibrium 
equation, 

Vp b = -PbV<A (45) 

where and pb are the pressure and mass density of the gas. At z <^ 100 the gas temperature 
is decoupled from the CMB, and its pressure evolves adiabatically (ignoring atomic or molecular 
cooling), 



Pb = ( Pb\ 
Pb \PbJ 



5/3 

(46) 



where a bar denotes the background conditions. We substitute equation ( |4q ) into (45) and get the 
solution, 

p b _f 1 2 M rryM 3 / 2 (4?) 



Pb V 5 k B f J 

where T = p~\,\im v / '{k-Qph) is the background gas temperature. If we define T v ; r = — ^fj,m p <p/kB 
as the virial temperature for a potential depth —(f), then the overdensity of the baryons at the 
virialization redshift is 

*-£-i-(i + §£)*-,. 

This solution is approximate for two reasons: (i) we assumed that the gas is stationary throughout 
the entire region and ignored the transitions to infall and the Hubble expansion at the interface 
between the collapsed object and the background intergalactic medium (henceforth IGM), and (ii) 
we ignored entropy production at the virialization shock surrounding the object. Nevertheless, the 
result should provide a better estimate for the minimum mass of collapsed baryonic objects than 
the Jeans mass does, since it incorporates the nonlinear potential of the dark matter. 

We may define the threshold for the collapse of baryons by the criterion that their mean 
overdensity, <5b, exceeds a value of 100, amounting to ^ 50% of the baryons that would assemble in 
the absence of gas pressure, according to the spherical top-hat collapse model (§|2.3|). Equation (H8j) 
then implies that T v ; r > 17.2 T. 
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As mentioned before, the gas temperature evolves at z <^ 160 according to the relation T ~ 
170[(1 + z)/100] 2 K. This implies that baryons are overdense by 5^ > 100 only inside halos with 
a virial temperature T vir ;> 2.9 x 10 3 [(1 + z)/100] 2 K. Based on the top-hat model (p^), this 
implies a minimum halo mass for baryonic objects of 

Mm ^, 0xw3 (^Y'' 2 (^Y 3l5 (l^r M S , (49) 



v 0.15 J ^0.022 y V 10 

where we set [i = 1.22 and consider sufficiently high redshifts so that fi^j ~ 1. This minimum mass is 
coincidentally almost identical to the naive Jeans mass calculation of linear theory in equation (El]) 
despite the fact that it incorporates shell crossing by the dark matter, which is not accounted for by 
linear theory. Unlike the Jeans mass, the minimum mass depends on the choice for an overdensity 
threshold [taken arbitrarily as 5^ > 100 in equation (|49|)]. To estimate the minimum halo mass 
which produces any significant accretion we set, e.g., 5b = 5, and get a mass which is lower than 
M m i n by a factor of 27. 

Of course, once the first stars and quasars form they heat the surrounding IGM by either 
outflows or radiation. As a result, the Jeans mass which is relevant for the formation of new 
objects changes (Ostriker & Gnedin 1997; Gnedin 2000a). The most dramatic change occurs when 
the IGM is photo-ionized and is consequently heated to a temperature of ~ (1-2) x 10 4 K. As we 



discuss in §6.5, this heating episode had a dramatic impact on galaxy formation. 



3.3. Molecular Chemistry, Photo-Dissociation, and Cooling 

Before metals are produced, the primary molecule which acquires sufficient abundance to affect 
the thermal state of the pristine cosmic gas is molecular hydrogen, H_2. The dominant H2 formation 
process is 

H + e~ -> H~ + hv, (50) 
H- + H H 2 + e- (51) 

where free electrons act as catalysts. The complete set of chemical reactions leading to the formation 
of H-2 is summarized in Table 1, together with the associated rate coefficients (see also Haiman, 
Thoul, & Loeb 1996; Abel et al. 1997; Galli & Palla 1998; and the review by Abel & Haiman 2000). 
Table 2 shows the same for deuterium mediated reactions. Due to the low gas density, the chemical 
reactions are slow and the molecular abundance is far from its value in chemical equilibrium. After 
cosmological recombination the fractional H2 abundance is small, ~ 10~ 6 relative to hydrogen 
by number (Lepp & Shull 1984; Shapiro, Giroux k Babul 1994). At redshifts z < 100, the gas 
temperature in most regions is too low for collisional ionization to be effective, and free electrons 
(over and above the residual electron fraction) are mostly produced through photoionization of 
neutral hydrogen by UV or X-ray radiation. 
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In objects with baryonic masses ;> 3 x 10 4 Mq, gravity dominates and results in the bottom-up 
hierarchy of structure formation characteristic of CDM cosmologies; at lower masses, gas pressure 
delays the collapse. The first objects to collapse are those at the mass scale that separates these 
two regimes. Such objects reach virial temperatures of several hundred degrees and can fragment 
into stars only through cooling by molecular hydrogen (e.g., Abel 1995; Tegmark et al. 1997). In 
other words, there are two independent minimum mass thresholds for star formation: the Jeans 
mass (related to accretion) and the cooling mass. For the very first objects, the cooling threshold 
is somewhat higher and sets a lower limit on the halo mass of ~ 5 x 10 4 M Q at z ~ 20. 

However, molecular hydrogen (H2) is fragile and can easily be photo-dissociated by photons 
with energies of 11. 26-13. 6eV, to which the IGM is transparent even before it is ionized. The 
photo-dissociation occurs through a two-step process, first suggested by Solomon in 1965 (compare 
Field et al. 1966) and later analyzed quantitatively by Stecher & Williams (1967). Haiman, Rees, 
& Loeb (1997) evaluated the average cross-section for this process between 11.26eV and 13.6eV, 
by summing the oscillator strengths for the Lyman and Werner bands of H2, and obtained a 
value of 3.71 x 10~ 18 cm 2 . They showed that the UV flux capable of dissociating H2 throughout 
the collapsed environments in the universe is lower by more than two orders of magnitude than 
the minimum flux necessary to ionize the universe. The inevitable conclusion is that soon after 
trace amounts of stars form, the formation of additional stars due to H2 cooling is suppressed. 
Further fragmentation is possible only through atomic line cooling, which is effective in objects 
with much higher virial temperatures, T v i r ^ 10 4 K. Such objects correspond to a total mass ^ 
1O 8 M [(1 + z)/10]~ 3 / 2 . Figure | illustrates this sequence of events by describing two classes of 
objects: those with T vir < 10 4 K (small dots) and those with T v ; r > 10 4 K (large dots). In the first 
stage (top panel), some low- mass objects collapse, form stars, and create ionized hydrogen (H II) 
bubbles around them. Once the UV background between 11.2-13.6eV reaches a specific critical 
level, H2 is photo-dissociated throughout the universe and the formation of new stars is delayed 
until objects with T YU ;> 10 4 K collapse (Haiman, Abel, & Rees 2000; Ciardi, Ferrara, &; Abel 2000; 
Ciardi et al. 2000). Machacek, Bryan & Abel (2000) have confirmed that the soft UV background 
can delay the cooling and collapse of low-mass halos (~ 1O 6 M ) based on analytical arguments 
and three-dimensional hydrodynamic simulations; they also determined the halo mass threshold 
for collapse for a range of UV fluxes. Omukai & Nishi (1999; see also Silk 1977) have argued that 
the photo-dissociation of H2 could be even more effective due to a small number of stars embedded 
within the gas clouds themselves. 

When considering the photo-dissociation of H2 before reionization, it is important to incor- 
porate the processed spectrum of the UV background at photon energies below the Lyman limit. 
Due to the absorption at the Lyman-series resonances this spectrum obtains the sawtooth shape 
shown in Figure |TI|. For any photon energy above Lya at a particular redshift, there is a limited 
redshift interval beyond which no contribution from sources is possible because the corresponding 
photons are absorbed through one of the Lyman-series resonances along the way. Consider, for 
example, an energy of 11 eV at an observed redshift z = 10. Photons received at this energy would 
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have to be emitted at the 12.1 eV Ly/3 line from z = 11.1. Thus, sources in the redshift interval 
10-11.1 could be seen at 11 eV, but radiation emitted by sources at z > 11.1 eV would have passed 
through the 12.1 eV energy at some intermediate redshift, and would have been absorbed. Thus, 
an observer viewing the universe at any photon energy above Lya would see sources only out to 
some horizon, and the size of that horizon would depend on the photon energy. The number of 
contributing sources, and hence the total background flux at each photon energy, would depend on 
how far this energy is from the nearest Lyman resonance. Most of the photons absorbed along the 
way would be re-emitted at Lya and then redshifted to lower energies. The result is a sawtooth 
spectrum for the UV background before reionization, with an enhancement below the Lya energy 
(see Haiman et al. 1997 for more details). Unfortunately, the direct detection of the redshifted 
sawtooth spectrum as a remnant of the reionization epoch is not feasible due to the much higher 
flux contributed by foreground sources at later cosmic times. 

The radiative feedback on H2 need not be only negative, however. In the dense interiors of gas 
clouds, the formation rate of H2 could be accelerated through the production of free electrons by 
X-rays. This effect could counteract the destructive role of H2 photo-dissociation (Haiman, Rees, 
& Loeb 1996). Haiman, Abel, & Rees (2000) have shown that if a significant (^ 10%) fraction 
of the early UV background is produced by massive black holes (mini-quasars) with hard spectra 
extending to photon energies ~ 1 keV, then the X-rays will catalyze H2 production and the net 
radiative feedback will be positive, allowing low mass objects to fragment into stars. These objects 
may greatly alter the topology of reionization (§|6.3|). However, if such quasars do not exist or if low 



mass objects are disrupted by supernova-driven winds (see §7.2), then most of the stars will form 
inside objects with virial temperatures ^ 10 4 K, where atomic cooling dominates. Figure 12 and 
Table 3 summarize the cooling rates as a function of gas temperature in high-redshift, metal-free 
objects. 
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Fig. 11. — The average spectrum during the initial phase of the reionization epoch. The upper 
panel shows that absorption by neutral hydrogen and helium suppresses the flux above 13.6eV 
up to the keV range. The lower panel shows a close-up of the sawtooth modulation due to line 
absorption below 13.6 eV. A constant comoving density of sources was assumed, with each source 
emitting a power-law continuum, which would result in the spectrum shown by the dashed lines if 
absorption were not taken into account. 
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Fig. 12. — Cooling rates as a function of temperature for a primordial gas composed of atomic hydro- 
gen and helium, as well as molecular hydrogen, in the absence of any external radiation. We assume 
a hydrogen number density nn = 0.045 cm -3 , corresponding to the mean density of virialized halos 
at z = 10. The plotted quantity A/n 2 H is roughly independent of density (unless nu 3> 10 cm -3 ), 
where A is the volume cooling rate (in erg/sec/cm 3 ). The solid line shows the cooling curve for an 
atomic gas, with the characteristic peaks due to collisional excitation of H I and He II. The dashed 
line (calculated using the code of Abel available at http:/ /logy.harvard.edu/tabel/PGas/cool.html ) 
shows the additional contribution of molecular cooling, assuming a molecular abundance equal to 
0.1% of n H . 
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Table 1. Reaction Rates for Hydrogen Species 





Reaction 


Rate Coefficient 
(cmV 1 ) 


Reference 


(1) 


H + e" H+ + 2e- 


5.85 x 10" n r 1/2 exp(-157, 809.1/T)(1 + TlJ 2 )- 1 


1 


(2) 


H+ + e~ -> H + hv 


8.40 x 10 _11 T- 1/2 T 3 _0 - 2 (1 + Tg- 7 )- 1 


1 


(3) 


H + eT -> H~ + hv 


See expression in reference 


2 


(4) 


H + H" -> H 2 + e" 


1.30 x 10~ 9 


1 


(5) 


HT + H+ -> 2H 


7.00 x 10~ 7 T~ 1/2 


1 


(6) 


H 2 + e- -> H + H" 


2.70 x 10" 8 T" 3/2 exp(-43,000/T) 


1 


(7) 


H 2 + H -> 3H 


See expression in reference 


1 


(8) 


H 2 + H+ H+ + H 


2.40 x 10~ 9 exp(-21,200/T) 


1 


(9) 


H 2 + e" -> 2H + e" 


4.38 x 10~ 10 exp(-102,000/T)T°' 35 


1 


(10) 


H" + e~ -> H + 2e~ 


4.00 x 10- 12 Texp(-8750/T) 


1 


(11) 


H" + H -» 2H + e" 


5.30 x 10" 20 Texp(-8750/T) 


1 


(12) 


H" + H+ -> H+ + e" 


See expression in reference 


1 



References. — (1) Haiman, Thoul, & Loeb 1996; (2) Abel, et al. 1997. 



Table 2. Reaction Rates for Deuterium Species 



Rate Coefficient 





Reaction 




(cm^- 1 ) 


Reference 


(1) 


D+ + e~ -> D + hv 


8.40 x 10" 


-11 T -1/2 T -0.2 (1+T 0. 7) -1 


1 


(2) 


D + H+ -> D+ + H 


3.70 x 10" 


10 T ' 2S exp(-43/T) 


3 


(3) 


D+ + H -> D + H+ 


3.70 x 10" 


-10y0.28 


3 


(4) 


D+ + H 2 H+ + HD 


2.10 x 10" 


-9 


3 


(5) 


HD + H+ -> H 2 + D+ 


1.00 x 10" 


9 exp(-464/T) 


3 



References. - 



(1) Haiman, Thoul, & Loeb 1996; (3) Galli & Palla 1998. 
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Table 3. Radiative Cooling Processes in the Primordial Gas 





Cooling due to 


Cooling rate 
(erg s" 1 cm~ 3 ) 


Reference 


(1) 


Molecular hydrogen 


See expression in reference 


1 


(2) 


Deuterium hydride (HD) 


See expression in reference 


2 


(3) 


Atomic H & He 


See expression in reference 


3 


(4) 


Compton scattering 


5.6 x 1(T 36 (1 + z) 4 n e (T - T C mb) 


4 



Note. — T is the gas temperature in K, T 3 = T/10 3 K, T 5 = T/10 5 K, T 6 = T/10 6 K, n e 
is the density of free electrons, z is the redshift, and Tomb = 2.73 (1 + z) K is the temperature 
of the CMB. 

References. - (1) Galli & Palla 1998; (2) Flower, Le Bourlot, Pineau des Forets, & 
Roueff 2000; (3) Cen 1992; Verner & Ferland 1996; Ferland et al. 1992; Voronov 1997 (4) 
Ikeuchi & Ostriker 1986. 
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4. Fragmentation of the First Gaseous Objects 
4.1. Star Formation 

Fragmentation into Stars 

As mentioned in the preface, the fragmentation of the first gaseous objects is a well-posed 
physics problem with well specified initial conditions, for a given power-spectrum of primordial 
density fluctuations. This problem is ideally suited for three-dimensional computer simulations, 
since it cannot be reliably addressed in idealized ID or 2D geometries. 

Recently, two groups have attempted detailed 3D simulations of the formation process of 
the first stars in a halo of ~ 10 6 Mq by following the dynamics of both the dark matter and 
the gas components, including H2 chemistry and cooling (Deuterium is not expected to play a 
significant role; Bromm 2000). Bromm et al. (1999) have used a Smooth Particle Hydrodynamics 
(SPH) code to simulate the collapse of a top-hat overdensity with a prescribed solid-body rotation 
(corresponding to a spin parameter A = 5%) and additional small perturbations with P(k) oc /c~ 3 
added to the top-hat profile. Abel et al. (2000) isolated a high-density filament out of a larger 
simulated cosmological volume and followed the evolution of its density maximum with exceedingly 
high resolution using an Adaptive Mesh Refinement (AMR) algorithm. 

The generic results of Bromm et al. (1999; see also Bromm 2000) are illustrated in Figure [l3|. 
The collapsing region forms a disk which fragments into many clumps. The clumps have a typical 
mass ~ 1O 2 -1O 3 M . This mass scale corresponds to the Jeans mass for a temperature of ~ 500K 
and the density ~ 10 4 cm -3 where the gas lingers because its cooling time is longer than its 
collapse time at that point (see Figure 0). This characteristic density is determined by the fact 
that hydrogen molecules reach local thermodynamic equilibrium at this density. At lower densities, 
each collision leads to an excited state and to radiative cooling, so the overall cooling rate is 
proportional to the collision rate, and the cooling time is inversely proportional to the gas density. 
Above the density of ~ 10 cm -3 , however, the relative occupancy of each excited state is fixed at 
the thermal equilibrium value (for a given temperature) , and the cooling time is nearly independent 
of density (e.g., Lepp & Shull 1983). Each clump accretes mass slowly until it exceeds the Jeans 
mass and collapses at a roughly constant temperature (i.e., isothermally) due to H2 cooling. The 
clump formation efficiency is high in this simulation due to the synchronized collapse of the overall 
top-hat perturbation. 

Bromm (2000, Chapter 7) has simulated the collapse of one of the above-mentioned clumps 
with ~ 1OOOM and demonstrated that it does not tend to fragment into sub-components. Rather, 
the clump core of ~ IOOMq free-falls towards the center leaving an extended envelope behind 
with a roughly isothermal density profile. At very high gas densities, three-body reactions become 
important in the chemistry of H2. Omukai & Nishi (1998) have included these reactions as well as 
radiative transfer and followed the collapse in spherical symmetry up to stellar densities. Radiation 
pressure from nuclear burning at the center is unlikely to reverse the infall as the stellar mass builds 
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Fig. 13. — Numerical results from Bromm et al. (1999), showing gas properties at z = 31.2 for 
a collapsing slightly inhomogeneous top-hat region with a prescribed solid-body rotation. Each 
point in the figure is a gas particle in the simulation, (a) Free electron fraction (by number) 
vs. hydrogen number density (in cm -3 ). At densities exceeding n ~ 10 3 cm~ 3 , recombination is 
very efficient, and the gas becomes almost completely neutral, (b) Molecular hydrogen fraction 
vs. number density. After a quick initial rise, the H2 fraction approaches the asymptotic value of 
/ ~ 10~ 3 , due to the H~ channel, (c) Gas temperature vs. number density. At densities below 
~ 1 cm~ 3 , the gas temperature rises because of adiabatic compression until it reaches the virial 
value of T v i r ~ 5000 K. At higher densities, cooling due to H2 drives the temperature down again, 
until the gas settles into a quasi- hydrostatic state at T ~ 500 K and n ~ 10 4 cm -3 . Upon further 
compression due to accretion and the onset of gravitational collapse, the gas shows a further modest 
rise in temperature, (d) Jeans mass (in Mq) vs. number density. The Jeans mass reaches a value 
of Mj ~ 1O 3 M for the quasi-hydrostatic gas in the center of the potential well, and reaches the 
resolution limit of the simulation, M res ~ 2OOM , for densities close to n = 10 8 cm~ 3 . 
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Fig. 14. — Gas and clump morphology at z = 28.9 in the simulation of Bromm et al. (1999). Top 
row: The remaining gas in the diffuse phase. Bottom row: Distribution of clumps. The numbers 
next to the dots denote clump mass in units of Mq. Left panels: Face-on view. Right panels: Edge- 
on view. The length of the box is 30 pc. The gas has settled into a flattened configuration with 
two dominant clumps of mass close to 20, OOOM . During the subsequent evolution, the clumps 
survive without merging, and grow in mass only slightly by accretion of surrounding gas. 
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up. These calculations indicate that each clump may end up as a single massive star; however, it 
is possible that angular momentum or nuclear burning may eventually halt the monolithic collapse 
and lead to further fragmentation. 



The Jeans mass (§ |3.1[ ), which is defined based on small fluctuations in a background of uniform 
density, does not strictly apply in the context of collapsing gas cores. We can instead use a slightly 
modified critical mass known as the Bonnor-Ebert mass (Bonnor 1956; Ebert 1955). For baryons 
in a background of uniform density pf,, perturbations are unstable to gravitational collapse in a 
region more massive than the Jeans mass 

1 / kT \ 3/2 

Mj = 2.9—1 . (52) 

yfpi \Gpm p J 

Instead of a uniform background, we consider a spherical, non-singular, isothermal, self-gravitating 
gas in hydrostatic equilibrium, i.e., a centrally-concentrated object which more closely resembles 
the gas cores found in the above-mentioned simulations. We consider a finite sphere in equilibrium 
with an external pressure. In this case, small fluctuations are unstable and lead to collapse if the 
sphere is more massive than the Bonnor-Ebert mass Mbej given by the same expression as equation 
d52| ) but with a different coefficient (1.2 instead of 2.9) and with p^ denoting in this case the gas 
(volume) density at the surface of the sphere. 

In their simulation, Abel et al. (2000) adopted the actual cosmological density perturbations 
as initial conditions. The simulation focused on the density peak of a filament within the IGM, and 
evolved it to very high densities (Figure [l^). Following the initial collapse of the filament, a clump 
core formed with ~ 200M©, amounting to only ~ 1% of the virialized gas mass. Subsequently due 
to slow cooling, the clump collapsed subsonically in a state close to hydrostatic equilibrium (see 



Figure 16). Unlike the idealized top-hat simulation of Bromm et al. (2000), the collapse of the 
different clumps within the filament is not synchronized. Once the first star forms at the center of 
the first collapsing clump, it is likely to affect the formation of other stars in its vicinity. 

If the clumps in the above simulations end up forming individual very massive stars, then 
these stars will likely radiate copious amounts of ionizing radiation (Carr, Bond, & Arnett 1984; 
Tumlinson & Shull 2000; Bromm et al. 2000) and expel strong winds. Hence, the stars will have 
a large effect on their interstellar environment, and feedback is likely to control the overall star 
formation efficiency. This efficiency is likely to be small in galactic potential wells which have 
a virial temperature lower than the temperature of photoionized gas, ~ 10 4 K. In such potential 
wells, the gas may go through only a single generation of star formation, leading to a "suicidal" 
population of massive stars. 

The final state in the evolution of these stars is uncertain; but if their mass loss is not too 
extensive, then they are likely to end up as black holes (Bond, Carr, & Arnett 1984; Fryer, Woosley, 
& Heger 2001). The remnants may provide the seeds of quasar black holes (Larson 1999). Some 
of the massive stars may end their lives by producing 7-ray bursts. If so then the broad-band 
afterglows of these bursts could provide a powerful tool for probing the epoch of reionization 
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Fig. 15. — Zooming in on the core of a star forming region with the Adaptive Mesh Refinement 
simulation of Abel et al. (2000). The panels show different length scales, decreasing clockwise by 
an order of magnitude between adjacent panels. Note the large dynamic range of scales which are 
being resolved, from 6 kpc (top left panel) down to 10,000 AU (bottom left panel). 
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Fig. 16. — Gas profiles from the simulation of Abel et al. (2000). The cell size on the finest 
grid corresponds to 0.024 pc, while the simulation box size corresponds to 6.4 kpc. Shown are 
spherically-averaged mass-weighted profiles around the baryon density peak shortly before a well 
defined fragment forms (z = 19.1). Panel (a) shows the baryonic number density (solid line), 
enclosed gas mass in solar mass (thin solid line with circles) , and the local Bonnor-Ebert mass Mbe 
(dashed line; see text). Panel (b) plots the molecular hydrogen fraction (by number) /h 2 (solid 
line) and the free electron fraction x (dashed line). The H2 cooling time, tu 2 , the time it takes a 
sound wave to travel to the center, t CTOSS , and the free-fall time ts = [37r/(32G/?)] 1 / 2 are given in 
panel (c). Panel (d) gives the temperature in K as a function of radius. The bottom panel gives 
the local sound speed, c s (solid line with circles), the root-mean-square radial velocities of the dark 
matter (dashed line) and the gas (dashed line with asterisks) as well as the root-mean-square gas 
velocity (solid line with square symbols). The vertical dotted line indicates the radius (~ 5 pc) at 
which the gas has reached its minimum temperature allowed by H2 cooling. The virial radius of 
the 5.6 x 1O 6 M halo is 106 pc. 



-43- 



(Lamb & Reichart 2000; Ciardi & Loeb 2000). There is no better way to end the dark ages than 
with 7-ray burst fireworks. 

Where are the first stars or their remnants located today? The very first stars formed in rare 
high-o" peaks and hence are likely to populate the cores of present-day galaxies (White & Springel 
1999). However, the star clusters which formed in low-cr peaks at later times are expected to behave 
similarly to the collisionless dark matter particles and populate galaxy halos (Loeb 1998). 



4-1.2. Emission Spectrum of Metal- Free Stars 

The evolution of metal-free (Population III) stars is qualitatively different from that of enriched 
(Population I and II) stars. In the absence of the catalysts necessary for the operation of the 
CNO cycle, nuclear burning does not proceed in the standard way. At first, hydrogen burning 
can only occur via the inefficient PP chain. To provide the necessary luminosity, the star has 
to reach very high central temperatures (T c ~ 10 81 K). These temperatures are high enough for 
the spontaneous turn-on of helium burning via the triple-a process. After a brief initial period of 
triple-a burning, a trace amount of heavy elements forms. Subsequently, the star follows the CNO 
cycle. In constructing main-sequence models, it is customary to assume that a trace mass fraction 
of metals (Z ~ 10 -9 ) is already present in the star (El Eid et al. 1983; Castellani et al. 1983). 

Figures |l7] and ^ show the luminosity L vs. effective temperature T for zero-age main sequence 
stars in the mass ranges of 2-90M Q (Figure 17) and 100-IOOOMq (Figure [HI). Note that above 



~ WOMq the effective temperature is roughly constant, T e g ~ 10 5 K, implying that the spectrum is 
independent of the mass distribution of the stars in this regime (Bromm et al. 2000). As is evident 
from these Figures (see also Tumlinson & Shull 2000), both the effective temperature and the 
ionizing power of metal-free (Pop III) stars are substantially larger than those of metal-rich (Pop 
I) stars. Metal-free stars with masses p> 20M@ emit between 10 47 and 10 48 H I and He I ionizing 
photons per second per solar mass of stars, where the lower value applies to stars of ~ 2OM and 
the upper value applies to stars of ^ 1OOM (see Tumlinson Sz Shull 2000 and Bromm et al. 2000 
for more details). These massive stars produce 10 -10 5 ionizing photons per stellar baryon over a 
lifetime of ~ 3 x 10 6 years [which is much shorter than the age of the universe, equation (pi]) in 
§2.1 1 . However, this powerful UV emission is suppressed as soon as the interstellar medium out of 



which new stars form is enriched by trace amounts of metals. Even though the collapsed fraction 
of baryons is small at the epoch of reionization, it is likely that most of the stars responsible for 
the reionization of the universe formed out of enriched gas. 

Will it be possible to infer the initial mass function (IMF) of the first stars from spectroscopic 
observations of the first galaxies? Figure [n] compares the observed spectrum from a Salpeter 
IMF (dN+/dM cx M~ 2 - 35 ) and a heavy IMF (with all stars more massive than IOOMq) for a 
galaxy at z s = 10. The latter case follows from the assumption that each of the dense clumps 
in the simulations described in the previous section ends up as a single star with no significant 
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Fig. 17. — Luminosity vs. effective temperature for zero-age main sequences stars in the mass range 
of 2-90M© (from Tumlinson & Shull 2000). The curves show Pop I (Z Q = 0.02, on the right) and 
Pop III stars (on the left) in the mass range 2-90 Mq. The diamonds mark decades in metallicity 
in the approach to Z = from 10~ 2 down to 10~ 5 at 2 Mq, down to 10 -10 at 15 Mq, and down 
to 10~ 13 at 90 Mq. The dashed line along the Pop III zero-age main sequence assumes pure H-He 
composition, while the solid line (on the left) marks the upper MS with Zq = 10 -10 for the M > 15 
M Q models. Squares mark the points corresponding to pre-enriched evolutionary models from El 
Eid et al. (1983) at 80 M & and from Castellani et al. (1983) at 25 M . 
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Fig. 18. — Same as Figure 17 but for very massive stars above 100M© (from Bromm et al. 2000). 
Left solid line: Pop III zero-age main sequence (ZAMS). Right solid line: Pop I ZAMS. In each 
case, stellar luminosity (in L & ) is plotted vs. effective temperature (in K). Diamond- shaped symbols: 
Stellar masses along the sequence, from IOOMq (bottom) to IOOOMq (top) in increments of IOOMq. 
The Pop III ZAMS is systematically shifted to higher effective temperature, with a value of ~ 10 5 
K which is approximately independent of mass. The luminosities, on the other hand, are almost 
identical in the two cases. 
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fragmentation or mass loss. The difference between the plotted spectra cannot be confused with 
simple reddening due to normal dust. Another distinguishing feature of the IMF is the expected 
flux in the hydrogen and helium recombination lines, such as Lya and He II 1640 A, from the 
interstellar medium surrounding these stars. We discuss this next. 




Fig. 19. — Comparison of the predicted flux from a Pop III star cluster at z s = 10 for a Salpeter 
IMF (Tumlinson & Shull 2000) and a massive IMF (Bromm et al. 2000). Plotted is the observed flux 
(in nJy per 10 6 M Q of stars) vs. observed wavelength (in [im) for a flat universe with = 0.7 and 
h = 0.65. Solid line: The case of a heavy IMF. Dotted line: The fiducial case of a standard Salpeter 
IMF. The cutoff below A & s = 1216 A (1 + z s ) = 1.34^m is due to Gunn-Peterson absorption. (The 
cutoff has been slightly smoothed here by the damping wing of the Lya line, with reionization 
assumed to occur at z = 7; see § |9,1.1| for details.) Clearly, for the same total stellar mass, the 
observable flux is larger by an order of magnitude for stars which are biased towards having masses 
£ 100M© . 



4-1-3. Emission of Recombination Lines from the First Galaxies 

The hard UV emission from a star cluster or a quasar at high redshift is likely reprocessed by 
the surrounding interstellar medium, producing very strong recombination lines of hydrogen and 
helium (Oh 1999; Tumlinson & Shull 2000; see also Baltz, Gnedin & Silk 1998). We define N ion 
to be the production rate per unit stellar mass of ionizing photons by the source. The emitted 
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luminosity L^ c per unit stellar mass in a particular recombination line is then estimated to be 

^U™c = PU™c^^^ion(l - Pcont)PU S ne j (53) 

where p^™ e is the probability that a recombination leads to the emission of a photon in the cor- 
responding line, v is the frequency of the line and Pcont an d P^ne are * ne esca P e probabilities for 
the ionizing photons and the line photons, respectively. It is natural to assume that the stellar 
cluster is surrounded by a finite H II region, and hence that Pcont i s close to zero (Wood & Loeb 
2000; Ricotti & Shull 2000). In addition, pf^ c is likely close to unity in the H II region, due to the 
lack of dust in the ambient metal-free gas. Although the emitted line photons may be scattered by 
neutral gas, they diffuse out to the observer and in the end survive if the gas is dust free. Thus, 
for simplicity, we adopt a value of unity for pf^ e (two-photon decay is generally negligible as a way 
of losing line photons in these environments). 

As a particular example we consider case B recombination which yields pf^ e of about 0.65 and 
0.47 for the Lya and He II 1640 A lines, respectively. These numbers correspond to an electron 
temperature of ~ 3 x 10 4 K and an electron density of ~ 10 2 — 10 3 cm -3 inside the H II region 
(Storey & Hummer 1995). For example, we consider the extreme and most favorable case of metal- 
free stars all of which are more massive than ~ IOOMq. In this case Lf^ e = 1.7 x 10 37 and 2.2 x 10 36 
erg s^Mq 1 for the recombination luminosities of Lya and He II 1640 A per stellar mass (Bromm 
et al. 2000). A cluster of 1O 6 M in such stars would then produce 4.4 and 0.6 xl0 9 L Q in the Lya 
and He II 1640 A lines. Comparably-high luminosities would be produced in other recombination 
lines at longer wavelengths, such as He II 4686 A and Ha (Oh 2000; Oh, Haiman, & Rees 2000). 

The rest-frame equivalent width of the above emission lines measured against the stellar con- 
tinuum of the embedded star cluster at the line wavelengths is given by 

/ rem \ 

W X = (^J , (54) 

where L\ is the spectral luminosity per unit wavelength of the stars at the line resonance. The 
extreme case of metal-free stars which are more massive than 100M Q yields a spectral luminosity 
per unit frequency L v = 2.7 'x 10 21 and 1.8 x 10 21 erg s _1 Hz _1 Mq 1 at the corresponding wavelengths 
(Bromm et al. 2000). Converting to L\, this yields rest-frame equivalent widths of W\ = 3100 A and 
1100 A for Lya and He II 1640 A, respectively. These extreme emission equivalent widths are more 
than an order of magnitude larger than the expectation for a normal cluster of hot metal-free stars 
with the same total mass and a Salpeter IMF under the same assumptions concerning the escape 
probabilities and recombination (Kudritzki et al. 2000). The equivalent widths are, of course, larger 
by a factor of (1 + z s ) in the observer frame. Extremely strong recombination lines, such as Lya 
and He II 1640 A, are therefore expected to be an additional spectral signature that is unique to 
very massive stars in the early universe. The strong recombination lines from the first luminous 
objects are potentially detectable with NGST (Oh, Haiman, & Rees 2000). 

High-redshift objects could also, in principle, be detected through their cooling radiation. 
However, a simple estimate of the radiated energy shows that it is very difficult to detect the 
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corresponding signal in practice. As it cools, the gas loses much of its gravitational binding energy, 
which is of order kBT v i T per baryon, with the virial temperature given by equation (26) in §2.3. 
Some fraction of this energy is then radiated as Lya photons. The typical galaxy halos around the 
reionization redshift have T v i r ~ 1 eV, and this must be compared to the nuclear energy output of 
7 MeV per baryon in stellar interiors. Clearly, for a star formation efficiency of ^ 1%, the stellar 
radiation is expected to be far more energetic than the cooling radiation. Both forms of energy 
should come out on a time-scale of order the dynamical time. Thus, even if the cooling radiation is 
concentrated in the Lya line, its detection is more promising for low redshift objects, while NGST 
will only be able to detect this radiation from the rare 4-cr halos (with masses ^> 10 1 1Mq) at z ~ 10 
(Haiman, Spaans, & Quataert 2000; Fardal et al. 2000). 



4.2. Black Hole Formation 

Quasars are more effective than stars in ionizing the intergalactic hydrogen because (i) their 
emission spectrum is harder, (ii) the radiative efficiency of accretion flows can be more than an 
order of magnitude higher than the radiative efficiency of a star, and (iii) quasars are brighter, and 
for a given density distribution in their host system, the escape fraction of their ionizing photons 
is higher than for stars. 

Thus, the history of reionization may have been greatly altered by the existence of massive 
black holes in the low-mass galaxies that populate the universe at high redshifts. For this reason, 
it is important to understand the formation of massive black holes (i.e., black holes with a mass far 
greater than a stellar mass). The problem of black hole formation is not a priori more complicated 
than the problem of star formation. Surprisingly, however, the amount of theoretical work on star 
formation far exceeds that on massive black hole formation. One of the reasons is that stars form 
routinely in our interstellar neighborhood where much data can be gathered, while black holes 
formed mainly in the distant past at great distances from our telescopes. As more information is 
gathered on the high-redshift universe, this state of affairs may begin to change. 

Here we adopt the view that massive black holes form out of gas and not through the dynamical 
evolution of dense stellar systems (see Rees 1984 for a review of the alternatives). To form a 
black hole inside a given dark matter halo, the baryons must cool. For most objects, this is only 
possible with atomic line cooling at virial temperatures T v ; r ;> 10 4 K and thus baryonic masses 
£ 10 7 Af Q [(l + z)/W] 3 / 2 . After losing their thermal pressure, the cold baryons collapse and form a 
thin disk on a dynamical time (Loeb & Rasio 1994). The basic question is then the following: what 
fraction of the cold baryons is able to sink to the very center of the potential well and form a massive 
black hole? Just as for star formation, the main barrier in this process is angular momentum. The 
centrifugal force opposes radial infall and keeps the gas in disks at a typical distance which is 
6-8 orders of magnitude larger than the Schwarzschild radius corresponding to the total gas mass. 
Eisenstein & Loeb (1995b) demonstrated that a small fraction of all objects have a sufficiently 
low angular momentum that the gas in them inevitably forms a compact semi-relativistic disk 
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that evolves to a black hole on a short viscous time-scale. These low-spin systems are born in 
special cosmological environments that exert unusually small tidal torques on them during their 
cosmological collapse. As long as the initial cooling time of the gas is short and its star formation 
efficiency is low, the gas forms the compact disk on a free-fall time. In most systems the baryons 
dominate gravity inside the scale length of the disk. Therefore, if the baryons in a low-spin system 
acquire a spin parameter which is only one sixth of the typical value, i.e., an initial rotation speed 
~ (16% x 0.05) x V c , then with angular momentum conservation they would reach rotational support 
at a radius r disk and circular velocity Vdisk such that Vdisk^disk ~ (16% x 0.05) V c fyi T , where r v ; r is 
the virial radius and V c the circular velocity of the halo. Using the relations: (GMh a i /V v i r ) ~ V c 2 , 
and [G^fc/fi^Mhaio/Tdisk] ~ ^disk' we & e t ^disk ~ 18V C . For T v i r ~ 10 4 K, the dark matter halo 
has a potential depth corresponding to a circular velocity of V c ~ 17 km s _1 , and the low-spin disk 
attains a characteristic rotation velocity of Vdisk ~ 300 km s~ x (sufficient to retain the gas against 
supernova-driven winds), a size <J 1 pc, and a viscous evolution time which is extremely short 
compared to the Hubble time. 

Low-spin dwarf galaxies populate the universe with a significant volume density at high red- 
shift; these systems are eventually incorporated into higher mass galaxies which form later. For 
example, a galactic bulge of ~ 1O 1O M0 in baryons forms out of ~ 10 3 building blocks of ~ 10 M 
each. In order to seed the growth of a quasar, it is sufficient that only one of these systems had 
formed a low-spin disk that produced a black hole progenitor. Note that if a low-spin object is 
embedded in an overdense region that eventually becomes a galactic bulge, then the black hole 
progenitor will sink to the center of the bulge by dynamical friction in less than a Hubble time 
(for a sufficiently high mass ^ 1O 6 M0 ; p. 428 of Binney & Tremaine) and seed quasar activ- 
ity. Based on the phase-space volume accessible to low-spin systems (oc j 3 ), we expect a fraction 
~ 6 -3 = 5 x 10 -3 of all the collapsed gas mass in the universe to be associated with low-spin 
disks (Eisenstein & Loeb 1995b). However, this is a conservative estimate. Additional angular 
momentum loss due to dynamical friction of gaseous clumps in dark matter halos (Navarro, Frenk, 
& White 1995) or bar instabilities in self-gravitating disks (Shlosman, Begelman, Sz Frank 1990) 
could only contribute to the black hole formation process. The popular paradigm that all galaxies 
harbor black holes at their center simply postulates that in all massive systems, a small fraction of 
the gas ends up as a black hole, but does not explain quantitatively why this fraction obtains its 
particular small value. The above scenario offers a possible physical context for this result. 

If the viscous evolution time is shorter than the cooling time and if the gas entropy is raised by 
viscous dissipation or shocks to a sufficiently high value, then the black hole formation process will 
go through the phase of a supermassive star (Shapiro & Teukolsky 1983, §17; see also Zel'dovich 
& Novikov 1971). The existence of angular momentum (Wagoner 1969) tends to stabilize the 
collapse against the instability which itself is due to general-relativistic corrections to the Newtonian 
potential (Shapiro & Teukolsky 1983, §17.4). However, shedding of mass and angular momentum 
along the equatorial plane eventually leads to collapse (Bisnovati-Kogan, Zel'dovich & Novikov 
1967; Loeb & Rasio 1994; Baumgarte & Shapiro 1999a). Since it is convectively unstable (Loeb & 
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Rasio 1994) and supported by radiation pressure, a supermassive star should radiate close to the 
Eddington limit (with modifications due to rotation; see Baumgarte & Shapiro 1999b) and generate 
a strong wind, especially if the gas is enriched with metals. The thermal+wind emission associated 
with the collapse of a supermassive star should be short-lived and could account for only a minority 
of all observed quasars. 

After the seed black hole forms, it is continually fed with gas during mergers. Mihos & 
Hernquist (1996) have demonstrated that mergers tend to deposit large quantities of gas at the 
centers of the merging galaxies, a process which could fuel a starburst or a quasar. If both of the 
merging galaxies contain black holes at their centers, dynamical friction will bring the black holes 
together. The final spiral-in of the black hole binary depends on the injection of new stars into orbits 
which allow them to extract angular momentum from the binary (Begelman, Blandford, & Rees 
1980). If the orbital radius of the binary shrinks to a sufficiently small value, gravitational radiation 
takes over and leads to coalescence of the two black holes. This will provide powerful sources for 



future gravitational wave detectors (such as the LISA project; see http: / /lisa.jpl. nasa.gov 



The fact that black holes are found in low-mass galaxies in the local universe implies that they 
are likely to exist also at high redshift. Local examples include the compact ellipticals M32 and 
NGC 4486B. In particular, van der Marel et al. (1997) infer a black hole mass of ~ 3.4 x 10 6 M Q in 
M32, which is a fraction ~ 8 x 10~ 3 of the stellar mass of the galaxy, ~ 4 x 10 8 M Q , for a central 
mass-to-light ratio of jy = 2. In NGC 4486B, Kormendy et al. (1997) infer a black hole mass of 
6 x 1O 8 M , which is a fraction ~ 9% of the stellar mass. 

Despite the poor current understanding of the black hole formation process, it is possible to 
formulate reasonable phenomenological prescriptions that fit the quasar luminosity function within 



the context of popular galaxy formation models. These prescription are described in § g.2.2 



5. Galaxy Properties 

5.1. Formation and Properties of Galactic Disks 

The formation of disk galaxies within hierarchical models of structure formation was first 
explored by Fall & Efstathiou (1980). More recently, the distribution of disk sizes was derived 
and compared to observations by Dalcanton, Spergel, & Summers (1997) and Mo, Mao, &; White 
(1998). Although these authors considered a number of detailed models, we adopt here the simple 
model of an exponential disk in a singular isothermal sphere halo. We consider a halo of mass M, 
virial radius r v i r , total energy E, and angular momentum J, for which the spin parameter is defined 
as 

A = J\E\ l l 2 G- l M~^l 2 . (55) 

The spin parameter simply expresses the halo angular momentum in a dimensionless form. The 
gas disk is assumed to collapse to a state of rotational support in the dark matter halo. If the disk 
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mass is a fraction of the halo mass and its angular momentum is a fraction of that of the 
halo, then the exponential scale radius of the disk is given by (Mo et al. 1998) 

R d = ^ (— ) Ar vir . (56) 
V2 \m d J 

The observed distribution of disk sizes suggests that the specific angular momentum of the 
disk is similar to that of the halo (e.g., Dalcanton et al. 1997; Mo et al. 1998), and so we assume 
that idl m d = 1- Although this result is implied by observed galactic disks, its origin in the disk 
formation process is still unclear. The formation of galactic disks has been investigated in a large 
number of numerical simulations (Navarro & Benz 1991; Evrard, Summers, & Davis 1994; Navarro, 
Frenk, & White 1995; Tissera, Lambas, & Abadi 1997; Navarro & Steinmetz 1997; Elizondo, et 
al. 1999). The overall conclusion is that the collapsing gas loses angular momentum to the dark 
matter halo during mergers, and the disks which form are much smaller than observed galactic 
disks. The most widely discussed solution for this problem is to prevent the gas from collapsing 
into a disk by injecting energy through supernova feedback (e.g. Eke, Efstathiou, & Wright 1999; 
Binney, Gerhard, & Silk 2001; Efstathiou 2000). However, some numerical simulations suggest that 
feedback may not adequately suppress the angular momentum losses (Navarro & Steinmetz 2000). 

With the assumption that jd/md = 1) the distribution of disk sizes is then determined by the 
Press-Schechter halo abundance and by the distribution of spin parameters [along with equation ( p4[ ) 
for r v ; r ]. The spin parameter distribution is approximately independent of mass, environment, 
and cosmological parameters, apparently a consequence of the scale-free properties of the early 
tidal torques between neighboring systems responsible for the spin of individual halos (Peebles 
1969; White 1984; Barnes & Efstathiou 1987; Heavens & Peacock 1988; Steinmetz & Bartelmann 
1995; Eisenstein & Loeb 1995a; Cole & Lacey 1996; Catelan & Theuns 1996). This distribution 
approximately follows a lognormal distribution in the vicinity of the peak, 



p(X)d\ = — -= exp 

0- A V27T 



In 2 (A/A) 
2al 



with A = 0.05 and a\ = 0.5 following Mo et al. (1998), who determined these values based on the 
N-body simulations of Warren et al. (1992). Although Mo et al. (1998) suggest a lower cutoff on A 
due to disk instability, it is unclear if halos with low A indeed cannot contain disks. If a dense bulge 
exists, it can prevent bar instabilities, or if a bar forms it may be weakened or destroyed when a 
bulge subsequently forms (Sellwood & Moore 1999). 



5.2. Phenomenological Prescription for Star Formation 



Schmidt (1959) put forth the hypothesis that the rate of star formation in a given region varies 
as a power of the gas density within that region. Thus, the star formation rate can be parameterized 
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as 

oc Pg , (58) 



where is the mass density of stars, and p g is the mass density of gas. Although Schmidt originally 
focused on different regions within our own Galaxy, this relation has since been used to interpret 
observations of the global star formation rates in different galaxies. 

One particular value of N is theoretically favored for self-gravitating disks (e.g., Larson 1992; 
Elmegreen 1994). The star formation rate can be written in the form 

< 59 > 

where e is an efficiency coefficient, and t c is a characteristic time for star formation. If t c is 

— 1/2 

proportional to the dynamical free-fall time, i.e., t c oc p g , then N = 1.5. However, observations 
yield estimates of surface densities E, not volume densities p. If the average gas scale height is 
roughly constant in different environments, then the same relation as equation (|58l ) should hold 
between the surface densities of stars and gas, with N = 1.5. 

Such a relation has, indeed, been observed to hold over a large range of physical conditions in 
galaxies. Synthetic models which include stellar evolution tracks and stellar atmosphere models are 
used to infer star formation rates using spectral observations of stellar populations. Star formation 
rates have been inferred in this way in the disks of normal spiral and irregular galaxies, most often 
using H a luminosities. Star formation also occurs in much denser environments in the nuclear 
regions of galaxies, where far- infrared luminosities are most useful for determining star formation 
rates. Thus, the relation between star formation and gas density has been measured over conditions 
ranging from the outskirts of normal disks to the central engines of infrared-luminous starburst 
galaxies. The result is a tight correlation in accordance with the Schmidt law, empirically given by 
(Kennicutt 1998) 



X SFR = (2.5 ± 0.7) x 10- 



v 1.4±0.15 



lM pc" 



M Q yr- 1 kpc- 2 , (60) 



where T*sfr and S 9as are the disk-averaged star formation rate and gas surface densities, respec- 
tively. This relation is observed to hold over almost five orders of magnitude in gas surface density. 



6. Radiative Feedback from the First Sources of Light 

6.1. Escape of Ionizing Radiation from Galaxies 

The intergalactic ionizing radiation field, a key ingredient in the development of reionization, 
is determined by the amount of ionizing radiation escaping from the host galaxies of stars and 
quasars. The value of the escape fraction as a function of redshift and galaxy mass remains a major 
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uncertainty in all current studies, and could affect the cumulative radiation intensity by orders of 
magnitude at any given redshift. Gas within halos is far denser than the typical density of the 
IGM, and in general each halo is itself embedded within an overdense region, so the transfer of the 
ionizing radiation must be followed in the densest regions in the universe. Numerical simulations 
of reionization are limited in their resolution of the densest regions and in the accuracy of their 
treatment of radiative transfer. 

The escape of ionizing radiation (hu > 13.6eV, A < 912 A) from the disks of present-day 
galaxies has been studied in recent years in the context of explaining the extensive diffuse ionized 
gas layers observed above the disk in the Milky Way (Reynolds et al. 1995) and other galaxies (e.g., 
Rand 1996; Hoopes, Walterbos, & Rand 1999). Theoretical models predict that of order 3-14% 
of the ionizing luminosity from O and B stars escapes the Milky Way disk (Dove & Shull 1994; 
Dove, Shull, & Ferrara 2000). A similar escape fraction of / esc = 6% was determined by Bland- 
Hawthorn & Maloney (1999) based on Ha measurements of the Magellanic Stream. From Hopkins 
Ultraviolet Telescope observations of four nearby starburst galaxies (Leitherer et al. 1995; Hurwitz, 
Jelinsky, &; Dixon 1997), the escape fraction was estimated to be in the range 3%< / CS c < 57%. If 
similar escape fractions characterize high-redshift galaxies, then stars could have provided a major 
fraction of the background radiation that reionized the IGM (e.g., Madau & Shull 1996; Madau 
1999). However, the escape fraction from high-redshift galaxies, which formed when the universe 
was much denser (p oc (1 + z) 3 ), may be significantly lower than that predicted by models meant to 
describe present-day galaxies. Current reionization calculations assume that galaxies are isotropic 
point sources of ionizing radiation and adopt escape fractions in the range 5% < / esc < 60% (see, 
e.g., Gnedin 2000a, Miralda-Escude et al. 2000). 

Clumping is known to have a significant effect on the penetration and escape of radiation from 
an inhomogeneous medium (e.g., Boisse 1990; Witt &; Gordon 1996, 2000; Neufeld 1991; Haiman 
& Spaans 1999; Bianchi et al. 2000). The inclusion of dumpiness introduces several unknown 
parameters into the calculation, such as the number and overdensity of the clumps, and the spatial 
correlation between the clumps and the ionizing sources. An additional complication may arise 
from hydrodynamic feedback, whereby part of the gas mass is expelled from the disk by stellar 
winds and supernovae (§0). 

Wood & Loeb (2000) used a three-dimensional radiation transfer code to calculate the steady- 
state escape fraction of ionizing photons from disk galaxies as a function of redshift and galaxy 
mass. The gaseous disks were assumed to be isothermal, with a sound speed c s ~ 10 km s , 
and radially exponential, with a scale-length based on the characteristic spin parameter and virial 
radius of their host halos. The corresponding temperature of ~ 10 4 K is typical for a gas which 
is continuously heated by photo-ionization from stars. The sources of radiation were taken to be 
either stars embedded in the disk, or a central quasar. For stellar sources, the predicted increase 
in the disk density with redshift resulted in a strong decline of the escape fraction with increasing 
redshift. The situation is different for a central quasar. Due to its higher luminosity and central 
location, the quasar tends to produce an ionization channel in the surrounding disk through which 
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much of its ionizing radiation escapes from the host. In a steady state, only recombinations in 
this ionization channel must be balanced by ionizations, while for stars there are many ionization 
channels produced by individual star-forming regions and the total recombination rate in these 
channels is very high. Escape fractions ^ 10% were achieved for stars at z ~ 10 only if ~ 90% of the 
gas was expelled from the disks or if dense clumps removed the gas from the vast majority 80%) 
of the disk volume (see Figure ^). This analysis applies only to halos with virial temperatures 
^ 10 4 K. Ricotti & Shull (2000) reached similar conclusions but for a quasi-spherical configuration 
of stars and gas. They demonstrated that the escape fraction is substantially higher in low-mass 
halos with a virial temperature <^ 10 4 K. However, the formation of stars in such halos depends on 
their uncertain ability to cool via the efficient production of molecular hydrogen (see § |3.3| ). 
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Fig. 20. — Escape fractions of stellar ionizing photons from a gaseous disk embedded within a 
1O 1O M halo which formed at z = 10 (from Wood & Loeb 2000). The curves show three different 
cases of dumpiness within the disk. The volume filling factor refers to either the ionizing emissivity, 
the gas clumps, or both, depending on the case. The escape fraction is substantial (^ 1%) only if 
the gas distribution is highly clumped. (Note: ISM is interstellar medium) 

The main uncertainty in the above predictions involves the distribution of the gas inside the 
host galaxy, as the gas is exposed to the radiation released by stars and the mechanical energy 
deposited by supernovae. Given the fundamental role played by the escape fraction, it is desirable 
to calibrate its value observationally. Recently, Steidel, Pettini, & Adelberger (2001) reported a 
preliminary detection of significant Lyman continuum flux in the composite spectrum of 29 Lyman 
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break galaxies (LBG) with redshifts in the range z = 3.40 ± 0.09. They co-added the spectra of 
these galaxies in order to be able to measure the low flux. Another difficulty in the measurement 
comes from the need to separate the Lyman-limit break caused by the interstellar medium from that 
already produced in the stellar atmospheres. After correcting for intergalactic absorption, Steidel 
et al. (2001) inferred a ratio between the emergent flux density at 1500A and 900A (rest frame) 
of 4.6 ± 1.0. Taking into account the fact that the stellar spectrum should already have an intrinsic 
Lyman discontinuity of a factor of ~ 3-5, but that only ~ 15-20% of the 1500A photons escape 
from typical LBGs without being absorbed by dust (Pettini et al. 1998a; Adelberger & Steidel 2000), 
the inferred 900A escape fraction is / csc ~ 10-20%. However, the observed blue spectrum suggests 
that these 29 particular LBGs may have a very low dust content, and the escape fraction in these 
galaxies may be 50% or higher (Haehnelt et al. 2001). Thus, although the galaxies in this sample 
were drawn from the bluest quartile of the LBG spectral energy distributions, the measurement 
implies that this quartile may itself dominate the hydrogen-ionizing background relative to quasars 
at z ~ 3. 



6.2. Propagation of Ionization Fronts in the IGM 

The radiation output from the first stars ionizes hydrogen in a growing volume, eventually 
encompassing almost the entire IGM within a single H II bubble. In the early stages of this 
process, each galaxy produces a distinct H II region, and only when the overall H II filling factor 
becomes significant do neighboring bubbles begin to overlap in large numbers, ushering in the 
"overlap phase" of reionization. Thus, the first goal of a model of reionization is to describe the 
initial stage, when each source produces an isolated expanding H II region. 

We assume a spherical ionized volume V, separated from the surrounding neutral gas by a 
sharp ionization front. In the case of a stellar ionizing spectrum, most ionizing photons are just 
above the hydrogen ionization threshold of 13.6 eV, where the absorption cross-section is high and 
a very thin layer of neutral hydrogen is sufficient to absorb all the ionizing photons. On the other 
hand, an ionizing source such as a quasar produces significant numbers of higher energy photons 
and results in a thicker transition region. 

In the absence of recombinations, each hydrogen atom in the IGM would only have to be 
ionized once, and the ionized proper volume V p would simply be determined by 

n H V p = N y , (61) 

where njj is the mean number density of hydrogen and N y is the total number of ionizing pho- 
tons produced by the source. However, the increased density of the IGM at high redshift implies 
that recombinations cannot be neglected. Indeed, in the case of a steady ionizing source (and 
neglecting the cosmological expansion), a steady-state volume would be reached corresponding to 
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the Stromgren sphere, with recombinations balancing ionizations: 

9 dN~ , , 

a B n 2 H V p = , (62) 

where the recombination rate depends on the square of the density and on the case B recombination 
coefficient olb = 2.6 x 1CP 13 cm 3 s" 1 for hydrogen at T = 10 4 K. The exact evolution for an 
expanding H II region, including a non-steady ionizing source, recombinations, and cosmological 
expansion, is given by (Shapiro & Giroux 1987) 

™{^-™V p ) = d J±- aB („l)v p . (63) 

In this equation, the mean density n# varies with time as l/a 3 (i). A critical feature of the physics 
of reionization is the dependence of recombination on the square of the density This means that if 
the IGM is not uniform, but instead the gas which is being ionized is mostly distributed in high- 
density clumps, then the recombination time is very short. This is often dealt with by introducing 
a volume-averaged clumping factor C (in general time-dependent), defined byQ 

2 \ 1-2 



C = {nij)/njj. (64) 

If the ionized volume is large compared to the typical scale of clumping, so that many clumps 
are averaged over, then C can be assumed to be approximately spatially uniform. In general, 
equation ( p3| ) can be solved by supplementing it with equation (Ell) and specifying C. Switching 
to the comoving volume V, the resulting equation is 

dV 1 dN^ C_ 

^ = <^r-°^ n ^' (65) 

where the present number density of hydrogen is 

n° H = 1.88 x 1(T 7 ( cm" 3 . (66) 
H ^0.022 J v ; 

This number density is lower than the total number density of baryons by a factor of ~ 0.76, 
corresponding to the primordial mass fraction of hydrogen. The solution for V(t) (generalized from 
Shapiro & Giroux 1987) around a source which turns on at t = t{ is 

ViQ-fJjrtgLS™*, (67) 
hi njj dt' 

where 

F(t',t) = -a B n° H J t t i ^ldt". (68) 



4 The recombination rate depends on the number density of electrons, and in using equation (^) we are neglecting 
the small contribution caused by partially or fully ionized helium. 
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At high redshift (equations (g) and (10) in §^]l]), and with the additional assumption of a constant 



C, the function F simplifies as follows. Defining 

f(t) = a(i)" 3/2 , (69) 

we derive 

F(t>, t) = C [f(t') - f(t)] = -0.262 [/(*') - f(jt)] , (70) 

where the last equality assumes C = 10 and our standard choice of cosmological parameters: 
J7 m = 0.3, = 0.7, and Qb = 0.045 (see the end of §|l]). Although this expression for F(t',t) is 
in general an accurate approximation at high redshift, in the particular case of the ACDM model 
(where Q m + $7 a = 1) we get the exact result by replacing equation (|69[) with 




m = \l- + ^r^- (7i) 

The size of the resulting H II region depends on the halo which produces it. Consider a halo 
of total mass M and baryon fraction f2f,/f2 m . To derive a rough estimate, we assume that baryons 
are incorporated into stars with an efficiency of /star = 10%, and that the escape fraction for the 
resulting ionizing radiation is also / esc = 10%. If the stellar IMF is similar to the one measured 
locally [Scalo 1998; equation @], then A 7 » 4000 ionizing photons are produced per baryon in 
stars (for a metallicity equal to 1/20 of the solar value). We define a parameter which gives the 
overall number of ionizations per baryon, 

A^ion = A 7 / s t ar / eS c • (72) 

If we neglect recombinations then we obtain the maximum comoving radius of the region which the 
halo of mass M can ionize, 

/ ■' N y V /3 _ f 3 N ion Q b MV /3 __, (N ion M x ! 



rmax_ ^<J -(^^"a^J ~ 675kpc [~WWM^J ' (73) 

for our standard set of parameters. Note that this radius is larger than the halo virial radius 
[equation (p4j)] by a factor of ~ 20, essentially independent of redshift and halo mass. The actual 
radius never reaches this size if the recombination time is shorter than the lifetime of the ionizing 
source. For an instantaneous starburst with the Scalo (1998) IMF [equation (|97|)], the production 
rate of ionizing photons can be approximated as (Haiman, personal communication) 

dN y a-lN~ f 1 iit<t s , 

— r 2 = x \ / t \- a (74) 

at a t s (r) otherwise, 

where A 7 = 4000, a = 4.5, and the most massive stars fade away with the characteristic time- 
scale t s = 3 x 10 6 yr. In Figure ^ we show the time evolution of the volume ionized by such a 
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source, with the volume shown in units of the maximum volume V mSuX which corresponds to r max in 
equation (^). We consider a source turning on at z = 10 (solid curves) or z = 15 (dashed curves), 
with three cases for each: no recombinations, C = 1, and C = 10, in order from top to bottom (Note 
that the result is independent of redshift in the case of no recombinations) . When recombinations 
are included, the volume rises and reaches close to Vm ax before dropping after the source turns off. 
At large t recombinations stop due to the dropping density, and the volume approaches a constant 
value (although V <C V max at large t if C = 10). 




t [year] 

Fig. 21. — Expanding H II region around an isolated ionizing source. The comoving ionized volume 
V is expressed in units of the maximum possible volume, V m3iX = 47rr^ iax /3 [with r max given in 
equation (f73|)1, and the time is measured after an instantaneous starburst which produces ionizing 
photons according to equation (|74|). We consider a source turning on at z = 10 (solid curves) or 
z = 15 (dashed curves), with three cases for each: no recombinations, C = 1, and C = 10, in order 
from top to bottom. The no-recombination curve is identical for the different source redshifts. 

We obtain a similar result for the size of the H II region around a galaxy if we consider a mini- 
quasar rather than stars. For the typical quasar spectrum (Elvis et al. 1994), if we assume a radiative 
efficiency of ~ 6% then roughly 11,000 ionizing photons are produced per baryon incorporated into 
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the black hole (Haiman, personal communication). The efficiency of incorporating the baryons in 
a galaxy into a central black hole is low 0.6% in the local universe, e.g. Magorrian et al. 1998; 
see also § |8.2.2 ), but the escape fraction for quasars is likely to be close to unity, i.e., an order of 
magnitude higher than for stars (see §|6.1|). Thus, for every baryon in galaxies, up to ~ 65 ionizing 
photons may be produced by a central black hole and ~ 40 by stars, although both of these numbers 
for -/V; on are highly uncertain. These numbers suggest that in either case the typical size of H II 
regions before reionization may be <J 1 Mpc or ~ 10 Mpc, depending on whether 10 8 Mq halos or 

1 1 



10 12 M(d halos dominate. 



6.3. Reionization of the IGM 

6.3.1. Hydrogen Reionization 

In this section we summarize recent progress, both analytic and numerical, made toward 
elucidating the basic physics of reionization and the way in which the characteristics of reionization 
depend on the nature of the ionizing sources and on other input parameters of cosmological models. 

The process of the reionization of hydrogen involves several distinct stages. The initial, "pre- 
overlap" stage (using the terminology of Gnedin 2000a) consists of individual ionizing sources 
turning on and ionizing their surroundings. The first galaxies form in the most massive halos at 
high redshift, and these halos are biased and are preferentially located in the highest-density regions. 



Thus the ionizing photons which escape from the galaxy itself (see §10]) must then make their way 
through the surrounding high-density regions, which are characterized by a high recombination 
rate. Once they emerge, the ionization fronts propagate more easily into the low-density voids, 
leaving behind pockets of neutral, high-density gas. During this period the IGM is a two-phase 
medium characterized by highly ionized regions separated from neutral regions by ionization fronts. 
Furthermore, the ionizing intensity is very inhomogeneous even within the ionized regions, with 
the intensity determined by the distance from the nearest source and by the ionizing luminosity of 
this source. 

The central, relatively rapid "overlap" phase of reionization begins when neighboring H II 
regions begin to overlap. Whenever two ionized bubbles are joined, each point inside their common 
boundary becomes exposed to ionizing photons from both sources. Therefore, the ionizing intensity 
inside H II regions rises rapidly, allowing those regions to expand into high-density gas which had 
previously recombined fast enough to remain neutral when the ionizing intensity had been low. 
Since each bubble coalescence accelerates the process of reionization, the overlap phase has the 
character of a phase transition and is expected to occur rapidly, over less than a Hubble time at 
the overlap redshift. By the end of this stage most regions in the IGM are able to see several 
unobscured sources, and therefore the ionizing intensity is much higher than before overlap and 
it is also much more homogeneous. An additional ingredient in the rapid overlap phase results 
from the fact that hierarchical structure formation models predict a galaxy formation rate that 
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rises rapidly with time at the relevant redshift range. This process leads to a state in which the 
low-density IGM has been highly ionized and ionizing radiation reaches everywhere except for gas 
located inside self-shielded, high-density clouds. This marks the end of the overlap phase, and this 
important landmark is most often referred to as the 'moment of reionization'. 

Some neutral gas does, however, remain in high-density structures which correspond to Lyman 
Limit systems and damped Lya systems seen in absorption at lower redshifts. The high-density 
regions are gradually ionized as galaxy formation proceeds, and the mean ionizing intensity also 
grows with time. The ionizing intensity continues to grow and to become more uniform as an 
increasing number of ionizing sources is visible to every point in the IGM. This "post-overlap" phase 
continues indefinitely, since collapsed objects retain neutral gas even in the present universe. The 
IGM does, however, reach another milestone (of limited significance) at z ~ 1.6, the breakthrough 
redshift (which is determined by the probability of intersecting Lyman limit systems; Madau, 
Haardt, & Rees 1999). Below this redshift, all ionizing sources are visible to each other, while 
above this redshift absorption by the Lya forest clouds implies that only sources in a small redshift 
range are visible to a typical point in the IGM. 

Semi-analytic models of the pre-overlap stage focus on the evolution of the H II filling factor, 
i.e., the fraction of the volume of the universe which is filled by H II regions. We distinguish between 
the naive filling factor F-& ii and the actual filling factor or porosity Qh ii- The naive filling factor 
equals the number density of bubbles times the average volume of each, and it may exceed unity 
since when bubbles begin to overlap the overlapping volume is counted multiple times. However, 
as explained below, in the case of reionization the linearity of the physics means that Fr n is a 
very good approximation to Qu n up to the end of the overlap phase of reionization. 

The model of individual H II regions presented in the previous section can be used to un- 
derstand the development of the total filling factor. Starting with equation (|65|), if we assume a 
common clumping factor C for all H II regions then we can sum each term of the equation over all 
bubbles in a given large volume of the universe, and then divide by this volume. Then V is replaced 
by the filling factor and iV 7 by the total number of ionizing photons produced up to some time t, 
per unit volume. The latter quantity equals the mean number of ionizing photons per baryon times 



the mean density of baryons n&. Following the arguments leading to equation (73), we find that if 
we include only stars then 

^ = Mongol , (75) 

n b 

where the collapse fraction F c0 \ is the fraction of all the baryons in the universe which are in 
galaxies, i.e., the fraction of gas which settles into halos and cools efficiently inside them. In 
writing equation (|75| ) we are assuming instantaneous production of photons, i.e., that the time- 
scale for the formation and evolution of the massive stars in a galaxy is short compared to the 
Hubble time at the formation redshift of the galaxy. In a model based on equation (jBlf), the 
near-equality between Fn n and Qu n results from the linearity of this equation. First, the total 
number of ionizations equals the total number of ionizing photons produced by stars, i.e., all 
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ionizing photons contribute regardless of the spatial distribution of sources; and second, the total 
recombination rate is proportional to the total ionized volume, regardless of its topology. Thus, 
even if two or more bubbles overlap the model remains an accurate approximation for Qh ii (at least 
until Qh ii becomes nearly equal to 1). Note, however, that there still are a number of important 
simplifications in the model, including the assumption of a homogeneous (though possibly time- 
dependent) clumping factor, and the neglect of feedback whereby the formation of one galaxy may 
suppress further galaxy formation in neighboring regions. These complications are discussed in 
detail below and in § |6.5| and §0. 

Under these assumptions we convert equation (|65|), which describes individual H II regions, to 
an equation which statistically describes the transition from a neutral universe to a fully ionized 
one (compare Madau et al. 1999 and Haiman & Loeb 1997): 

dQn ii iVjon dF co \ C _ 

~dT = 076 ~ir - a ^ n »Q* n , (76) 

where we assumed a primordial mass fraction of hydrogen of 0.76. The solution (in analogy with 
equation (p7\)) is 

where F(t',t) is determined by equations (|68[)-([7l|). 



A simple estimate of the collapse fraction at high redshift is the mass fraction (given by 



equation (31) in the Press-Schechter model) in halos above the cooling threshold, which is the 
minimum mass of halos in which gas can cool efficiently. Assuming that only atomic cooling is 
effective during the redshift range of reionization (§ |3.3|) , the minimum mass corresponds roughly 



to a halo of virial temperature T v i r = 10 4 K, which can be converted to a mass using equation (26). 
With this prescription we derive (for N- lon = 40) the reionization history shown in Figure ^ for 
the case of a constant clumping factor C. The solid curves show Qh ii as a function of redshift for 
a clumping factor C = (no recombinations), C = 1, C = 10, and C = 30, in order from left to 
right. Note that if C ~ 1 then recombinations are unimportant, but if C ;> 10 then recombinations 
significantly delay the reionization redshift (for a fixed star-formation history). The dashed curve 
shows the collapse fraction F co \ in this model. For comparison, the vertical dotted line shows the 
z = 5.8 observational lower limit (Fan et al. 2000) on the reionization redshift. 

Clearly, star-forming galaxies in CDM hierarchical models are capable of ionizing the universe 
at z ~ 6-15 with reasonable parameter choices. This has been shown by a number of theo- 
retical, semi-analytic calculations (Fukugita & Kawasaki 1994; Shapiro, Giroux, & Babul 1994; 
Kamionkowski, Spergel, & Sugiyama 1994; Tegmark, Silk, & Blanchard 1994; Haiman & Loeb 
1997; Valageas & Silk 1999; Chiu & Ostriker 2000; Ciardi et al. 2000) as well as numerical simula- 
tions (Cen & Ostriker 1993; Gnedin & Ostriker 1997; Gnedin 2000a). Similarly, if a small fraction 
(^ 1%) °f tli e gas in each galaxy accretes onto a central black hole, then the resulting mini-quasars 
are also able to reionize the universe, as has also been shown using semi-analytic models (Fukugita 
& Kawasaki 1994; Haiman & Loeb 1998; Valageas & Silk 1999). Note that the prescription whereby 
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Fig. 22. — Semi-analytic calculation of the reionization of the IGM (for N[ on = 40), showing the 
redshift evolution of the filling factor Qh ii- Solid curves show Qn n for a clumping factor C = 
(no recombinations), C = 1, C = 10, and C = 30, in order from left to right. The dashed curve 
shows the collapse fraction F co \, and the vertical dotted line shows the z = 5.8 observational lower 
limit (Fan et al. 2000) on the reionization redshift. 
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a constant fraction of the galactic mass accretes onto a central black hole is based on local obser- 
vations (see § g.2.2| ) which indicate that z = galaxies harbor central black holes of mass equal 



to ~ 0.2-0.6% of their bulge mass. Although the bulge constitutes only a fraction of the total 
baryonic mass of each galaxy, the higher gas-to-stellar mass ratio in high-redshift galaxies, as well 
as their high merger rates compared to their low-redshift counterparts, suggest that a fraction of 
a percent of the total gas mass in high-redshift galaxies may have contributed to the formation of 
quasar black holes. 

Although many models yield a reionization redshift around 7-12, the exact value depends on 
a number of uncertain parameters affecting both the source term and the recombination term in 
equation (|76|). The source parameters include the formation efficiency of stars and quasars and the 
escape fraction of ionizing photons produced by these sources. The formation efficiency of low-mass 
galaxies may also be reduced by feedback from galactic outflows. These parameters affecting the 
sources are discussed elsewhere in this review (see § |5.2j , |3.2,2| , |P| , and [?]) . Even when the clumping 
is inhomogeneous, the recombination term in equation (^) is generally valid if C is defined as in 
equation (|64]), where we take a global volume average of the square of the density inside ionized 
regions (since neutral regions do not contribute to the recombination rate). The resulting mean 
clumping factor depends on the density and clustering of sources, and on the distribution and 
topology of density fluctuations in the IGM. Furthermore, the source halos should tend to form in 
overdense regions, and the clumping factor is affected by this cross-correlation between the sources 
and the IGM density. 

Valageas & Silk (1999) and Chiu & Ostriker (2000) calculated the clumping factor semi- 
analytically by averaging over the IGM on the one hand and virialized halos on the other hand, 
with the average weighed according to the gas fraction in halos. The semi- analytic methods used 
in these two detailed calculations of reionization have different advantages: Valageas &: Silk (1999) 
included a model for clumping and absorption by Lya clouds, but Chiu & Ostriker (2000) used a 
generally more realistic two-phase model with separate ionized and neutral regions. Miralda-Escude, 
Haehnelt, & Rees (2000) went further in their modeling of the clumping factor by attempting to 
account for the geometry of ionized regions. They presented a simple model for the distribution of 
density fluctuations, and more generally they discussed the implications of inhomogeneous clumping 
during reionization. They noted that as ionized regions grow, they more easily extend into low- 
density regions, and they tend to leave behind high-density concentrations, with these neutral 
islands being ionized only at a later time. They therefore argued that, since at high-redshift the 
collapse fraction is low, most of the high-density regions, which would dominate the clumping factor 
if they were ionized, will in fact remain neutral and occupy only a tiny fraction of the total volume. 
Thus, the development of reionization through the end of the overlap phase should occur almost 
exclusively in the low-density IGM, and the effective clumping factor during this time should be ~ 1, 



making recombinations relatively unimportant (see Figure 22). Only in the post-reionization phase, 
Miralda-Escude et al. (2000) argued, do the high density clouds and filaments become gradually 
ionized as the mean ionizing intensity further increases. 
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The complexity of the process of reionization is illustrated by the recent numerical simulation 
by Gnedin (2000a) of stellar reionization (in ACDM with VL m = 0.3). This simulation uses a 
formulation of radiative transfer which relies on several rough approximations; although it does 
not include the effect of shadowing behind optically-thick clumps, it does include for each point 
in the IGM the effects of an estimated local optical depth around that point, plus a local optical 
depth around each ionizing source. This simulation helps to understand the advantages of the 
various theoretical approaches, while pointing to the complications which are not included in the 



simple models. Figures 23 and [2j, taken from Figure 3 in Gnedin (2000a), show the state of the 
simulated universe just before and just after the overlap phase, respectively. They show a thin (15 
h~ l comoving kpc) slice through the box, which is 4 h~ l Mpc on a side. The simulation achieves a 
spatial resolution of l/i -1 kpc, and uses 128 3 each of dark matter particles and baryonic particles 
(with each baryonic particle having a mass of 5 x 10 5 Mq). The figures show the redshift evolution 
of the ionizing intensity averaged over the entire volume, J21 (upper right panel), and, visually, the 
spatial distribution of three quantities: the neutral hydrogen fraction (upper left panel), the gas 
density (lower left panel), and the gas temperature (lower right panel). Note the obvious features 
around the edges, resulting from the periodic boundary conditions assumed in the simulation (e.g., 
the left and right edges match identically). Also note that the intensity J21 is defined as the 
radiation intensity at the Lyman limit, expressed in units of 10 -21 erg cm -2 s _1 sr" ^Hz -1 . For a 
given source emission, the intensity inside H II regions depends on absorption and radiative transfer 
through the IGM (e.g., Haardt & Madau 1996; Abel & Haehnelt 1999) 

Figure ^3] shows the two-phase IGM at z = 7.7, with ionized bubbles emanating from many 
independent sources, although there is one main concentration (located at the right edge of the 
image, vertically near the center; note the periodic boundary conditions). The bubbles are shown 
expanding into low density regions and beginning to overlap at the center of the image. The 
topology of ionized regions is clearly complex: While the ionized regions are analogous to islands 
in an ocean of neutral hydrogen, the islands themselves contain small lakes of dense neutral gas. 
One aspect which has not been included in theoretical models of clumping is clear from the figure. 
The sources themselves are located in the highest density regions (these being the sites where 
the earliest galaxies form) and must therefore ionize the gas in their immediate vicinity before 
the radiation can escape into the low density IGM. For this reason, the effective clumping factor 
is of order 100 in the simulation and also, by the overlap redshift, roughly ten ionizing photons 
have been produced per baryon. As emphasized by Gnedin (2000a), some of these numbers are 
resolution dependent, since the clumping factor accounts only for absorption by gas at the highest 
resolvable density. A higher-resolution simulation would have higher-density gas clumps and - 
depending on the geometry of those clumps — a higher or possibly lower clumping factor than the 



low-resolution simulation. Figure 24 shows that by z = 6.7 the low density regions have all become 
highly ionized along with a rapid increase in the ionizing intensity. The only neutral islands left 
are the highest density regions which are not near the sources (compare the two panels on the 
left). However, we emphasize that the quantitative results of this simulation must be considered 
preliminary, since the effects of increased resolution and a more accurate treatment of radiative 
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transfer are yet to be explored. Methods are being developed for incorporating a more complete 
treatment of radiative transfer into three dimensional cosmological simulations (e.g., Abel, Norman, 
& Madau 1999; Razoumov & Scott 1999). 




Fig. 23. — Visualization at z = 7.7 of a numerical simulation of reionization, adopted from Figure 3c 
of Gnedin (2000a). The panels display the logarithm of the neutral hydrogen fraction (upper left), 
the gas density in units of the cosmological mean (lower left), and the gas temperature in Kelvin 
(lower right). These panels show a two-dimensional slice of the simulation (not a two-dimensional 
projection). Also shown is the redshift evolution of the logarithm of the ionizing intensity averaged 
over the entire simulation volume (upper right). Note the periodic boundary conditions. 

Gnedin, Ferrara, & Zweibel (2000) investigated an additional effect of reionization. They 
showed that the Biermann battery in cosmological ionization fronts inevitably generates coherent 
magnetic fields of an amplitude ~ 10 -19 Gauss. These fields form as a result of the breakout of 
the ionization fronts from galaxies and their propagation through the H I filaments in the IGM. 
Although the fields are too small to directly affect galaxy formation, they could be the seeds for 
the magnetic fields observed in galaxies and X-ray clusters today. 

If quasars contribute substantially to the ionizing intensity during reionization then several 
aspects of reionization are modified compared to the case of pure stellar reionization. First, the 
ionizing radiation emanates from a single, bright point-source inside each host galaxy, and can 
establish an escape route (H II funnel) more easily than in the case of stars which are smoothly 
distributed throughout the galaxy (§|6.1|). Second, the hard photons produced by a quasar penetrate 
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Fig. 24. — Visualization at z = 6.7 of a numerical simulation of reionization, adopted from Figure 3c 
of Gnedin (2000a). The panels display the logarithm of the neutral hydrogen fraction (upper left), 
the gas density in units of the cosmological mean (lower left), and the gas temperature in Kelvin 
(lower right). These panels show a two-dimensional slice of the simulation (not a two-dimensional 
projection). Also shown is the redshift evolution of the logarithm of the ionizing intensity averaged 
over the entire simulation volume (upper right). Note the periodic boundary conditions. 
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deeper into the surrounding neutral gas, yielding a thicker ionization front. Finally, the quasar X- 
rays catalyze the formation of H2 molecules and allow stars to keep forming in very small halos 

(8ED- 

Oh (2000) showed that star-forming regions may also produce significant X-rays at high red- 
shift. The emission is due to inverse Compton scattering of CMB photons off relativistic electrons 
in the ejecta, as well as thermal emission by the hot supernova remnant. The spectrum expected 
from this process is even harder than for typical quasars, and the hard photons photoionize the 
IGM efficiently by repeated secondary ionizations. The radiation, characterized by roughly equal 
energy per logarithmic frequency interval, would produce a uniform ionizing intensity and lead to 
gradual ionization and heating of the entire IGM. Thus, if this source of emission is indeed effective 
at high redshift, it may have a crucial impact in changing the topology of reionization. Even if stars 
dominate the emission, the hardness of the ionizing spectrum depends on the initial mass function. 
At high redshift it may be biased toward massive, efficiently ionizing stars (see § [4.1.1| ), but this 
remains very much uncertain. 

Semi-analytic as well as numerical models of reionization depend on an extrapolation of hierar- 
chical models to higher redshifts and lower-mass halos than the regime where the models have been 
compared to observations. These models have the advantage that they are based on the current 
CDM paradigm which is supported by a variety of observations of large-scale structure, galaxy clus- 
tering, and the CMB. The disadvantage is that the properties of high-redshift galaxies are derived 
from those of their host halos by prescriptions which are based on low-redshift observations, and 
these prescriptions will only be tested once abundant data is available on galaxies which formed 
during the reionization era. An alternative approach to analyzing the possible ionizing sources 
which brought about reionization is to extrapolate from the observed populations of galaxies and 
quasars at currently accessible redshifts. This has been attempted, e.g., by Madau et al. (1999) 
and Miralda-Escude et al. (2000). The general conclusion is that a high-redshift source population 
similar to the one observed at z = 3-4 would produce roughly the needed ionizing intensity for 
reionization. A precise conclusion, however, remains elusive because of the same kinds of uncer- 
tainties as those found in the models based on CDM: The typical escape fraction, and the faint 
end of the luminosity function, are both not well determined even at z = 3-4, and in addition the 
clumping factor at high redshift must be known in order to determine the importance of recombi- 
nations. Future direct observations of the source population at redshifts approaching reionization 
may help resolve some of these questions. 



6.3.2. Helium Reionization 

The sources that reionized hydrogen very likely caused the single reionization of helium from 
He I to He II. Neutral helium is ionized by photons of 24.6 eV or higher energy, and its recombination 
rate is roughly equal to that of hydrogen. On the other hand, the ionization threshold of He II 
is 54.4 eV, and fully ionized helium recombines ^ 5 times faster than hydrogen. This means that 
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for both quasars and galaxies, the reionization of He II should occur later than the reionization of 
hydrogen, even though the number of helium atoms is smaller than hydrogen by a factor of 13. 
The lower redshift of He II reionization makes it more accessible to observations and allows it to 
serve in some ways as an observational preview of hydrogen reionization. 

The Lya absorption by intergalactic He II (at wavelength 304A) has been observed in four 
quasars at redshifts 2.4 < z < 3.2 (Jakobsen et al. 1994; Davidsen et al. 1996; Hogan et al. 1997; 
Reimers et al. 1997; Anderson et al. 1999; Heap et al. 2000). The results are consistent among the 
different quasars, and we illustrate them here with one particular spectrum. In Figure [2^, adopted 
from Figure 4 of Heap et al. (2000), we show a portion of the spectrum of the z = 3.286 quasar 
Q 0302 — 003, obtained with the Space Telescope Imaging Spectrograph on-board the Hubble Space 
Telescope. The observed spectrum (solid line) is compared to a simulated spectrum (gray shading) 
based on the H I Lya forest observed in the same quasar. In deriving the simulated spectrum, Heap 
et al. assumed a ratio of He II to H I column densities of 100, and pure turbulent line broadening. 
The wavelength range shown in the figure corresponds to He II Lya in the redshift range 2.8-3.3. 

The observed flux shows a clear break shortward of the quasar emission line at an observed 
A = 1300 A. Relatively near the quasar, at A =1285-1300A, a shelf of relatively high transmission is 
likely evidence of the 'proximity effect', in which the emission from the quasar itself creates a highly 
ionized local region with a reduced abundance of absorbing ions. In the region at A =1240-1280 
A (z = 3.08-3.21), on the other hand, the very low flux level implies an average optical depth of 
r ~ 4.5-5 due to He II. Another large region with average r ~ 4, a region spanning ~ 100 comoving 
Mpc along the line of sight, is evident at A =1180-1210A (z = 2.88-2.98). The strong continuous 
absorption in these large regions, and the lack of correlation with the observed H I Lya forest, is 
evidence for a He II Gunn-Peterson absorption trough due to the diffuse IGM. It also suggests a 
rather soft UV background with a significant stellar contribution, i.e., a background that ionizes 
the diffuse hydrogen much more thoroughly than He II. Significant emission is observed in between 
the two regions of constant high absorption. A small region around 1216A is contaminated by 
geo-coronal Lya, but the emission at 1230-1235A apparently corresponds to a real, distinct gap 
in the He II abundance, which could be caused by a local source photo-ionizing a region of radius 
~ 10 comoving Mpc. The region at A =1150-1175A (z = 2.78-2.86) shows a much higher overall 
transmission level than the regions at slightly higher redshift. Heap et al. measure an average 
r = 1.9 in this region, and note that the significant correlation of the observed spectrum with 
the simulated one suggests that much of the absorption is due to a He II Lya forest while the 
low-density IGM provides a relatively low opacity in this region. The authors conclude that the 
observed data suggest a sharp opacity break occurring between z = 3.0 and 2.9, accompanied by a 
hardening of the UV ionizing background. However, even the relatively high opacity at z ;> 3 only 
requires ~ 0.1% of helium atoms not to be fully ionized, in a region at the mean baryon density. 
Thus, the overlap phase of full helium reionization may have occurred significantly earlier, with the 
ionizing intensity already fairly uniform but still increasing with time at z ~ 3. 

The properties of helium reionization have been investigated numerically by a number of au- 
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thors. Zheng & Davidsen (1995) modeled the He II proximity effect, and a number of authors 
(Miralda-Escude et al. 1996; Croft et al. 1997; Zhang et al. 1998) used numerical simulations to 
show that the observations generally agree with cold dark matter models. They also found that 
helium absorption particularly tests the properties of under-dense voids which produce much of 
the He II opacity but little opacity in H I. According to the semi-analytic model of inhomogeneous 
reionization of Miralda-Escude, Haehnelt, & Rees (2000; see also § |6.3.1 ), the total emissivity of ob- 



served quasars at redshift 3 suffices to completely reionize helium before z = 3. They find that the 
observations at z ~ 3 can be reproduced if a population of low-luminosity sources, perhaps galaxies, 
has ionized the low-density IGM up to an overdensity of around 12 relative to the cosmological 
mean, with luminous quasars creating the observed gaps of transmitted flux. 

The conclusion that an evolution of the ionization state of helium has been observed is also 
strengthened by several indirect lines of evidence. Songaila & Cowie (1996) and Songaila (1998) 
found a rapid increase in the Si IV/C IV ratio with decreasing redshift at z = 3, for intermediate 
column density hydrogen Lya absorption lines. They interpreted this evolution as a sudden hard- 
ening below z = 3 of the spectrum of the ionizing background. Boksenberg et al. (1998) also found 
an increase in the Si IV/C IV ratio, but their data implied a much more gradual increase from 
z = 3.8 to z = 2.2. 

The full reionization of helium due to a hard ionizing spectrum should also heat the IGM to 
20,000 K or higher, while the IGM can only reach ~10,000 K during a reionization of hydrogen alone 
(although a temperature of ~15,000 K may be reached due to Compton heating by the hard X-ray 
background: Madau & Efstathiou 1999). This increase in temperature can serve as an observational 
probe of helium reionization, and it should also increase the suppression of dwarf galaxy formation 
(§ |6.5[ ). The temperature of the IGM can be measured by searching for the smallest line- widths 
among hydrogen Lya absorption lines (Schaye et al. 1999). In general, bulk velocity gradients 
contribute to the line width on top of thermal velocities, but a lower bound on the width is set 
by thermal broadening, and the narrowest lines can be used to measure the temperature. Several 
different measurements (Ricotti et al. 2000; Schaye et al. 2000; Bryan & Machacek 2000; McDonald 
et al. 2000) have found a nearly isothermal IGM at a temperature of ~20,000 K at z = 3, higher than 
expected in ionization equilibrium and suggestive of photo-heating due to ongoing reionization of 
helium. However, the measurement errors remain too large for a firm conclusion about the redshift 
evolution of the IGM temperature or its equation of state. 

Clearly, the reionization of helium is already a rich phenomenological subject. Our knowledge 
will benefit from measurements of increasing accuracy, made toward many more lines of sight, 
and extended to higher redshift. New ways to probe helium will also be useful. For example, 
Miralda-Escude (2000) has suggested that continuum He II absorption in soft X-rays can be used 
to determine the He II fraction along the line of sight, although the measurement requires an 
accurate subtraction of the Galactic contribution to the absorption, based on the Galactic H I 
column density as determined by 21 cm maps. 
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6.4. Photo-evaporation of Gaseous Halos After Reionization 



The end of the reionization phase transition resulted in the emergence of an intense UV back- 
ground that filled the universe and heated the IGM to temperatures of ~ 1-2 x 10 4 K (see the 
previous section). After ionizing the rarefied IGM in the voids and filaments on large scales, the 
cosmic UV background penetrated the denser regions associated with the virialized gaseous halos 
of the first generation of objects. A major fraction of the collapsed gas had been incorporated by 
that time into halos with a virial temperature <^ 10 4 K, where the lack of atomic cooling prevented 
the formation of galactic disks and stars or quasars. Photoionization heating by the cosmic UV 
background could then evaporate much of this gas back into the IGM. The photo-evaporating halos, 
as well as those halos which did retain their gas, may have had a number of important consequences 
just after reionization as well as at lower redshifts. 

In this section we focus on the process by which gas that had already settled into virialized halos 
by the time of reionization was evaporated back into the IGM due to the cosmic UV background. 
This process was investigated by Barkana & Loeb (1999) using semi-analytic methods and idealized 
numerical calculations. They first considered an isolated spherical, centrally-concentrated dark 
matter halo containing gas. Since most of the photo-evaporation occurs at the end of overlap, 
when the ionizing intensity builds up almost instantaneously, a sudden illumination by an external 
ionizing background may be assumed. Self-shielding of the gas implies that the halo interior sees 
a reduced intensity and a harder spectrum, since the outer gas layers preferentially block photons 
with energies just above the Lyman limit. It is useful to parameterize the external radiation field 
by a specific intensity per unit frequency, v, 



where vl is the Lyman limit frequency, and J21 is the intensity at vl expressed in units of 
10~ 21 erg cm~ 2 s _1 sr _1 Hz _1 . The intensity is normalized to an expected post-reionization value 
of around unity for the ratio of ionizing photon density to the baryon density. Different power laws 
can be used to represent either quasar spectra (a ~ 1.8) or stellar spectra (a ~ 5). 

Once the gas is heated throughout the halo, some fraction of it acquires a sufficiently high 
temperature that it becomes unbound. This gas expands due to the resulting pressure gradi- 
ent and eventually evaporates back into the IGM. The pressure gradient force (per unit volume) 
A;bV '(T 'p I pm p ) competes with the gravitational force of pGM/r 2 . Due to the density gradient, the 
ratio between the pressure force and the gravitational force is roughly equal to the ratio between 
the thermal energy ~ k-p,T and the gravitational binding energy ~ pm p GM jr (which 
at the virial radius r v ; r ) per particle. Thus, if the kinetic energy exceeds the potential energy (or 
roughly if T > T v ; r ), the repulsive pressure gradient force exceeds the attractive gravitational force 
and expels the gas on a dynamical time (or faster for halos with T S> T v j r ). 

The left panel of Figure ^ (adopted from Figure 3 of Barkana & Loeb 1999) shows the fraction 
of gas within the virial radius which becomes unbound after reionization, as a function of the total 
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halo circular velocity, with halo masses at z = 8 indicated at the top. The two pairs of curves 
correspond to spectral index a = 5 (solid) or a = 1.8 (dashed). In each pair, a calculation which 
assumes an optically-thin halo leads to the upper curve, but including radiative transfer and self- 
shielding modifies the result to the one shown by the lower curve. In each case self-shielding lowers 
the unbound fraction, but it mostly affects only a neutral core containing ~ 30% of the gas. Since 
high energy photons above the Lyman limit penetrate deep into the halo and heat the gas efficiently, 
a flattening of the spectral slope from a = 5 to a = 1.8 raises the unbound gas fraction. This Figure 
is essentially independent of redshift if plotted in terms of circular velocity, but the conversion to 
a corresponding mass does vary with redshift. The characteristic circular velocity where most of 
the gas is lost is ~ 10-15 km s _1 , but clearly the effect of photo-evaporation is gradual, going from 
total gas removal down to no effect over a range of a factor of ~ 100 in halo mass. 

Given the values of the unbound gas fraction in halos of different masses, the Press-Schechter 
mass function (§ |2.4| ) can be used to calculate the total fraction of the IGM which goes through 
the process of accreting onto a halo and then being recycled into the IGM at reionization. The 
low-mass cutoff in this sum over halos is given by the lowest mass halo in which gas has assembled 
by the reionization redshift. This mass can be estimated by the linear Jeans mass Mj in equation 
@ in {§□]. The Jeans mass does not in general precisely equal the limiting mass for accretion (see 
the discussion in the next section). Indeed, at a given redshift some gas can continue to fall into 
halos of lower mass than the Jeans mass at that redshift. On the other hand, the larger Jeans mass 
at higher redshifts means that a time-averaged Jeans mass may be more appropriate, as indicated 
by the filtering mass. In practice, the Jeans mass is sufficiently accurate since at z ~ 10-20 it agrees 
well with the values found in the numerical spherical collapse calculations of Haiman, Thoul, &: 
Loeb (1996). 

The right panel of Figure |26] (adopted from Figure 7 of Barkana & Loeb 1999) shows the total 
fraction of gas in the universe which evaporates from halos at reionization, versus the reionization 
redshift. The solid line assumes a spectral index a = 1.8, and the dotted line assumes a = 5, 
showing that the result is insensitive to the spectrum. Even at high redshift, the amount of 
gas which participates in photo-evaporation is significant, which suggests a number of possible 
implications as discussed below. The gas fraction shown in the figure represents most (~ 60-80% 
depending on the redshift) of the collapsed fraction before reionization, although some gas does 
remain in more massive halos. 

The photo-evaporation of gas out of large numbers of halos may have interesting implications. 
First, gas which falls into halos and is expelled at reionization attains a different entropy than if 
it had stayed in the low-density IGM. The resulting overall reduction in the entropy is expected 
to be small — the same as would be produced by reducing the temperature of the entire IGM by 
a factor of ~ 1.5 — but localized effects near photo-evaporating halos may be more significant. 
Furthermore, the resulting ~ 20 km s -1 outflows induce small-scale fluctuations in peculiar veloc- 
ity and temperature. These outflows are usually well below the resolution limit of most numerical 
simulations, but some outflows were resolved in the simulation of Bryan et al. (1998). The evap- 
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Fig. 25. — Lya absorption by intergalactic He II. This spectrum of the z = 3.286 quasar Q 0302 — 
003, adopted from Figure 4 of Heap et al. (2000), was obtained using the Space Telescope Imaging 
Spectrograph. The observed spectrum (solid line) is compared to a simulated spectrum (gray 
shading) based on the H I Lya forest observed in the same quasar. In deriving the simulated 
spectrum, Heap et al. assumed a ratio of He II to H I column densities of 100, and pure turbulent 
line broadening. 
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Fig. 26. — Effect of photo-evaporation on individual halos and on the overall halo population. The 
left panel shows the unbound gas fraction (within the virial radius) versus total halo circular velocity 
or mass, adopted from Figure 3 of Barkana & Loeb (1999). The two pairs of curves correspond to 
spectral index a = 5 (solid) or a = 1.8 (dashed), in each case at z = 8. In each pair, assuming 
an optically-thin halo leads to the upper curve, while the lower curve shows the result of including 
radiative transfer and self shielding. The right panel shows the total fraction of gas in the universe 
which evaporates from halos at reionization, versus the reionization redshift, adopted from Figure 7 
of Barkana & Loeb (1999). The solid line assumes a spectral index a = 1.8, and the dotted line 
assumes a = 5. 
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orating halos may consume a significant number of ionizing photons in the post-overlap stage of 
reionization (e.g., Haiman, Abel, & Madau 2000), but a definitive determination requires detailed 
simulations which include the three-dimensional geometry of source halos and sink halos. 

Although gas is quickly expelled out of the smallest halos, photo-evaporation occurs more 
gradually in larger halos which retain some of their gas. These surviving halos initially expand but 
they continue to accrete dark matter and to merge with other halos. These evaporating gas halos 
could contribute to the high column density end of the Lya forest (Bond, Szalay, & Silk 1988). Abel 
& Mo (1998) suggested that, based on the expected number of surviving halos, a large fraction of the 
Lyman limit systems at z ~ 3 may correspond to mini-halos that survived reionization. Surviving 
halos may even have identifiable remnants in the present universe, as discussed in §[0|. These ideas 
thus offer the possibility that a population of halos which originally formed prior to reionization 
may correspond almost directly to several populations that are observed much later in the history of 
the universe. However, the detailed dynamics of photo-evaporating halos are complex, and detailed 
simulations are required to confirm these ideas. Photo-evaporation of a gas cloud has been followed 
in a two dimensional simulation with radiative transfer, by Shapiro & Raga (2000). They found 
that an evaporating halo would indeed appear in absorption as a damped Lya system initially, and 
as a weaker absorption system subsequently. Future simulations will clarify the contribution to 
quasar absorption lines of the entire population of photo-evaporating halos. 



6.5. Suppression of the Formation of Low Mass Galaxies 

At the end of overlap, the cosmic ionizing background increased sharply, and the IGM was 
heated by the ionizing radiation to a temperature ;> 10 4 K. Due to the substantial increase in 
the IGM temperature, the intergalactic Jeans mass increased dramatically, changing the minimum 
mass of forming galaxies (Rees 1986; Efstathiou 1992; Gnedin & Ostriker 1997; Miralda-Escude & 
Rees 1998). 

Gas infall depends sensitively on the Jeans mass. When a halo more massive than the Jeans 
mass begins to form, the gravity of its dark matter overcomes the gas pressure. Even in halos below 
the Jeans mass, although the gas is initially held up by pressure, once the dark matter collapses 
its increased gravity pulls in some gas (Haiman, Thoul, & Loeb 1996). Thus, the Jeans mass is 
generally higher than the actual limiting mass for accretion. Before reionization, the IGM is cold 
and neutral, and the Jeans mass plays a secondary role in limiting galaxy formation compared 
to cooling. After reionization, the Jeans mass is increased by several orders of magnitude due 
to the photoionization heating of the IGM, and hence begins to play a dominant role in limiting 
the formation of stars. Gas infall in a reionized and heated universe has been investigated in 
a number of numerical simulations. Thoul & Weinberg (1996) inferred, based on a spherically- 
symmetric collapse simulation, a reduction of ~ 50% in the collapsed gas mass due to heating, for 
a halo of circular velocity V c ~ 50 km s _1 at z = 2, and a complete suppression of infall below 
V c ~ 30 km s _1 . Kitayama & Ikeuchi (2000) also performed spherically-symmetric simulations 
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but included self-shielding of the gas, and found that it lowers the circular velocity thresholds by 
~ 5 km s _1 . Three dimensional numerical simulations (Quinn, Katz, & Efstathiou 1996; Weinberg, 
Hernquist, & Katz 1997; Navarro & Steinmetz 1997) found a significant suppression of gas infall in 
even larger halos (V c ~ 75 km s -1 ), but this was mostly due to a suppression of late infall at z <^ 2. 

When a volume of the IGM is ionized by stars, the gas is heated to a temperature Tigm ~ 10 4 
K. If quasars dominate the UV background at reionization, their harder photon spectrum leads to 
Tigm ~ 2 x 10 4 K. Including the effects of dark matter, a given temperature results in a linear 
Jeans mass (see § |3.1| ) corresponding to a halo circular velocity of 



Vj = 81' TlGM 
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kins" 1 , (79) 



where we used equation (25) and assumed fi = 0.6. In halos with V c > Vj, the gas fraction in 
infalling gas equals the universal mean of 0,b/il, m , but gas infall is suppressed in smaller halos. 
Even for a small dark matter halo, once it collapses to a virial overdensity of A c /f2^ relative to 
the mean, it can pull in additional gas. A simple estimate of the limiting circular velocity, below 
which halos have essentially no gas infall, is obtained by substituting the virial overdensity for the 
mean density in the definition of the Jeans mass. The resulting estimate is 

Hun = 34 ( - v /2 km s -i. (80) 



1.5 x 10 4 K 

This value is in rough agreement with the numerical simulations mentioned above (see also the 
related discussion in §3.2). 



Although the Jeans mass is closely related to the rate of gas infall at a given time, it does not 
directly yield the total gas residing in halos at a given time. The latter quantity depends on the 
entire history of gas accretion onto halos, as well as on the merger histories of halos, and an accurate 
description must involve a time-averaged Jeans mass. Gnedin (2000b) showed that the gas content 
of halos in simulations is well fit by an expression which depends on the filtering mass, a particular 
time-averaged Jeans mass (Gnedin & Hui 1998; see also § |3.ip . Gnedin (2000b) calculated the Jeans 
and filtering masses using the mean temperature in the simulation to define the sound speed, and 
found the following fit to the simulation results: 

M g = ^ , (si) 

9 [l + (2 1 /3_i)M c /M] 3 

where M g is the average gas mass of all objects with a total mass M, fb = £lb/Q>m is the universal 
baryon fraction, and the characteristic mass Mc is the total mass of objects which on average 
retain 50% of their gas mass. The characteristic mass was well fit by the filtering mass at a range 
of redshifts from z = 4uptoz~15. 



7. Feedback from Galactic Outflows 
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7.1. Propagation of Supernova Outflows in the IGM 



Star formation is accompanied by the violent death of massive stars in supernova explosions. 
In general, if each halo has a fixed baryon fraction and a fixed fraction of the baryons turns into 
massive stars, then the total energy in supernova outflows is proportional to the halo mass. The 
binding energies of both the supernova ejecta and of all the gas in the halo are proportional to the 
halo mass squared. Thus, outflows are expected to escape more easily out of low-mass galaxies, and 
to expel a greater fraction of the gas from dwarf galaxies. At high redshifts, most galaxies form in 
relatively low-mass halos, and the high halo merger rate leads to vigorous star formation. Thus, 
outflows may have had a great impact on the earliest generations of galaxies, with consequences that 
may include metal enrichment of the IGM and the disruption of dwarf galaxies. In this subsection 
we present a simple model for the propagation of individual supernova shock fronts in the IGM. 
We discuss some implications of this model, but we defer to the following subsection the brunt of 
the discussion of the cosmological consequences of outflows. 

For a galaxy forming in a given halo, the supernova rate is related to the star formation rate. 
In particular, for a Scalo (1998) initial stellar mass function, if we assume that a supernova is 
produced by each M > 8M & star, then on average one supernova explodes for every 126 M Q of star 
formation, expelling an ejecta mass of ~ 3M including ~ 1M of heavy elements. We assume 
that the individual supernovae produce expanding hot bubbles which merge into a single overall 
region delineated by an outwardly moving shock front. We assume that most of the baryons in the 
outflow lie in a thin shell, while most of the thermal energy is carried by the hot interior. The total 
ejected mass, which is lifted out of the halo by the outflow, equals a fraction _f gas of the total halo 
gas mass. The ejected mass includes some of the supernova ejecta itself. We let /eject denote the 
fraction of the supernova ejecta that winds up in the outflow (with /eject < 1 since some metals 
may be deposited in the disk and not ejected). Since at high redshift most of the halo gas is likely 
to have cooled onto a disk, we assume that the mass carried by the outflow remains constant until 
the shock front reaches the halo virial radius. We assume an average supernova energy of lO 51 ^! 
erg, a fraction / w i n d of which remains in the outflow after the outflow escapes from the disk. The 
outflow must overcome the gravitational potential of the halo, which we assume to have a Navarro, 
Frenk, & White (1997) density profile [NFW; see equation (§|) in §(U)]. Since the entire shell 
mass must be lifted out of the halo, we include the total shell mass as well as the total injected 
energy at the outset. This assumption is consistent with the fact that the burst of star formation in 
a halo is typically short compared to the total time for which the corresponding outflow expands. 

The escape of an outflow from an NFW halo depends on the concentration parameter cn of 
the halo. Simulations by Bullock et al. (2000) indicate that the concentration parameter decreases 
with redshift, and their results may be extrapolated to our regime of interest (i.e., to smaller halo 
masses and higher redshifts) by assuming that 
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Although we calculate below the dynamics of each outflow in detail, it is also useful to estimate 
which halos can generate large-scale outflows by comparing the kinetic energy of the outflow to the 
potential energy needed to completely escape (i.e., to infinite distance) from an NFW halo. We 
thus find that the outflow can escape from its originating halo if the circular velocity is below a 
critical value given by 

rr onn /^5l/wmd(??/0-l) _! 

V crit = 200 \ — r km s , (83) 

V /gas 9{CN) 

where the efficiency r] is the fraction of baryons incorporated in stars, and 

9{x) = (i + x)Hi + x )- x ■ (84) 

Note that the contribution to / gas of the supernova ejecta itself is 0.024ry/ e j ect , so the ejecta mass 
is usually negligible unless / gas 1%. Equation (|83|) can also be used to yield the maximum gas 
fraction / gas which can be ejected from halos, as a function of their circular velocity. Although this 
equation is most general, if we assume that the parameters f gas and / w j n d are independent of M 
and z then we can normalize them based on low-redshift observations. If we specify cn ~ 10 (with 
#(10) = 6.1) at z = 0, then setting £51 = 1 and 77 = 10% yields the required energy efficiency as a 
function of the ejected halo gas fraction: 



/wind — 1-5/, 
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A value of V a it ~ 100 km s _1 is suggested by several theoretical and observational arguments which 
are discussed in the next subsection. However, these arguments are not conclusive, and V cr n may 
differ from this value by a large factor, especially at high redshift (where outflows are observationally 
unconstrained at present). Note the degeneracy between / gas and / w i n( j which remains even if V CT it 
is specified. Thus, if V CI it ~ 100 km s~ x then a high efficiency / w i n d ~ 1 is required to eject most of 
the gas from all halos with V c < K:rit> but only / w j n d ~ 10% is required to eject 5-10% of the gas. 
The evolution of the outflow does depend on the value of / w i n d and not just the ratio / w ind//gas> 
since the shell accumulates material from the IGM which eventually dominates over the initial mass 
carried by the outflow. 

We solve numerically for the spherical expansion of a galactic outflow, elaborating on the basic 
approach of Tegmark, Silk, & Evrard (1993). We assume that most of the mass m carried along 
by the outflow lies in a thin, dense, relatively cool shell of proper radius R. The interior volume, 
while containing only a fraction /j n t < 1 of the mass m, carries most of the thermal energy in a 
hot, isothermal plasma of pressure pi n t and temperature T. We assume a uniform exterior gas, 
at the mean density of the universe (at each redshift), which may be neutral or ionized, and may 
exert a pressure p ex t as indicated below. We also assume that the dark matter distribution follows 
the NFW profile out to the virial radius, and is at the mean density of the universe outside the 
halo virial radius. Note that in reality an overdense distribution of gas as well as dark matter may 
surround each halo due to secondary infall. 
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The shell radius R in general evolves as follows: 

d 2 R „ 2p fdR r \ dm Gm /„.„, 1 \ 8 „ n 

where the right-hand-side includes forces due to pressure, sweeping up of additional mass, gravity, 
and a cosmological constant, respectively^. The shell is accelerated by internal pressure and de- 
celerated by external pressure, i.e., dp = pint — Pext- In the gravitational force, M(R) is the total 
enclosed mass, not including matter in the shell, and im is the effective contribution of the shell 
mass in the thin-shell approximation (Ostriker & McKee 1988). The interior pressure is determined 
by energy conservation, and evolves according to (Tegmark et al. 1993): 

dpint = L _ pint dR 
dt ~ 2ttR 3 R dt ' 1 ' 

where the luminosity L incorporates heating and cooling terms. We include in L the supernova 
luminosity L sn (during a brief initial period of energy injection), cooling terms L coo \, ionization 
Li on , and dissipation Ldi ss - For simplicity, we assume ionization equilibrium for the interior plasma, 
and a primordial abundance of hydrogen and helium. We include in L coo i all relevant atomic 
cooling processes in hydrogen and helium, i.e., collisional processes, Bremsstrahlung emission, and 
Compton cooling off the CMB. Compton scattering is the dominant cooling process for high-redshift 
outflows. We include in Li on only the power required to ionize the incoming hydrogen upstream, 
at the energy cost of 13.6 eV per hydrogen atom. The interaction between the expanding shell and 
the swept-up mass dissipates kinetic energy. The fraction of this energy which is re-injected into 
the interior depends on complex processes occurring near the shock front, including turbulence, 
non-equilibrium ionization and cooling, and so (following Tegmark et al. 1993) we let 

L ^ = \f^^§- HR ) > ( 88 ) 
where we set fd = 1 and compare below to the other extreme of fd = 0. 

In an expanding universe, it is preferable to describe the propagation of outflows in terms of 
comoving coordinates since, e.g., the critical result is the maximum comoving size of each outflow, 
since this size yields directly the total IGM mass which is displaced by the outflow and injected 
with metals. Specifically, we apply the following transformation (Shandarin 1980): 

di = a~ 2 dt, R = a' x R, p = a 5 p, p = a 3 p . (89) 

For = 0, Voit (1996) obtained (with the time origin t = at redshift z\): 

i- 2 



yj\ + Q m zi ~ v 7 ! + tt m z , (90) 



5 The last term, which is due to the cosmological constant, is an effective repulsion which arises in the Newtonian 
limit of the full equations of General Relativity. 
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while for m + = 1 there is no simple analytic expression. We set f5 = i?/f v i r , in terms of the 
virial radius r v ; r [equation fl24|)] of the source halo. We define ag as the ratio of the shell mass 
m to ^Trpbfy- 1T , where pb = Pb(z = 0) is the mean baryon density of the Universe at z = 0. More 
generally, we define 

as[P) ~ R3 \ 1 + ( a i - 1) /fp otherwise. lJij 

Here we assumed, as noted above, that the shell mass is constant until the halo virial radius is 
reached, at which point the outflow begins to sweep up material from the IGM. We thus derive the 
following equations: 

3 p %RHfo m 5((3) if/3<l 

(92) 




d 2 R 

di2 1 A p {dk ^ %RHfo m 5(p) + %RHfo b a s {l3) otherwise, 



along with 



3 < 93 > 

In the evolution equation for R, for j3 < 1 we assume for simplicity that the baryons are distributed 
in the same way as the dark matter, since in any case the dark matter halo dominates the overall 
gravitational potential. For /3 > 1, however, we correct (via the last term on the right-hand side) 
for the presence of mass in the shell, since at (3 3> 1 this term may become important. The 
(3 > 1 equation also includes the braking force due to the swept-up IGM mass. The enclosed mean 
overdensity for the NFW profile [equation (^§1) in §( |2.3D 1 surrounded by matter at the mean density 
is 

{ A e ln(l+c N /3)-c N /3/(l+c N /3) . f R -. 
ln(l+ CN )- CN /(l+ CN ) up^-t 
, \ 1 (94) 
[Tt~ l )W otherwise. 

The physics of supernova shells is discussed in Ostriker & McKee (1988) along with a number 
of analytical solutions. The propagation of cosmological blast waves has also been computed by 
Ostriker & Cowie (1981), Bertschinger (1985) and Carr & Ikeuchi (1985). Voit (1996) derived an 
exact analytic solution to the fluid equations which, although of limited validity, is nonetheless 
useful for understanding roughly how the outflow size depends on several of the parameters. The 
solution requires an idealized case of an outflow which at all times expands into a homogeneous 
IGM. Peculiar gravitational forces, and the energy lost in escaping from the host halo, are neglected, 
cooling and ionization losses are also assumed to be negligible, and the external pressure is not 
included. The dissipated energy is assumed to be retained, i.e., fd is set equal to unity. Under 
these conditions, the standard Sedov- Taylor self-similar solution (Sedov 1946, 1959; Taylor 1950) 
generalizes to the cosmological case as follows (Voit 1996): 

fl = ^y 5 fa/6 > (95) 
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where £ = 2.026 and Eq = Eq/(1 + z\) 2 in terms of the initial (i.e., at t = t = and z = z\) energy 
Eq. Numerically, the comoving radius is 



/ 0.022 E \ x / 5 / 10 j \ 2/5 



R = 280 ——5- " kpc . (96) 

\n b h 2 10 56 er g y yi + zi 10 10 yr/ 

In solving the equations described above, we assume that the shock front expands into a pre- 
ionized region which then recombines after a time determined by the recombination rate. Thus, the 
external pressure is included initially, it is turned off after the pre-ionized region recombines, and 
it is then switched back on at a lower redshift when the universe is reionized. When the ambient 
IGM is neutral and the pressure is off, the shock loses energy to ionization. In practice we find 
that the external pressure is unimportant during the initial expansion, although it is generally 
important after reionization. Also, at high redshift ionization losses are much smaller than losses 
due to Compton cooling. In the results shown below, we assume an instantaneous reionization at 
z = 9. 



Figure 27 shows the results for a starting redshift z = 15, for a halo of mass 5.4 x 10 7 Mq, 
stellar mass 8.OxlO 5 M , comoving t v [ t — 12 kpc, and circular velocity — 20 km/s. We show the 
shell comoving radius in units of the virial radius of the source halo (top panel) , and the physical 
peculiar velocity of the shock front (bottom panel). Results are shown (solid curve) for the standard 
set of parameters f- m t = 0.1, fd = 1, / w i n d = 75%, and / gas = 50%. For comparison, we show several 
cases which adopt the standard parameters except for no cooling (dotted curve), no reionization 
(short-dashed curve), fd = (long-dashed curve), or / w j n d = 15% and / gas = 10% (dot-short dashed 
curve). When reionization is included, the external pressure halts the expanding bubble. We freeze 
the radius at the point of maximum expansion (where dR/dt = 0), since in reality the shell will 
at that point begin to spread and fill out the interior volume due to small-scale velocities in the 
IGM. For the chosen parameters, the bubble easily escapes from the halo, but when f wm d and / gas 
are decreased the accumulated IGM mass slows down the outflow more effectively. In all cases the 
outflow reaches a size of 10-20 times r v ; r , i.e., 100-200 comoving kpc. If all the metals are ejected 
(i.e., /eject = l)j then this translates to an average metallicity in the shell of ~ l-5xl0 -3 in units 
of the solar metallicity (which is 2% by mass). The asymptotic size of the outflow varies roughly 
as /^ d , as predicted by the simple solution in equation (p5[), but the asymptotic size is rather 
insensitive to / gas (at a fixed / w i n d) since the outflow mass becomes dominated by the swept-up 
IGM mass once R ^ 4f v j r . With the standard parameter values (i.e., those corresponding to the 



solid curve), Figure |27| also shows (dot-long dashed curve) the Voit (1996) solution of equation 
(|95|). The Voit solution behaves similarly to the no-reionization curve at low redshift, although it 
overestimates the shock radius by ~ 30%, and the overestimate is greater compared to the more 
realistic case which does include reionization. 



Figure 28 shows different curves than Figure 27 but on an identical layout. A single curve 
starting at z = 15 (solid curve) is repeated from Figure [27], and it is compared here to outflows 
with the same parameters but starting at z = 20 (dotted curve), z = 10 (short-dashed curve), and 
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z = 5 (long-dashed curve). A V c = 20 km/s halo, with a stellar mass equal to 1.5% of the total 
halo mass, is chosen at the three higher redshifts, but at z = 5 a V c = 42 km/s halo is assumed. 
Because of the suppression of gas infall after reionization (§ |6.5| ), we assume that the z = 5 outflow 
is produced by supernovae from a stellar mass equal to only 0.3% of the total halo mass (with 
a similarly reduced initial shell mass), thus leading to a relatively small final shell radius. The 
main conclusion from both Figures is the following: In all cases, the outflow undergoes a rapid 
initial expansion over a fractional redshift interval 5z/z ~ 0.2, at which point the shell has slowed 
down to ~ 10 km/s from an initial 300 km/s. The rapid deceleration is due to the accumulating 
IGM mass. External pressure from the reionized IGM completely halts all high-redshift outflows, 
and even without this effect most outflows would only move at ~ 10 km/s after the brief initial 
expansion. Thus, it may be possible for high-redshift outflows to pollute the Lyman alpha forest 
with metals without affecting the forest hydrodynamically at z <^ 4. While the bulk velocities 
of these outflows may dissipate quickly, the outflows do sweep away the IGM and create empty 
bubbles. The resulting effects on observations of the Lyman alpha forest should be studied in detail 
(some observational signatures of feedback have been suggested recently by Theuns, Mo, &; Schaye 
2001). 

Barkana & Loeb (2001, in preparation) derive the overall filling factor of supernova bubbles 
based on this formalism. In the following subsection we survey previous analytic and numerical 
work on the collective astrophysical effects of galactic outflows. 

7.2. Effect of Outflows on Dwarf Galaxies and on the IGM 

Galactic outflows represent a complex feedback process which affects the evolution of cosmic gas 
through a variety of phenomena. Outflows inject hydrodynamic energy into the interstellar medium 
of their host galaxy. As shown in the previous subsection, even a small fraction of this energy suffices 
to eject most of the gas from a dwarf galaxy, perhaps quenching further star formation after the 
initial burst. At the same time, the enriched gas in outflows can mix with the interstellar medium 
and with the surrounding IGM, allowing later generations of stars to form more easily because 
of metal-enhanced cooling. On the other hand, the expanding shock waves may also strip gas in 
surrounding galaxies and suppress star formation. 

Dekel & Silk (1986) attempted to explain the different properties of diffuse dwarf galaxies in 
terms of the effect of galactic outflows (see also Larson 1974; Vader 1986, 1987). They noted the 
observed trends whereby lower-mass dwarf galaxies have a lower surface brightness and metallicity, 
but a higher mass-to-light ratio, than higher mass galaxies. They argued that these trends are 
most naturally explained by substantial gas removal from an underlying dark matter potential. 
Galaxies lying in small halos can eject their remaining gas after only a tiny fraction of the gas has 
turned into stars, while larger galaxies require more substantial star formation before the resulting 
outflows can expel the rest of the gas. Assuming a wind efficiency / W i n d ~ 100%, Dekel & Silk 
showed that outflows in halos below a circular velocity threshold of V cr n ~ 100 km/s have sufficient 
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Fig. 27. — Evolution of a supernova outflow from a z = 15 halo of circular velocity V c = 20 km/s. 
Plotted are the shell comoving radius in units of the virial radius of the source halo (top panel), 
and the physical peculiar velocity of the shock front (bottom panel). Results are shown for the 
standard parameters f- m t = 0.1, = 1, f w - m d = 75%, and / gas = 50% (solid curve). Also shown for 
comparison are the cases of no cooling (dotted curve), no reionization (short-dashed curve), fd = 
(long-dashed curve), or / w i n( j = 15% and / gas = 10% (dot-short dashed curve), as well as the simple 
Voit (1996) solution of equation ( |95| ) for the standard parameter set (dot-long dashed curve). In 
cases where the outflow halts, we freeze the radius at the point of maximum expansion. 
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Fig. 28. — Evolution of supernova outflows at different redshifts. The top and bottom panels are 
arranged similarly to Figure 27. The z = 15 outflow (solid curve) is repeated from Figure ^7[ and 
it is compared here to outflows with the same parameters but starting at z = 20 (dotted curve), 
z = 10 (short-dashed curve), and z = 5 (long-dashed curve). A V c = 20 km/s halo is assumed 
except for z = 5, in which case a V c = 42 km/s halo is assumed to produce the outflow (see text). 
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energy to expel most of the halo gas. Furthermore, cooling is very efficient for the characteristic 
gas temperatures associated with V cr it <^ 100 km/s halos, but it becomes less efficient in more 
massive halos. As a result, this critical velocity is expected to signify a dividing line between bright 
galaxies and diffuse dwarf galaxies. Although these simple considerations may explain a number of 
observed trends, many details are still not conclusively determined. For instance, even in galaxies 
with sufficient energy to expel the gas, it is possible that this energy gets deposited in only a small 
fraction of the gas, leaving the rest almost unaffected. 

Since supernova explosions in an inhomogeneous interstellar medium lead to complicated hy- 
drodynamics, in principle the best way to determine the basic parameters discussed in the previous 
subsection (/ W i n d, /gas, and /eject) is through detailed numerical simulations of individual galaxies. 
Mac Low & Ferrara (1999) simulated a gas disk within az = dark matter halo. The disk was 
assumed to be azimuthal and initially smooth. They represented supernovae by a central source of 
energy and mass, assuming a constant luminosity which is maintained for 50 million years. They 
found that the hot, metal-enriched ejecta can in general escape from the halo much more easily 
than the colder gas within the disk, since the hot gas is ejected in a tube perpendicular to the 
disk without displacing most of the gas in the disk. In particular, most of the metals were ex- 
pelled except for the case with the most massive halo considered (with 10 9 Mq in gas) and the 
lowest luminosity (10 37 erg/s, or a total injection of 2 x 10 52 erg). On the other hand, only a 
small fraction of the total gas mass was ejected except for the least massive halo (with 10 6 Mq in 
gas), where a luminosity of 10 38 erg/s or more expelled most of the gas. We note that beyond the 
standard issues of numerical resolution and convergence, there are several difficulties in applying 
these results to high-redshift dwarf galaxies. Clumping within the expanding shells or the ambient 
interstellar medium may strongly affect both the cooling and the hydrodynamics. Also, the effect 
of distributing the star formation throughout the disk is unclear since in that case several charac- 
teristics of the problem will change; many small explosions will distribute the same energy over a 
larger gas volume than a single large explosion [as in the Sedov- Taylor solution (Sedov 1946, 1959; 
Taylor 1950); see, e.g., equation (|95|) 1, and the geometry will be different as each bubble tries to 
dig its own escape route through the disk. Also, high-redshift disks should be denser by orders of 
magnitude than z = disks, due to the higher mean density of the universe at early times. Thus, 
further numerical simulations of this process are required in order to assess its significance during 
the reionization epoch. 

Some input on these issues also comes from observations. Martin (1999) showed that the 
hottest extended X-ray emission in galaxies is characterized by a temperature of ~ 10 6 ' 7 K. This 
hot gas, which is lifted out of the disk at a rate comparable to the rate at which gas goes into new 
stars, could escape from galaxies with rotation speeds of ^ 130 km/s. However, these results are 
based on a small sample that includes only the most vigorous star-forming local galaxies, and the 
mass-loss rate depends on assumptions about the poorly understood transfer of mass and energy 
among the various phases of the interstellar medium. 

Many authors have attempted to estimate the overall cosmological effects of outflows by com- 
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bining simple models of individual outflows with the formation rate of galaxies, obtained via semi- 
analytic methods (Couchman & Rees 1986; Blanchard, Valls-Gabaud, & Mamon 1992; Tegmark, 
et al. 1993; Voit 1996; Nath & Trentham 1997; Prunet & Blanchard 2000; Ferrara, Pettini, & 
Shchekinov 2000; Scannapieco & Broadhurst 2000) or numerical simulations (Gnedin & Ostriker 
1997; Gnedin 1998; Cen & Ostriker 1999; Aguirre et al. 2000a). The main goal of these calculations 
is to explain the characteristic metallicities of different environments as a function of redshift. For 
example, the IGM is observed to be enriched with metals at redshifts z <^ 5. Identification of C IV, 
Si IV and O VI absorption lines which correspond to Lya absorption lines in the spectra of high- 
redshift quasars has revealed that the low-density IGM has been enriched to a metal abundance 
(by mass) of Zi GM ~ lO- 2 - 5 ( ±a5 )Z , where Z & = 0.019 is the solar metallicity (Meyer & York 1987; 
Tytler et al. 1995; Songaila & Cowie 1996; Lu et al 1998; Cowie & Songaila 1998; Songaila 1997; 
Ellison et al. 2000; Prochaska & Wolfe 2000). The metal enrichment has been clearly identified 
down to H I column densities of ~ 10 14,5 cm -2 . The detailed comparison of cosmological hydrody- 
namic simulations with quasar absorption spectra has established that the forest of Lya absorption 
lines is caused by the smoothly-fluctuating density of the neutral component of the IGM (Cen et 
al. 1994; Zhang, Anninos & Norman 1995; Hernquist et al. 1996). The simulations show a strong 
correlation between the H I column density and the gas overdensity <5 gas (e.g., Dave et al. 1999), 
implying that metals were dispersed into regions with an overdensity as low as <5 gas ~ 3 or possibly 
even lower. 

In general, dwarf galaxies are expected to dominate metal enrichment at high-redshift for 
several reasons. As noted above and in the previous subsection, outflows can escape more easily 
out of the potential wells of dwarfs. Also, at high redshift, massive halos are rare and dwarf halos 
are much more common. Finally, as already noted, the Sedov- Taylor solution (Sedov 1946, 1959; 
Taylor 1950) [or equation (^5|)] implies that for a given total energy and expansion time, multiple 
small outflows (i.e., caused by explosions with a small individual energy release) fill large volumes 
more effectively than would a smaller number of large outflows. Note, however, that the strong 
effect of feedback in dwarf galaxies may also quench star formation rapidly and reduce the efficiency 
of star formation in dwarfs below that found in more massive galaxies. 

Cen & Ostriker (1999) showed via numerical simulation that metals produced by supernovae 
do not mix uniformly over cosmological volumes. Instead, at each epoch the highest density regions 
have much higher metallicity than the low-density IGM. They noted that early star formation occurs 
in the most overdense regions, which therefore reach a high metallicity (of order a tenth of the solar 
value) by z ~ 3, when the IGM metallicity is lower by 1-2 orders of magnitude. At later times, 
the formation of high-temperature clusters in the highest-density regions suppresses star formation 
there, while lower-density regions continue to increase their metallicity. Note, however, that the 
spatial resolution of the hydrodynamic code of Cen & Ostriker is a few hundred kpc, and anything 
occurring on smaller scales is inserted directly via simple parametrized models. Scannapieco & 
Broadhurst (2000) implemented expanding outflows within a numerical scheme which, while not a 
full gravitational simulation, did include spatial correlations among halos. They showed that winds 
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from low-mass galaxies may also strip gas from nearby galaxies (see also Scannapieco, Ferrara, 
& Broadhurst 2000), thus suppressing star formation in a local neighborhood and substantially 
reducing the overall abundance of galaxies in halos below a mass of ~ 1O 1O M0. Although quasars 
do not produce metals, they may also affect galaxy formation in their vicinity via energetic outflows 
(Efstathiou & Rees 1988; Babul & White 1991; Silk & Rees 1998; Natarajan, Sigurdsson, & Silk 
1998). 

Gnedin & Ostriker (1997) and Gnedin (1998) identified another mixing mechanism which, they 
argued, may be dominant at high redshift (z ;> 4). In a collision between two proto- galaxies, the 
gas components collide in a shock and the resulting pressure force can eject a few percent of the 
gas out of the merger remnant. This is the merger mechanism, which is based on gravity and 
hydrodynamics rather than direct stellar feedback. Even if supernovae inject most of their metals 
in a local region, larger-scale mixing can occur via mergers. Note, however, that Gnedin's (1998) 
simulation assumed a comoving star formation rate at z ;> 5 of ~ IMq per year per comoving Mpc 3 , 
which is 5-10 times larger than the observed rate at redshift 3-4 (§|8.1|). Aguirre et al. (2000a) used 
outflows implemented in simulations to conclude that winds of ~ 300 km/s at z <^ 6 can produce 
the mean metallicity observed at z ~ 3 in the Lyct forest. Aguirre et al. (2000b) explored another 
process, where metals in the form of dust grains are driven to large distances by radiation pressure, 
thus producing large-scale mixing without displacing or heating large volumes of IGM gas. The 
success of this mechanism depends on detailed microphysics such as dust grain destruction and 
the effect of magnetic fields. The scenario, though, may be directly testable because it leads to 
significant ejection only of elements which solidify as grains. 

Feedback from galactic outflows encompasses a large variety of processes and influences. The 
large range of scales involved, from stars or quasars embedded in the interstellar medium up to 
the enriched IGM on cosmological scales, make possible a multitude of different, complementary 
approaches, promising to keep galactic feedback an active field of research. 

8. Properties of the Expected Source Population 

8.1. The Cosmic Star Formation History 

One of the major goals of the study of galaxy formation is to achieve an observational de- 
termination and a theoretical understanding of the cosmic star formation history. By now, this 
history has been sketched out to a redshift z ~ 4 (see, e.g., the compilation of Blain et al. 1999a). 
This is based on a large number of observations in different wavebands. These include various 
ultraviolet /optical/near- infrared observations (Madau et al. 1996; Gallego et al. 1996; Lilly et al. 
1996; Connolly et al. 1997; Treyer et al. 1998; Tresse & Maddox 1998; Pettini et al. 1998a,b; Cowie, 
Songaila & Barger 1999; Gronwall 1999; Glazebrook et al. 1999; Yan et al. 1999; Flores et al. 1999; 
Steidel et al. 1999). At the shortest wavelengths, the extinction correction is likely to be large 
(a factor of ~ 5) and is still highly uncertain. At longer wavelengths, the star formation history 
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has been reconstructed from submillimeter observations (Blain et al. 1999b; Hughes et al. 1998) 
and radio observations (Cram 1998). In the submillimeter regime, a major uncertainty results 
from the fact that only a minor portion of the total far infrared emission of galaxies comes out in 
the observed bands, and so in order to estimate the star formation rate it is necessary to assume 
a spectrum based, e.g., on a model of the dust emission (see the discussion in Chapman et al. 
2000). In general, estimates of the star formation rate (hereafter SFR) apply locally-calibrated 
correlations between emission in particular lines or wavebands and the total SFR. It is often not 
possible to check these correlations directly on high-redshift populations, and together with the 
other uncertainties (extinction and incompleteness) this means that current knowledge of the star 
formation history must be considered to be a qualitative sketch only. Despite the relatively early 
state of observations, a wealth of new observatories in all wavelength regions promise to greatly 
advance the field. In particular, NGST will be able to detect galaxies and hence determine the 
star formation history out to z ^ 10. 

Hierarchical models have been used in many papers to match observations on star formation 
at z <J 4 (see, e.g. Baugh et al. 1998; Kauffmann & Chariot 1998; Somerville & Primack 1998, and 
references therein) . In this section we focus on theoretical predictions for the cosmic star formation 
rate at higher redshifts. The reheating of the IGM during reionization suppressed star formation 
inside the smallest halos (§ |6.5[) . Reionization is therefore predicted to cause a drop in the cosmic 
SFR. This drop is accompanied by a dramatic fall in the number counts of faint galaxies. Barkana 
& Loeb (2000b) argued that a detection of this fall in the faint luminosity function could be used 
to identify the reionization redshift observationally. 

A model for the SFR can be constructed based on the extended Press-Schechter theory. The 
starting point is the abundance of dark matter halos, obtained using the Press-Schechter model. 
The abundance of halos evolves with redshift as each halo gains mass through mergers with other 
halos. If dp[Mi,t\ — > M,t] is the probability that a halo of mass Mi at time t± will have merged 
to form a halo of mass between M and M + dM at time t > t\, then in the limit where t\ tends to 
t we obtain an instantaneous merger rate d 2 p[M 1 -> M,t]/(dM dt). This quantity was evaluated 
by Lacey & Cole [1993, their equation (2.18)], and it is the basis for modeling the rate of galaxy 
formation. 

Once a dark matter halo has collapsed and virialized, the two requirements for forming new 
stars are gas infall and cooling. We assume that by the time of reionization, photo-dissociation of 



molecular hydrogen (see §0) has left only atomic transitions as an avenue for efficient cooling. Be- 
fore reionization, therefore, galaxies can form in halos down to a circular velocity of V c ~ 17 km s _1 , 
where this limit is set by cooling. On the other hand, when a volume of the IGM is ionized by 
stars or quasars, the gas is heated and the increased pressure suppresses gas infall into halos with a 
circular velocity below V c ~ 80 km s" 1 , halting infall below V c ~ 30 km s _1 (§ |6.5| ). Since the sup- 
pression acts only in regions that have been heated, the reionization feedback on galaxy formation 
depends on the fraction of the IGM which is ionized at each redshift. In order to include a gradual 
reionization in the model, we take the simulations of Gnedin (2000a) as a guide for the redshift 
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interval of reionization. 

In general, new star formation in a given galaxy can occur either from primordial gas or from 
recycled gas which has already undergone a previous burst of star formation. The former occurs 
when a massive halo accretes gas from the IGM or from a halo which is too small to have formed 
stars. The latter occurs when two halos, in which a fraction of the gas has already turned to 
stars, merge and trigger star formation in the remaining gas. Numerical simulations of starbursts 
in interacting z = galaxies (e.g., Mihos & Hernquist 1994; 1996) found that a merger triggers 
significant star formation in a halo even if it merges with a much less massive partner. Preliminary 
results (Somerville 2000, private communication) from simulations of mergers at z ~ 3 find that 
they remain effective at triggering star formation even when the initial disks are dominated by gas. 
Regardless of the mechanism, we assume that feedback limits the star formation efficiency, so that 
only a fraction r/ of the gas is turned into stars. 

Given the SFR and the total number of stars in a halo of mass M, the luminosity and spectrum 
can be derived from an assumed stellar initial mass function. We assume an initial mass function 
which is similar to the one measured locally. If n(M) is the total number of stars with masses less 
than M, then the stellar initial mass function, normalized to a total mass of IMq, is (Scalo 1998) 



by the dense Lya forest. We do not, however, include dust extinction, which could be significant 
in some individual galaxies despite the low mean metallicity expected at high redshift. 

Much of the star formation at high redshift is expected to occur in low mass, faint galaxies, and 
even NGST may only detect a fraction of the total SFR. A realistic estimate of this fraction must 
include the finite resolution of the instrument as well as its detection limit for faint sources (Barkana 
& Loeb 2000a). We characterize the instrument's resolution by a minimum circular aperture of 
angular diameter 9 a . We describe the sensitivity of NGST by F£ s , the minimum spectral fluxQ, 
averaged over wavelengths 0.6-3.5/im, required to detect a point source (i.e., a source which is much 
smaller than 9 a ). For an extended source of diameter 9 S ^> 9 a , we assume that the signal-to-noise 
ratio can be improved by using a larger aperture, with diameter 9 S . The noise amplitude scales 
as the square root of the number of noise (sky) photons, or the square root of the corresponding 
sky area. Thus, the total flux needed for detection of an extended source at a given signal-to-noise 
threshold is larger than F£ s by a factor of 9 s /9 a . We adopt a simple interpolation formula between 




(97) 



where Mi = M/Mq. We assume a metallicity Z = 0.001, and use the stellar population model 
results of Leitherer et al. (1999) []. We also include a Lya cutoff in the spectrum due to absorption 



6 Model spectra of star-forming galaxies were obtained from http://www.stsci.edu/science/starburst99/ 
7 Note that F„ B is the total spectral flux of the source, not just the portion contained within the aperture. 
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the regimes of point-like and extended sources, and assume that a source is detectable if its flux is 
at least y/l + {d a /6 a )' 2 F£ s . 

We combine this result with a model for the distribution of disk sizes at each value of halo 
mass and redshift (§[5.11). We adopt a value of F^ s = 0.25 nJy [], assuming a deep 300-hour 
integration on an 8-meter NGST and a spectral resolution of 10:1. This resolution should suffice 
for a ~ 10% redshift measurement, based on the Lya cutoff. We also choose the aperture diameter 
to be 6 a = / /06, close to the expected NGST resolution at 2//m. 

Figure ^ shows our predictions for the star formation history of the universe, adopted from 
Figure 1 of Barkana & Loeb (2000b) with slight modifications (in the initial mass function and the 
values of the cosmological parameters). Letting z re i on denote the redshift at the end of overlap, 
we show the SFR for z TC i on = 6 (solid curves), z rc \ on = 8 (dashed curves), and 2: re i on = 10 (dotted 
curves). In each pair of curves, the upper one is the total SFR, and the lower one is the fraction 
detectable with NGST. The curves assume a star formation efficiency rj = 10% and an IGM 
temperature Tjgm = 2 x 10 4 K. Although photoionization directly suppresses new gas infall after 
reionization, it does not immediately affect mergers which continue to trigger star formation in gas 
which had cooled prior to reionization. Thus, the overall suppression is dominated by the effect 
on star formation in primordial (unprocessed) gas. The contribution from merger-induced star 
formation is comparable to that from primordial gas at z < z re i on , and it is smaller at z > z re i n- 
However, the recycled gas contribution to the detectable SFR is dominant at the highest redshifts, 
since the brightest, highest mass halos form in mergers of halos which themselves already contain 
stars. Thus, even though most stars at z > z re i on form out of primordial, zero-metallicity gas, a 
majority of stars in detectable galaxies may form out of the small gas fraction that has already 
been enriched by the first generation of stars. 



Points with error bars in Figure 29 are observational estimates of the cosmic SFR per comoving 
volume at various redshifts (as compiled by Blain et al. 1999a). We choose r\ = 10% to obtain a 
rough agreement between the models and these observations at z ~ 3-4. An efficiency of order 
this value is also suggested by observations of the metallicity of the Lya forest at z = 3 (Haiman 
& Loeb 1999b). The SFR curves are roughly proportional to the value of r]. Note that in reality r\ 
may depend on the halo mass, since the effect of supernova feedback may be more pronounced in 
small galaxies (§0). Figure 29 shows a sharp rise in the total SFR at redshifts higher than z rc [ on . 



Although only a fraction of the total SFR can be detected with NGST, the detectable SFR displays 
a definite signature of the reionization redshift. However, current observations at lower redshifts 
demonstrate the observational difficulty in measuring the SFR directly. The redshift evolution of 
the faint luminosity function provides a clearer, more direct observational signature. We discuss 
this topic next. 



3 We obtained the flux limit using the NGST calculator at jittp: / /www. ngst.stsci.edu /nms / main / 
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Fig. 29. — Redshift evolution of the SFR (in Mq per year per comoving Mpc 3 ), adopted from Fig- 
ure 1 of Barkana & Loeb (2000b) with slight modifications. Points with error bars are observational 
estimates (compiled by Blain et al. 1999a). Also shown are model predictions for a reionization 
redshift z Te i on = 6 (solid curves), 2: re i on = 8 (dashed curves), and z Te { on = 10 (dotted curves), with 
a star formation efficiency rj = 10%. In each pair of curves, the upper one is the total SFR, and 
the lower one is the fraction detectable with NGST at a limiting point source flux of 0.25 nJy We 
assume the ACDM model (with parameters given at the end of §[!]) . 
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8.2. Number Counts 

8.2.1. Galaxies 

As shown in the previous section, the cosmic star formation history should display a signature 
of the reionization redshift. Much of the increase in the star formation rate beyond the reionization 
redshift is due to star formation occurring in very small, and thus faint, galaxies. This evolution 
in the faint luminosity function constitutes the clearest observational signature of the suppression 
of star formation after reionization. 

Figure [30] shows the predicted redshift distribution in ACDM (with parameters given at the 
end of §|l|) of galaxies observed with NGST. The plotted quantity is dN/dz, where N is the number 
of galaxies per NGST field of view (4' x 4'). The model predictions are shown for a reionization 
redshift z r eion = 6 (solid curve), z r eion = 8 (dashed curve), and z rc i on = 10 (dotted curve), with a 
star formation efficiency rj = 10%. All curves assume a point-source detection limit of 0.25 nJy. 
This plot is updated from Fig ure 7 of Barkana & Loeb (2000a) in that redshifts above -z r eion 
included. 

Clearly, thousands of galaxies are expected to be found at high redshift. This will allow 
a determination of the luminosity function at many redshift intervals, and thus a measurement 
of its evolution. As the redshift is increased, the luminosity function is predicted to gradually 



change shape during the overlap era of reionization. Figure 31 shows the predicted evolution of the 
luminosity function for various values of z Te i on . This Figure is adopted from Figure 2 of Barkana 
& Loeb (2000b) with modifications (in the initial mass function, the values of the cosmological 
parameters, and the plot layout). All curves show d 2 N/(dz dlnF^ s ), where N is the total number 
of galaxies in a single field of view of NGST, and F^ s is the limiting point source flux at 0.6- 
3.5/um for NGST. Each panel shows the result for a reionization redshift z TC i on = 6 (solid curve), 
-Zreion = 8 (dashed curve), and z re i n = 10 (dotted curve). Figure [H| shows the luminosity function 
as observed at z = 5 (upper left panel) and (proceeding clockwise) at z = 7, z = 9, and z = 11. 
Although our model assigns a fixed luminosity to all halos of a given mass and redshift, in reality 
such halos would have some dispersion in their merger histories and thus in their luminosities. We 
thus include smoothing in the plotted luminosity functions. Note the enormous increase in the 
number density of faint galaxies in a pre-reionization universe. Observing this dramatic increase 
toward high redshift would constitute a clear detection of reionization and of its major effect on 
galaxy formation. 



8.2.2. Quasars 

Dynamical studies indicate that massive black holes exist in the centers of most nearby galaxies 
(Richstone et al. 1998; Kormendy & Ho 2000; Kormendy 2000, and references therein). This leads 
to the profound conclusion that black hole formation is a generic consequence of galaxy formation. 



- 91 - 




5 10 15 

Redshift 



Fig. 30. — Predicted redshift distribution of galaxies observed with NGST, adopted and modified 
from Figure 7 of Barkana & Loeb (2000a). The distribution in the ACDM model (with parameters 
given at the end of §|l|), with a star formation efficiency r/ = 10%, is shown for a reionization 
redshift z re ion = 6 (solid curve), z re \ on = 8 (dashed curve), and 2: re i on = 10 (dotted curve). The 
plotted quantity is dN/dz, where N is the number of galaxies per NGST field of view. All curves 
assume a limiting point source flux of 0.25 nJy. 
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Fig. 31. — Predicted luminosity function of galaxies at a fixed redshift, adopted and modified from 
Figure 2 of Barkana & Loeb (2000b). With rj = 10%, the curves show d 2 N/(dz dlnF? 8 ), where 
iV is the total number of galaxies in a single field of view of NGST, and F^ s is the limiting point 
source flux averaged over 0.6-3.5//m for NGST. The luminosity function is shown at z = 5, z = 7, 
z = 9, and z = 11, with redshift increasing clockwise starting with the upper left panel. Each case 
assumes the ACDM model (with parameters given at the end of §[l]) and a reionization redshift 
"Zreion = 6 (solid curves), z re ion = 8 (dashed curves), or z re ion = 10 (dotted curves). The expected 
NGST detection limit is shown by the vertical dashed line. 
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The suggestion that massive black holes reside in galaxies and power quasars dates back to the 
sixties (Zel'dovich 1964; Salpeter 1964; Lynden-Bell 1969). Efstathiou & Rees (1988) pioneered the 
modeling of quasars in the modern context of galaxy formation theories. The model was extended 
by Haehnelt & Rees (1993) who added more details concerning the black hole formation efficiency 
and lightcurve. Haiman h Loeb (1998) and Haehnelt, Natarajan, & Rees (1998) extrapolated 
the model to high redshifts. All of these discussions used the Press-Schechter theory to describe 
the abundance of galaxy halos as a function of mass and redshift. More recently, Kauffmann & 
Haehnelt (2000; also Haehnelt & Kauffmann 2000) embedded the description of quasars within 
semi-analytic modeling of galaxy formation, which uses the extended Press-Schechter formalism to 
describe the merger history of galaxy halos. 

In general, the predicted evolution of the luminosity function of quasars is constrained by the 
need to match the observed quasar luminosity function at redshifts z <J 5, as well as data from the 
Hubble Deep Field (HDF) on faint point-sources. Prior to reionization, we may assume that quasars 
form only in galaxy halos with a circular velocity ^ 10 km s" 1 (or equivalently a virial temperature 
^ 10 4 K), for which cooling by atomic transitions is effective. After reionization, quasars form only 
in galaxies with a circular velocity ;> 50 km s _1 , for which substantial gas accretion from the warm 
(~ 10 4 K) IGM is possible. The limits set by the null detection of quasars in the HDF are consistent 
with the number counts of quasars which are implied by these thresholds (Haiman, Madau, & Loeb 
1999). 

For spherical accretion of ionized gas, the bolometric luminosity emitted by a black hole has a 
maximum value beyond which radiation pressure prevents gas accretion. This Eddington luminosity 
(Eddington 1926) is derived by equating the radiative repulsive force on a free electron to the 
gravitational attractive force on an ion in the plasma, 

La T GM hh fi e m p 



4 7rr 2 c r 2 

where ot = 6.65 x 10 -25 cm 2 is the Thomson cross-section, [i e m p is the average ion mass per 
electron, and Mbh is the black hole mass. Since both forces scale as r~ 2 , the limiting Eddington 
luminosity is independent of radius r in the Newtonian regime, and for gas of primordial composition 
is given by, 

= ^cGM hh ^m p = iq46 / M^_\ 

a T \10 8 M Q J & 1 ; 

Generically, the Eddington limit applies to within a factor of order unity also to simple accretion 
flows in a non-spherical geometry (Frank, King, Sz Raine 1992). 

The total luminosity of a black hole is related to its mass accretion rate by the radiative 
efficiency, e, 

L = eM hh c 2 . (100) 

For accretion through a thin Keplerian disk onto a Schwarzschild (non-rotating) black hole, e = 
5.7%, while for a maximally rotating Kerr black hole, e = 42% (Shapiro & Teukolsky 1983, p. 
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429). The thin disk configuration, for which these high radiative efficiencies are attainable, exists 
for L disk £ 0.5L E (Laor & Netzer 1989) . 

The accretion time can be defined as 



M bh \o.ij \l e j 

This time is comparable to the dynamical time inside the central kpc of a typical galaxy, t^ yn ~ 
(1 kpc/100 km s" 1 ) = 10 7 yr. As long as its fuel supply is not limited and e is constant, a black 
hole radiating at the Eddington limit will grow its mass exponentially with an e-folding time equal 
to r. The fact that r is much shorter than the age of the universe even at high redshift implies 
that black hole growth is mainly limited by the feeding rate Mbh(0> or by the total fuel reservoir, 
and not by the Eddington limit. 

The "simplest model" for quasars involves the following three assumptions (Haiman & Loeb 
1998): 

(i) A fixed fraction of the baryons in each "newly formed" galaxy ends up making a central black 
hole. 

(ii) Each black hole shines at its maximum (Eddington) luminosity for a universal amount of time. 

(iii) All black holes share the same emission spectrum during their luminous phase. 

Note that these assumptions relate only to the most luminous phase of the black hole accretion 
process, and they may not be valid during periods when the radiative efficiency or the mass accretion 
rate is very low. Such periods are not of interest here since they do not affect the luminosity function 
of bright quasars, which is the observable we wish to predict. The first of the above assumptions is 
reasonable as long as the fraction of virialized baryons in the universe is much smaller than unity; it 
does not include a separate mechanism for fueling black hole growth during mergers of previously- 
formed galaxies, and thus, under this assumption, black holes would not grow in mass once most of 
the baryons were virialized. The second hypothesis is motivated by the fact that for a sufficiently 
high fueling rate (which may occur in the early stage of the collapse/merger of a galaxy), quasars 
are likely to shine at their maximum possible luminosity. The resulting luminosity should be close 
to the Eddington limit over a period of order r. The third assumption can be implemented by 
incorporating the average quasar spectrum measured by Elvis et al. (1994). 

At high redshifts the number of "newly formed" galaxies can be estimated based on the time- 
derivative of the Press-Schechter mass function, since the collapsed fraction of baryons is small and 
most galaxies form out of the unvirialized IGM. Haiman & Loeb (1998, 1999a) have shown that 
the above simple prescription provides an excellent fit to the observed evolution of the luminosity 
function of bright quasars between redshifts 2.6 < z < 4.5 (see the analytic description of the 
existing data in Pei 1995). The observed decline in the abundance of bright quasars (Schneider, 
Schmidt, & Gunn 1991; Pei 1995) results from the deficiency of massive galaxies at high redshifts. 
Consequently, the average luminosity of quasars declines with increasing redshift. The required 
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ratio between the mass of the black hole and the total baryonic mass inside a halo is / -^gas — 
10~ 3 ' 2 r2 m /r2fc = 5.5 x 10~ 3 , comparable to the typical value of ~ 2-6 x 10~ 3 found for the ratio of 
black hole mass to spheroid mass in nearby elliptical galaxies (Magorrian et al. 1998; Kormendy 



2000). The required lifetime of the bright phase of quasars is ~ 10 yr. Figure 32 shows the most 
recent prediction of this model (Haiman & Loeb 1999a) for the number counts of high-redshift 
quasars, taking into account the above-mentioned thresholds for the circular velocities of galaxies 
before and after reionization0. 

We do, however, expect a substantial intrinsic scatter in the ratio Mbh/M gas . Observationally, 
the scatter around the average value of log 10 (Mbh/-^) is 0.3 (Magorrian et al. 1998), while the 
standard deviation in log 10 (Mbh/Mg as ) has been found to be a ~ 0.5. Such an intrinsic scatter 
would flatten the predicted quasar luminosity function at the bright end, where the luminosity 
function is steeply declining. However, Haiman & Loeb (1999a) have shown that the flattening 
introduced by the scatter can be compensated for through a modest reduction in the fitted value 
for the average ratio between the black hole mass and halo mass by ~ 50% in the relevant mass 
range (10 s M £ M bh £ 10 10 M ). 

In reality, the relation between the black hole and halo masses may be more complicated than 
linear. Models with additional free parameters, such as a non-linear (mass and redshift dependent) 
relation between the black hole and halo mass, can also produce acceptable fits to the observed 
quasar luminosity function (Haehnelt et al. 1998). The nonlinearity in a relation of the type 
Mbh -^haio w ith a > 1, may be related to the physics of the formation process of low-luminosity 
quasars (Haehnelt et al. 1998; Silk & Rees 1998), and can be tuned so as to reproduce the black hole 
reservoir with its scatter in the local universe (Cattaneo, Haehnelt, & Rees 1999). Recently, a tight 
correlation between the masses of black holes and the velocity dispersions of the bulges in which 
they reside, a, was identified in nearby galaxies. Ferrarese & Merritt (2000; see also Merritt &: 
Ferrarese 2001) inferred a correlation of the type Mbh oc (j 4 - 72 ^ - 36 ^ based on a selected sample of a 
dozen galaxies with reliable M b h estimates, while Gebhardt et al. (2000a, b) have found a somewhat 
shallower slope, Mbh oc cr 3 - 75 ^ - 3 ) based on a significantly larger sample. A non-linear relation of 
Mbh oc cr 5 oc Mj^j has been predicted by Silk & Rees (1998) based on feedback considerations, 
but the observed relation also follows naturally in the standard semi-analytic models of galaxy 
formation (Haehnelt & Kauffmann 2000). 

Figure |32] shows the predicted number counts in the "simplest model" described above (Haiman 
& Loeb 1999a), normalized to a 5' x 5' field of view. Figure |3^ shows separately the number per 
logarithmic flux interval of all objects with redshifts z > 5 (thin lines), and z > 10 (thick lines). 
The number of detectable sources is high; NGST will be able to probe of order 100 quasars at 
z > 10, and ~ 200 quasars at z > 5 per 5' x 5' field of view. The bright-end tail of the number 
counts approximately follows a power law, with dN/dF u oc F~ 2 ^ . The dashed lines show the 
corresponding number counts of "star-clusters" , assuming that each halo shines due to a starburst 



9 Note that the post-reionization threshold was not included in the original discussion of Haiman & Loeb (1998). 
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that converts a fraction of 2% (long-dashed) or 20% (short-dashed) of the gas into stars. 

Similar predictions can be made in the X-ray regime. Figure |33] shows the number counts of 
high-redshift X-ray quasars in the above "simplest model". This model fits the X-ray luminosity 
function of quasars at z ~ 3.5 as observed by ROSAT (Miyaji, Hasinger, & Schmidt 2000), using 
the same parameters necessary to fit the optical data (Pei 1995). Deep optical or infrared follow- 
ups on deep images taken with the Chandra X-ray satellite (CXO; see, e.g., Mushotzky et al. 2000; 
Barger et al. 2001; Giacconi et al. 2000) may be used to test these predictions in the relatively near 
future. 

The "simplest model" mentioned above predicts that black holes and stars make comparable 
contributions to the ionizing background prior to reionization. Consequently, the reionization of 
hydrogen and helium is predicted to occur roughly at the same epoch. A definitive identification of 
the He II reionization redshift will provide another powerful test of this model. Further constraints 
on the lifetime of the active phase of quasars may be provided by future measurements of the 
clustering properties of quasars (Haehnelt et al. 1998; Martini &; Weinberg 2001; Haiman & Hui 
2000). 



8. 2. 3. Supernovae 

The detection of galaxies and quasars becomes increasingly difficult at a higher redshift. This 
results both from the increase in the luminosity distance and the decrease in the average galaxy 
mass with increasing redshift. It therefore becomes advantageous to search for transient sources, 
such as supernovae or 7-ray bursts, as signposts of high-redshift galaxies (Miralda-Escude & Rees 
1997). Prior to or during the epoch of reionization, such sources are likely to outshine their host 
galaxies. 



The metals detected in the IGM (see §7^) signal the existence of supernova (SN) explosions at 
redshifts z ^ 5. Since each SN produces an average of ~ IMq of heavy elements (Woosley & Weaver 
1995), the inferred metallicity of the IGM, Ziqm, implies that there should be a supernova at z ;> 5 
for each ~ 1.7 x 10 4 M Q x (Zigm/IO -2,5 ^) -1 of baryons in the universe. We can therefore estimate 
the total supernova rate, on the entire sky, necessary to produce these metals at z ~ 5. Consider all 
SNe which are observed over a time interval At on the whole sky. Due to the cosmic time dilation, 
they correspond to a narrow redshift shell centered at the observer of proper width cAt/(l + z). In 
a flat f2 m = 0.3 cosmology, the total mass of baryons in a narrow redshift shell of width cAt/ (1 + z) 
around z = 5 is ~ [4vr(l + z)-' 3 (1.8c/ H ) 2 cAt] x [n b (3H$ /8ttG)(1 + z) 3 ] = 4.9 c 3 Vt b At/G. Hence, for 
h = 0.7 the total supernova rate across the entire sky at z ^ 5 is estimated to be (Miralda-Escude 
& Rees 1997), 



:v - * u,N (.i^y v, " l - (H,2i 



or roughly one SN per square arcminute per year. 
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Fig. 32. — Infrared number counts of quasars (averaged over the wavelength interval of 1-3.5/im) 
based on the "simplest quasar model" of Haiman & Loeb (1999b). The solid curves refer to quasars, 
while the long/short dashed curves correspond to star clusters with low/high normalization for the 
star formation efficiency. The curves labeled "5" or "10" show the cumulative number of objects 
with redshifts above z = 5 or 10. 
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Fig. 33. — Total number of quasars with redshift exceeding z = 5, z = 7, and z = 10 as a function 
of observed X-ray flux in the CXO detection band (from Haiman & Loeb 1999a). The numbers are 
normalized per 17' x 17' area of the sky. The solid curves correspond to a cutoff in circular velocity 
for the host halos of v c i TC > 50 km s \ the dashed curves to a cutoff of v c i TC > 100 kms" 1 . The 
vertical dashed line shows the CXO sensitivity for a 5a detection of a point source in an integration 
time of 5 x 10 5 seconds. 
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The actual SN rate at a given observed flux threshold is determined by the star formation rate 
per unit comoving volume as a function of redshift and the initial mass function of stars (Madau, 
della Valle, & Panagia 1998; Woods & Loeb 1998; Sullivan et al. 2000). To derive the relevant 
expression for flux-limited observations, we consider a general population of transient sources which 
are standard candles in peak flux and are characterized by a comoving rate per unit volume R{z). 
The observed number of new events per unit time brighter than flux F v at observed wavelength A 
for such a population is given by 

N(F U ;X)= R{z)(l + z)- 1 (dV c /dz)dz, (103) 

J o 

where z ma , x (F u ,\) is the maximum redshift at which a source will appear brighter than F v at 
an observed wavelength A = c/u, and dV c is the cosmology-dependent comoving volume element 
corresponding to a redshift interval dz. The above integrand includes the (1 + z) reduction in the 
apparent rate due to the cosmic time dilation. 

Figure shows the predicted SN rate as a function of limiting flux in various bands (Woods & 
Loeb 1998), based on the comoving star formation rate as a function of redshift that was determined 
empirically by Madau (1997). The actual star formation rate may be somewhat higher due to 
corrections for dust extinction (for a recent compilation of current data, see Blain & Natarajan 
2000). The dashed lines correspond to Type la SNe and the dotted lines to Type II SNe. For 
comparison, the solid lines indicate two crude estimates for the rate of 7-ray burst afterglows, 
which are discussed in detail in the next section. 



Equation ( |103j ) is appropriate for a threshold experiment, one which monitors the sky contin- 
uously and triggers when the detected flux exceeds a certain value, and hence identifies the most 
distant sources only when they are near their peak flux. For search strategies which involve taking a 
series of "snapshots" of a field and looking for variations in the flux of sources in successive images, 
one does not necessarily detect most sources near their peak flux. In this case, the total number of 
events (i.e., not per unit time) brighter than F v at observed wavelength A is given by 



00 



N(F U ; A) = / R(z%(z;F u ,X)(dV c /dz)dz, (104) 







where t*(z;F u , A) is the rest-frame duration over which an event will be brighter than the limiting 
flux F v at redshift z. This is a naive estimate of the so-called "control time" ; in practice, the effective 
duration over which an event can be observed is shorter, owing to the image subtraction technique, 
host galaxy magnitudes, and a number of other effects which reduce the detection efficiency (Pain 



et al. 1996). Figure 35 shows the predicted number counts of SNe as a function of limiting flux for 



the parameters used in Figure |34| (Woods & Loeb 1998). 

Supernovae also produce dust which could process the emission spectrum of galaxies. Al- 
though produced in galaxies, the dust may be expelled together with the metals out of galaxies 
by supernova-driven winds. Loeb & Haiman (1997) have shown that if each supernova produces 
~ 0.3 Mq of galactic dust, and some of the dust is expelled together with metals out of the shallow 
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Fig. 34. — Predicted cumulative rate N(> F v ) per year per square degree of supernovae at four 
wavelengths, corresponding to the K, R, B, and U bands (from Woods & Loeb 1998). The broken 
lines refer to different supernova types, namely SNe la (dashed curves) and SNe II (dotted curves). 
For comparison, the solid curves show estimates for the rates of gamma-ray burst (GRB) afterglows; 
the lower solid curve assumes the best-fit rate and luminosity for GRB sources which trace the star 
formation history (Wijers et al. 1998), while the upper solid curve assumes the best-fit values for a 
non-evolving GRB population. 
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Fig. 35. — Cumulative number counts N(> F v ) per square degree (from Woods & Loeb 1998). The 
notation is the same as in Figure 153. 
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potential wells of the early dwarf galaxies, then the optical depth for extinction by intergalactic 
dust may reach a few tenths at z ~ 10 for observed wavelengths of ~ 0.5-1 /xm [see Todini &: 
Ferrara (2000) for a detailed discussion on the production of dust in primordial Type II SNe]. The 
opacity in fact peaks in this wavelength band since at z ~ 10 it corresponds to rest-frame UV, 
where normal dust extinction is most effective. In these estimates, the amplitude of the opacity is 
calibrated based on the observed metallicity of the IGM at z <; 5. The intergalactic dust absorbs 
the UV background at the reionization epoch and re-radiates it at longer wavelengths. The flux 
and spectrum of the infrared background which is produced at each redshift depends sensitively on 
the distribution of dust around the ionizing sources, since the deviation of the dust temperature 
from the microwave background temperature depends on the local flux of UV radiation that it is 
exposed to. For reasonable choices of parameters, dust could lead to a significant spectral distortion 
of the microwave background spectrum that could be measured by a future spectral mission, going 
beyond the upper limit derived by the COBE satellite (Fixsen et al. 1996). 

The metals produced by supernovae may also yield strong molecular line emission. Silk & 
Spaans (1997) pointed out that the rotational line emission of CO by starburst galaxies is enhanced 
at high redshift due to the increasing temperature of the cosmic microwave background, which 
affects the thermal balance and the level populations of the atomic and molecular species. They 
found that the future Millimeter Array (MMA) could detect a starburst galaxy with a star formation 
rate of ~ 30 M yr _1 equally well at z = 5 and z = 30 because of the increasing cosmic microwave 
background temperature with redshift. Line emission may therefore be a more powerful probe of 
the first bright galaxies than continuum emission by dust. 

8.2.4- Gamma Ray Bursts 

The past decade has seen major observational breakthroughs in the study of Gamma Ray Burst 
(GRB) sources. The Burst and Transient Source Experiment (BATSE) on board the Compton 
Gamma Ray Observatory (Meegan et al. 1992) showed that the GRB population is distributed 
isotropically across the sky, and that there is a deficiency of faint GRBs relative to a Euclidean 
distribution. These were the first observational clues indicating a cosmological distance scale for 
GRBs. The localization of GRBs by X-ray observations with the BeppoSAX satellite (Costa et al. 
1997) allowed the detection of afterglow emission at optical (e.g., van Paradijs et al. 1997, 2000) and 
radio (e.g., Frail et al. 1997) wavelengths up to more than a year following the events (Fruchter et 
al. 1999; Frail et al. 2000). The afterglow emission is characterized by a broken power-law spectrum 
with a peak frequency that declines with time. The radiation is well-fitted by a model consisting 
of synchrotron emission from a decelerating blast wave (Blandford k, McKee 1976), created by the 
GRB explosion in an ambient medium, with a density comparable to that of the interstellar medium 
of galaxies (Waxman 1997; Sari, Piran, & Narayan 1998; Wijers & Galama 1999; Meszaros 1999; 
but see also Chevalier & Li 2000). The detection of spectral features, such as metal absorption lines 
in some optical afterglows (Metzger et al. 1997) and emission lines from host galaxies (Kulkarni et 



- 103 - 



al. 2000), allowed an unambiguous identification of cosmological distances to these sources. 

The nature of the central engine of GRBs is still unknown. Since the inferred energy release, 
in cases where it can be securely calibrated (Freedman & Waxman 2001; Frail et al. 2000), is 
comparable to that in a supernova, ~ 10 51 erg, most popular models relate GRBs to stellar remnants 
such as neutron stars or black holes (Eichler et al. 1989; Narayan, Paczyhski, & Piran 1992; 
Paczyhski 1991; Usov 1992; Mochkovitch et al. 1993; Paczyhski 1998; MacFadyen & Woosley 1999). 
Recently it has been claimed that the late evolution of some rapidly declining optical afterglows 
shows a component which is possibly associated with supernova emission (e.g., Bloom et al. 1999; 
Reichart 1999). If the supernova association is confirmed by detailed spectra of future afterglows, 
the GRB phenomenon will be linked to the terminal evolution of massive stars. 

Any association of GRBs with the formation of single compact stars implies that the GRB 
rate should trace the star formation history of the universe (Totani 1997; Sahu et al. 1997; Wijers 
et al. 1998; but see Krumholz, Thorsett & Harrison 1998). Owing to their high brightness, GRB 
afterglows could in principle be detected out to exceedingly high redshifts. Just as for quasars, 
the broad-band emission of GRB afterglows can be used to probe the absorption properties of 
the IGM out to the reionization epoch at redshift z ~ 10. Lamb & Reichart (2000) extrapolated 
the observed gamma-ray and afterglow spectra of known GRBs to high redshifts and emphasized 
the important role that their detection could play in probing the IGM (see also Miralda-Escude 
1998). In particular, the broad-band afterglow emission can be used to probe the ionization and 
metal-enrichment histories of the intervening IGM during the epoch of reionization. 

Ciardi & Loeb (2000) showed that unlike other sources (such as galaxies and quasars), which 
fade rapidly with increasing redshift, the observed infrared flux from a GRB afterglow at a fixed 
observed age is only a weak function of its redshift (Figure |36| ). A simple scaling of the long- 
wavelength spectra and the temporal evolution of afterglows with redshift implies that at a fixed 
time-lag after the GRB in the observer's frame, there is only a mild change in the observed flux 
at infrared or radio wavelengths with increasing redshift. This results in part from the fact that 
afterglows are brighter at earlier times, and that a given observed time refers to an earlier intrinsic 
time in the source frame as the source redshift increases. The "apparent brightening" of GRB 
afterglows with redshift could be further enhanced by the expected increase with redshift of the 
mean density of the interstellar medium of galaxies (Wood & Loeb 2000). Figure |37] shows the 
expected number counts of GRB afterglows, assuming that the GRB rate is proportional to the star 
formation rate and that the characteristic energy output of GRBs is ~ 10 52 erg and is isotropic. The 
figure implies that at any time there should be of order ~ 15 GRBs with redshifts z ^ 5 across the 
sky which are brighter than ~ 100 nJy at an observed wavelength of ~ 2/um. The infrared spectrum 
of these sources could be measured with NGST as a follow-up on their early X-ray localization with 
7-ray or X-ray detectors. Prior to reionization, the spectrum of GRB afterglows could reveal the 
long sought-after Gunn-Peterson trough (Gunn & Peterson 1965) due to absorption by the neutral 
IGM. 
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The predicted GRB rate and flux are subject to uncertainties regarding the beaming of the 
emission. The beaming angle may vary with observed time due to the decline with time of the 
Lorentz factor j(t) of the emitting material. As long as the Lorentz factor is significantly larger 
than the inverse of the beaming angle (i.e., 7 ^ 8 the afterglow flux behaves as if it were 
emitted by a spherically-symmetric fireball with the same explosion energy per unit solid angle. 
However, the lightcurve changes as soon as 7 declines below 9~ l , due to the lateral expansion of 
the jet (Rhoads 1997, 1999a, b; Panaitescu & Meszaros 1999). Finally, the isotropization of the 
energy ends when the expansion becomes sub-relativistic, at which point the remnant recovers 
the spherically-symmetric Sedov- Taylor solution (Sedov 1946, 1959; Taylor 1950) with the total 
remaining energy. When 7 ~ 1, the emission occurs from a roughly spherical fireball with the 
effective explosion energy per solid angle reduced by a factor of (27r# 2 /47r) relative to that at early 
times, representing the fraction of sky around the GRB source which is illuminated by the initial 
two (opposing) jets of angular radius 6 (see Ciardi &; Loeb 2000 for the impact of this effect on 
the number counts). The calibration of the GRB event rate per comoving volume, based on the 
number counts of GRBs (Wijers et al. 1998), is inversely proportional to this factor. 

The main difficulty in using GRBs as probes of the high-redshift universe is that they are rare, 
and hence their detection requires surveys which cover a wide area of the sky. The simplest strategy 
for identifying high-redshift afterglows is through all-sky surveys in the 7-ray or X-ray regimes. In 
particular, detection of high-redshift sources will become feasible with the high trigger rate provided 



by the forthcoming Swift satellite, to be launched in 2003 (see http:/ /swift.gsfc.nasa.gov/ , for more 
details). Swift is expected to localize ~300 GRBs per year, and to re-point within 20-70 seconds 
its on-board X-ray and UV-optical instrumentation for continued afterglow studies. The high- 
resolution GRB coordinates obtained by Swift will be transmitted to Earth within ~50 seconds. 
Deep follow-up observations will then be feasible from the ground or using the highly-sensitive 
infrared instruments on board NGST. Swift will be sufficiently sensitive to trigger on the 7-ray 
emission from GRBs at redshifts z ;> 10 (Lamb & Reichart 2000). 



8.3. Distribution of Disk Sizes 



Given the distribution of disk sizes at each value of halo mass and redshift (§ |5.lQ and the 
number counts of galaxies (§ |8,2.1 ), we derive the predicted size distribution of galactic disks. Note 



that although frequent mergers at high redshift may disrupt these disks and alter the morphologies 
of galaxies, the characteristic sizes of galaxies will likely not change dramatically. We show in 
Figure |3^ [an updated version of Figure 6 of Barkana Sz Loeb (2000a)] the distribution of galactic 
disk sizes at various redshifts, in the ACDM model (with parameters given at the end of §[l]). Given 
6 in arcseconds, each curve shows the fraction of the total number counts contributed by sources 
larger than 9. The diameter 9 is measured out to one exponential scale length. We show three pairs 
of curves, at z = 2, z = 5 and z = 10 (from right to left). Each pair includes the distribution for 
all galaxies (dashed line), and for galaxies detectable by NGST (solid line) with a limiting point 
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source flux of 0.25 nJy and with an efficiency rj = 10% assumed for the galaxies. The vertical 
dotted line indicates the expected NGST resolution of 0'.'06. 
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Fig. 36. — Observed flux from a 7-ray burst afterglow as a function of redshift (from Ciardi & 
Loeb 2000). The two sets of curves refer to a photon frequency v = 6 x 10 14 Hz (X b s = 5000 A, 
thin lines) and v = 1.5 x 10 14 Hz (X b s = 2/j,m, thick lines). Each set shows different observed 
times after the GRB trigger; from top to bottom: 1 hour (solid line), 1 day (dotted) and 10 days 
(dashed). The sharp suppression for 5000 A at z ^ 4.5 is due to IGM absorption. 

Among detectable galaxies, the typical diameter decreases from / /22 at z = 2 to C/10 at z = 5 
and / .'05 at z = 10. Note that in ACDM (with Q m = 0.3) the angular diameter distance (in units 
of c/Hq) actually decreases from 0.40 at z = 2 to 0.30 at z = 5 and 0.20 at z = 10. Galaxies are still 
typically smaller at the higher redshifts because a halo of a given mass is denser, and thus smaller, 
at higher redshift, and furthermore the typical halo mass is larger at low redshift due to the growth 
of cosmic structure with time. At z = 10, the distribution of detectable galaxies is biased, relative 
to the distribution of all galaxies, toward large galaxies, since NGST can only detect the brightest 
galaxies. The brightest galaxies tend to lie in the most massive and therefore largest halos, and 
this trend dominates over the higher detection threshold needed for an extended source compared 
to a point source (§ |8.1|) . 

Clearly, the angular resolution of NGST will be sufficiently high to resolve most galaxies. For 
example, NGST should resolve approximately 35% of z = 10 galaxies, 80% of z = 5 galaxies, 
and all but 1% of z = 2 galaxies. This implies that the shapes of these high-redshift galaxies 
can be studied with NGST. It also means that the high resolution of NGST is crucial in making 
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Fig. 37. — Predicted number of GRB afterglows per square degree with observed flux greater than 
F, at several different observed wavelengths (from Ciardi & Loeb 2000). From right to left, the 
observed wavelength equals 10 cm, 1 mm, 2 /jm and 5000 A. 



the majority of sources on the sky useful for weak lensing studies (although a mosaic of images is 
required for good statistics; see also the following subsection). 



8.4. Gravitational Lensing 

Detailed studies of gravitational lenses have provided a wealth of information on galaxies, both 
through modeling of individual lens systems (e.g., Schneider, Ehlers, & Falco 1992; Blandford & 
Narayan 1992) and from the statistical properties of multiply imaged sources (e.g., Turner, Ostriker, 
& Gott 1984; Maoz & Rix 1993; Kochanek 1996). 

The ability to observe large numbers of high-redshift objects promises to greatly extend gravi- 
tational lensing studies. Due to the increased path length along the line of sight to the most distant 
sources, their probability for being lensed is expected to be the highest among all possible sources. 
Sources at z > 10 will often be lensed by z > 2 galaxies, whose masses can then be determined with 
lens modeling. Similarly, the shape distortions (or weak lensing) caused by foreground clusters of 
galaxies will be used to determine the mass distributions of less massive and higher redshift clusters 
than currently feasible. In addition, it will be fruitful to exploit the magnification of the sources to 
resolve and study more distant galaxies than otherwise possible. 

These applications have been explored by Schneider & Kneib (1998), who investigated weak 
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Fig. 38. — Distribution of galactic disk sizes at various redshifts, in the ACDM model (with 
parameters given at the end of §|]), adopted and modified from Figure 6 of Barkana & Loeb (2000a). 
Given 9 in arcseconds, each curve shows the fraction of the total number counts contributed by 
sources larger than 0. The diameter 6 is measured out to one exponential scale length. We show 
three pairs of curves, at z = 2, z = 5 and z = 10 (from right to left). Each pair includes the 
distribution for all galaxies (dashed line), and for galaxies detectable by NGST (solid line) with 
a limiting point source flux of 0.25 nJy and with an efficiency rj = 10% assumed for the galaxies. 
The vertical dotted line indicates the expected NGST resolution of / /06. 
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lensing, and by Barkana & Loeb (2000a), who focused on strong lensing. Schneider & Kneib (1998) 
noted that the ability of NGST to take deeper exposures than is possible with current instruments 
will increase the observed density of sources on the sky, particularly of those at high redshifts. The 
large increase (by ~ 2 orders of magnitude over current surveys) may allow such applications as 
a detailed weak lensing mapping of substructure in clusters. Obviously, the source galaxies must 
be well resolved to allow an accurate shape measurement. Barkana & Loeb (2000a) estimated the 
size distribution of galactic disks (see § |B.3| ) and showed that with its expected ~ 0706 resolution, 
NGST should resolve most galaxies even at z ~ 10. 

The probability for strong gravitational lensing depends on the abundance of lenses, their 
mass profiles, and the angular diameter distances among the source, the lens and the observer. The 
statistics of existing lens surveys have been used at low redshifts to constrain the cosmological con- 
stant (for the most detailed work see Kochanek 1996, and references therein), although substantial 
uncertainties remain regarding the luminosity function of early-type galaxies and their dark matter 
content. Given the early stage of observations of the redshift evolution of galaxies and their dark 
halos, a theoretical approach based on the Press-Schechter mass function can be used to estimate 
the lensing rate. This approach has been used in the past for calculating lensing statistics at low 
redshifts, with an emphasis on lenses with image separations above 5" (Narayan & White 1988; 
Kochanek 1995; Maoz et al. 1997; Nakamura & Suto 1997; Phillips, Browne, & Wilkinson 2001; 
Ofek et al. 2001) or on the lensing rates of supernovae (Porciani &; Madau 2000; Marri et al. 2000). 

The probability for producing multiple images of a source at a redshift zs, due to gravitational 
lensing by lenses with density distributed as in a singular isothermal sphere, is obtained by inte- 
grating over lens redshift zl the differential optical depth (Turner, Ostriker, & Gott 1984; Fukugita 
et al. 1992) 

dr = 16vr 3 n (£f (1 + z L f (^f^)' |f ^ , (105) 

in terms of the comoving density of lenses n, velocity dispersion a, look-back time t, and angular 
diameter distances D among the observer, lens and source. More generally we replace na 4 by 

{n^)=J d ^^a\M,z L )dM, (106) 

where dn/dM is the Press-Schechter halo mass function. It is assumed that <j(M, z) = V C (M, z)/y/2 
and that the circular velocity V C (M, z) corresponding to a halo of a given mass is given by equa- 
tion (H). 

The ACDM model (with f2 m = 0.3) yields a lensing optical depth (Barkana & Loeb 2000a) of 
~ 1% for sources at zs = 10. The fraction of lensed sources in an actual survey is enhanced, however, 
by the so-called magnification bias. At a given observed flux level, unlensed sources compete with 
lensed sources that are intrinsically fainter. Since fainter galaxies are more numerous, the fraction 
of lenses in an observed sample is larger than the optical depth value given above. The expected 
slope of the luminosity function of the early sources (§|T^) suggests an additional magnification 
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bias of order 5, bringing the fraction of lensed sources at z$ = 10 to ~ 5%. The lensed fraction 
decreases to ~ 3% at z = 5. With the magnification bias estimated separately for each source 
population, the expected number of detected multiply-imaged sources per field of view of NGST 
(which we assume to be 4' x 4') is roughly 5 for z > 10 quasars, 10 for z > 5 quasars, 10 for z > 10 
galaxies, and 100 for z > 5 galaxies. 

High-redshift sources will tend to be lensed by galaxies at relatively high redshifts. In Figure |39| 
(adopted from Figure 2 of Barkana & Loeb 2000a) we show the lens redshift probability density 
p(zl), defined so that the fraction of lenses between zl and zl + dzi is p{z^)dz^. We consider a 
source at zs = 5 (solid curve) or at z$ = 10 (dashed curve). The curves peak around zl = 1, but 
in each case a significant fraction of the lenses are above redshift 2: 20% for zs = 5 and 36% for 
z s = 10. 




Fig. 39. — Distribution of lens redshifts for a fixed source redshift, for Press-Schechter halos in 
ACDM with £l m = 0.3 (adopted from Figure 2 of Barkana &: Loeb 2000a). Shown for a source at 
zs = 5 (solid curve) and for zs = 10 (dashed curve). The probability density p{zi) is shown, where 
the fraction of lenses between zl and zl + dzi is p^z^jdzi- 

The multiple images of lensed high-redshift sources should be easily resolvable. Indeed, image 
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separations are typically reduced by a factor of only 2-3 between zs = 2 and zs = 10, with the 
reduction almost entirely due to redshift evolution in the characteristic mass of the lenses. With 
a typical separation of 0.5-1" for z$ = 10, a large majority of lenses should be resolved given the 
NGST resolution of ~ 0"06. 

Lensed sources may be difficult to detect if their images overlap the lensing galaxy, and if the 
lensing galaxy has a higher surface brightness. Although the surface brightness of a background 
source will typically be somewhat lower than that of the foreground lens (Barkana & Loeb 2000a), 
the lensed images should be detectable since (i) the image center will typically be some distance 
from the lens center, of order half the image separation, and (ii) the younger stellar population and 
higher redshift of the source will make its colors different from those of the lens galaxy, permitting 
an easy separation of the two in multi-color observations. These helpful features are evident in 
the currently known systems which feature galaxy-galaxy strong lensing. These include two four- 
image 'Einstein cross' gravitational lenses and other lens candidates discovered by Ratnatunga et al. 
(1999) in the Hubble Space Telescope Medium Deep Survey, and a lensed three-image arc detected 
in the Hubble Deep Field South and studied in detail by Barkana et al. (1999). 



9. Observational Probes of the Epoch of Reionization 
9.1. Spectral Methods of Inferring the Reionization Redshift 

9.1.1. Cosmology with Lya Photons 

The scattering cross-section of the Lya resonance line by neutral hydrogen is given by (§23 of 
Peebles 1993) 

_ 3A|A| {v/v a f 

aa[U) Svr ^{v-v a Y + {A? a /A){v/v a r 1 j 

where A a = (8TT 2 e 2 f a /3m e c\ 2 a ) = 6.25 x 10 8 s" 1 is the Lya (2p -► Is) decay rate, f a = 0.4162 
is the oscillator strength, and A Q = 1216A and v a = (c/X a ) = 2.47 x 10 15 Hz are the wavelength 
and frequency of the Lya line. The term in the numerator is responsible for the classical Rayleigh 
scattering. 

We consider a source at a redshift z s beyond the redshift of reionization^] z re ion, and the 
corresponding scattering optical depth of a uniform, neutral IGM of hydrogen density tih,o(1 + Z ) S 
between the source and the reionization redshift. The optical depth is a function of the observed 
wavelength A b s , 

T(A obs ) = / dz— n H ,o(l + z) 6 a a [i/ ob s(l + *)] , (108) 



10 We define the reionization redshift to be the redshift at which the individual H II regions overlapped and most of 
the IGM volume was ionized. In most realistic scenarios, this transition occurs rapidly on a time-scale much shorter 



than the age of the universe (see §3.3.1). This is mainly due to the short distances between neighboring sources. 
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At wavelengths longer than Lya at the source, the optical depth obtains a small value; these 
photons redshift away from the line center along its red wing and never resonate with the line core 
on their way to the observer. Considering only the regime in which \v — v a \ 3> A Q , we may ignore 
the second term in the denominator of equation ( |107| ). This leads to an analytical result for the 
red damping wing of the Gunn-Peterson trough (Miralda-Escude 1998) 

A 



r(A 



obs J 



ill 2 1 



A 



3/2 
obs 



Ws) " '([(1 + *eioa)/(l + ^)]A ob i s ; 



for A obs > 1 



where t s is given in equation (|T|), and we also define 

Aobs 



A 



obs 



(1 + z s )X c 



and 



I(x) 



Q Q 

■' +l x 7/2+^5/2 + 3x 3/ 2+9x l/2 



1 



7 



In 



1 + x 1 ' 2 



X 



1/2 



(110) 



fill) 



(112) 



At wavelengths corresponding to the Lya resonance between the source redshift and the 
reionization redshift, (1 + ^reion)A Q < A ODS < (1 + z s )X a , the optical depth is given by equa- 
tion (|l]). Since r s ~ 10 5 , the flux from the source is entirely suppressed in this regime. Similarly, 
the Ly/3 resonance produces another trough at wavelengths (1 + z re \ on )\p < A < (1 + z s )\p, 
where \p = (27/32)A a = 1026 A, and the same applies to the higher Lyman series lines. If 
(1 + z s ) > 1.18(1 + -Zreion) then the Lya and the Ly/3 resonances overlap and no flux is transmitted 



in-between the two troughs (see Figure 40). The same holds for the higher Lyman-series resonances 



down to the Lyman limit wavelength of A c = 912A. 

At wavelengths shorter than A c , the photons are absorbed when they photoionize atoms of 
hydrogen or helium. The bound-free absorption cross-section from the ground state of a hydrogenic 
ion with nuclear charge Z and an ionization threshold huQ, is given by (Osterbrock 1974), 
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For neutral hydrogen, Z = 1 and z^h,o = (c/A c ) = 3.29 x 10 Hz (/i^h,o = 13.60 eV); for singly- 
ionized helium, Z = 2 and z^He n,o = 1-31 x 10 16 Hz (/iz^He n o = 54.42 eV). The cross-section for 
neutral helium is more complicated; when averaged over its narrow resonances it can be fitted to 
an accuracy of a few percent up to hv = 50 keV by the fitting function (Verner et al. 1996) 

-3 188 

^fe/,He i(v) = 9.492 x 10~ lfci cm 2 x (x - l) 2 + 4.158 y 1 """ ' ' " 



-1.953 



(l + 0.825y 



1/4 



(H5) 



- 112 - 



where x = [(z^/3.286 x 10 15 Hz) - 0.4434], y = x 2 + 4.563, and the threshold for ionization is 
!/ He I = 5.938 x 10 15 Hz (fa/He i,o = 24.59 eV). 

For rough estimates, the average photoionization cross-section for a mixture of hydrogen and 
helium with cosmic abundances can be approximated in the range of 54 < hu £ 10 3 eV as at,/ ~ 
(To(^/^h,o) _3 > where ao ~ 6 x 10~ 17 cm 2 (Miralda-Escude 2000). The redshift factor in the cross- 
section then cancels exactly the redshift evolution of the gas density and the resulting optical depth 
depends only on the elapsed cosmic time, t(z re i on ) — t(z s ). At high redshifts (equations (||) and (|l0|) 
in §^4]) this yields, 

cdt 



Tb/(A obs ) = / dz—n (l + zfa hi {v ohs (l + z 
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The bound-free optical depth only becomes of order unity in the extreme UV to soft X-rays, around 
hv ~ 0.1 keV, a regime which is unfortunately difficult to observe due to Galactic absorption 
(Miralda-Escude 2000). 

A sketch of the overall spectrum of a source slightly above the reionization redshift, i.e., with 
1 < [(1 + z s )/(l + Zreion)] < 1.18, is shown in Figure ^TJ. The transmitted flux between the 
Gunn-Peterson troughs due to Lya and Ly/3 absorption is suppressed by the Lya forest in the 
post-reionization epoch. Transmission of flux due to H II bubbles in the pre-reionization epoch is 
expected to be negligible (Miralda-Escude 1998). The redshift of reionization can be inferred in 
principle from the spectral shape of the red damping wing (Miralda-Escude & Rees 1998; Miralda- 
Escude 1998) or from the transmitted flux between the Lyman series lines (Haiman & Loeb 1999a). 
However, these signatures are complicated in reality by damped Lya systems along the line of 
sight or by the inhomogeneity or peculiar velocity field of the IGM in the vicinity of the source. 
Moreover, bright sources, such as quasars, tend to ionize their surrounding environment (Wood & 
Loeb 2000) and the resulting H II region in the IGM could shift the Lya trough substantially (Cen 
& Haiman 2000; Madau & Rees 2000). 

The inference of the Lya transmission properties of the IGM from the observed spectrum 
of high-redshift sources suffers from uncertainties about the precise emission spectrum of these 
sources, and in particular the shape of their Lya emission line. The first galaxies and quasars are 
expected to have pronounced recombination lines of hydrogen and helium due to the lack of dust 
in their interstellar medium (see § 4.1. 3| for more details). Lines such as H a or the He II 1640 A line 



should reach the observer unaffected by the intervening IGM, since their wavelength is longer than 
that of the Lya transition which dominates the IGM opacity (Oh 1999). However, as described 
above, the situation is different for the Lya line photons from the source. As long as z s > z Te i OJi , 
the intervening neutral IGM acts like a fog and obscures the view of the Lya line itself [in contrast 
to the situation with sources at z s < z re \ on , where most of the intervening IGM is ionized and 
only photons more energetic than Lya are suppressed by the Lya forest (see Figure §)]. Photons 
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which are emitted at the Lya line center have an initial scattering optical depth of ~ 10 5 in the 
surrounding medium. 

The Lya line photons are not destroyed but instead are absorbed and re-emitted^. Due to the 
Hubble expansion of the IGM around the source, the frequency of the photons is slightly shifted 
by the Doppler effect in each scattering event. As a result, the photons diffuse in frequency to 
the red side of the Lya resonance. Eventually, when their net frequency redshift is sufficiently 
large, they escape and travel freely towards the observer (see Figure |Il]). As a result, the source 
creates a faint Lya halo on the skyP|. The well-defined radiative transfer problem of a point source 
of Lya photons embedded in a uniform, expanding neutral IGM was solved by Loeb & Rybicki 
(1999). The Lya halo can be simply characterized by the frequency redshift relative to the line 



center, (y— u a ), which is required in order to make the optical depth from the source [equation (110)] 
equal to unity. At high redshifts, the leading term in equation (|110|) yields 



,!2 ,, / Mbh \ fl + Z s \ 3/2 



v± = 8.85 x 10^ Hz x °- T = ——^ . (117) 

This is the frequency interval over which the damping wing affects the source spectrum. A frequency 
shift of = 8.85 x 10 12 Hz relative to the line center corresponds to a fractional shift of {v-i,/v a ) = 
(v/c) = 3.6 x 10~ 3 or a Doppler velocity of v ~ 10 3 km s" 1 . The halo size is then defined by the 
corresponding proper distance from the source at which the Hubble velocity provides a Doppler 
shift of this magnitude, 

r. = 1.1 (^) Mpc. (118) 



O m /0.3. 

Typically, the Lya halo of a source at z s ~ 10 occupies an angular radius of ~ 15" on the sky and 



yields an asymmetric line profile as shown in Figures |41| and |42|. The scattered photons are highly 
polarized and so the shape of the halo would be different if viewed through a polarization filter 
(Rybicki & Loeb 1999). 



Detection of the diffuse Lya halos around bright high-redshift sources (which are sufficiently 
rare so that their halos do not overlap) would provide a unique tool for probing the distribution 
and the velocity field of the neutral IGM before the epoch of reionization. The Lya sources 



11 At the redshifts of interest, z s ~ 10, the low densities and lack of chemical enrichment of the IGM make the 
destruction of Lya photons by two-photon decay or dust absorption unimportant. 

12 The photons absorbed in the Gunn-Peterson trough are also re-emitted by the IGM around the source. However, 
since these photons originate on the blue side of the Lya resonance, they travel a longer distance from the source, 
compared to the Lya line photons, before they escape to the observer. The Gunn-Peterson photons are therefore 
scattered from a larger and hence dimmer halo around the source. The Gunn-Peterson halo is made even dimmer 
relative to the Lya line halo by the fact that the luminosity of the source per unit frequency is often much lower in 
the continuum than in the Lya line. 
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DETERMINING THE REIONIZATION REDSHIFT 




Fig. 40. — Sketch of the expected spectrum of a source at a redshift z s slightly above the reionization 
redshift z ve i on . The transmitted flux due to H II bubbles in the pre-reionization era and the Lya 
forest in the post-reionization era is exaggerated for illustration. 
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Fig. 41. — Loeb-Rybicki halos: Scattering of Lya line photons from a galaxy embedded in the 
neutral IGM prior to reionization. The line photons diffuse in frequency due to the Hubble expan- 
sion of the surrounding medium and eventually redshift out of resonance and escape to infinity. A 
distant observer sees a Lya halo surrounding the source, along with a characteristically asymmetric 
line profile. The observed line should be broadened and redshifted by about one thousand km s _1 
relative to other lines (such as H a ) emitted by the galaxy. 
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Fig. 42. — Monochromatic photon luminosity of a Lya halo as a function of frequency redshift, 
v = {va — v)lv±. The observed spectral flux of photons F(y) (in photons cm -2 s _1 Hz -1 ) from 
the entire Lya halo is F(y) = (Z(j/)/47r<i 2 j )(7V Q /^)(l + z s ) 2 where iV Q is the production rate of 
Lya photons by the source (in photons s _1 ), v = vv±/{l + z s ), and dj, is the standard luminosity 
distance to the source (from Loeb & Rybicki 1999). 
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serve as lamp posts which illuminate the surrounding H I fog. On sufficiently large scales where 
the Hubble flow is smooth and the gas is neutral, the Lya brightness distribution can be used 
to determine the cosmological mass densities of baryons and matter. Due to their low surface 
brightness, the detection of Lya halos through a narrow-band filter is much more challenging than 
direct observation of their sources at somewhat longer wavelengths. However, NGST might be able 
to detect the Lya halos around sources as bright as the quasar discovered by Fan et al. (2000) at 
z = 5.8 or the galaxy discovered by Hu et al. (1999) at z = 5.74, even if these sources were moved 
out to z ~ 10 (see §4 in Loeb & Rybicki 1999). The disappearance of Lya halos below a certain 
redshift can be used to determine ^reion* 



9.1.2. 21 cm Tomography of the Reionization Epoch 

The ground state of hydrogen exhibits hyperfine splitting involving the spins of the proton 
and the electron. The state with parallel spins (the triplet state) has a slightly higher energy 
than the state with anti-parallel spins (the singlet state). The 21 cm line associated with the 
spin-flip transition from the triplet to the singlet state is often used to detect neutral hydrogen in 
the local universe. At high redshift, the occurrence of a neutral pre-reionization IGM offers the 
prospect of detecting the first sources of radiation and probing the reionization era by mapping 
the 21 cm emission from neutral regions. While its energy density is estimated to be only a 1% 
correction to that of the CMB, the redshifted 21 cm emission should display angular structure as 
well as frequency structure due to inhomogeneities in the gas density field (Hogan & Rees 1979; 
Scott & Rees 1990), hydrogen ionized fraction, and spin temperature (Madau, Meiksin, & Rees 
1997). Some of the resulting signatures during the pre-overlap phase of reionization (§ |6.3.1 ) and 



during the overlap phase are discussed by Tozzi et al. (2000) and Shaver et al. (1999), respectively. 
Also, the 21 cm signatures have been explored in a numerical simulation by Gnedin & Ostriker 
(1997). Indeed, a full mapping of the distribution of H I as a function of redshift is possible in 
principle. Although detecting the presence of the largest H II regions may be within the reach 
of proposed instruments such as the Square Kilometer Array (hereafter SKA; see Taylor & Braun 
1999), these instruments may not have sufficient sensitivity at the sub-arcminute resolution that 
would be necessary for a detailed mapping. Moreover, serious technical challenges and problems 
due to foreground contamination must be overcome even for an initial detection of the reionization 
signal. 

The basic physics of the hydrogen spin transition is determined as follows (for a more detailed 
treatment, see Madau et al. 1997). The ground-state hyperfine levels of hydrogen tend to thermalize 
with the CMB background, making the IGM unobservable. If other processes shift the hyperfine 
level populations away from thermal equilibrium, then the gas becomes observable against the CMB 
in emission or in absorption. The relative occupancy of the spin levels is usually described in terms 
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of the hydrogen spin temperature T$, defined by 



" 3exp(-^ (1.19) 



n I T S 

where no and ri\ are the singlet and triplet hyperfine levels in the atomic ground state (n = 1), 
and = 0.07 K is defined by /cbT* = E21, where the energy of the 21 cm transition is E21 = 
5.9 x 10~ 6 eV, corresponding to a frequency of 1420 MHz. In the presence of the CMB alone, 
the spin states reach thermal equilibrium with = Xcmb = 2.73(1 + z) K on a time-scale of 
7*/(Tcmb^-io) — 3 x 10 5 (1 + z)~ l yr, where A w = 2.9 x 10~ 15 s _1 is the spontaneous decay rate 
of the hyperfine transition. This time-scale is much shorter than the age of the universe at all 
redshifts after cosmological recombination. 

The IGM is observable when the kinetic temperature Tk of the gas differs from Tomb and 
an effective mechanism couples T$ to Tk- Although collisional de-excitation of the triplet level 
(Purcell & Field 1956) is a possible mechanism, in the low-density IGM the dominant mechanism 
is scattering by Lya photons (Wouthuysen 1952; Field 1958). Continuum UV photons produced by 
early radiation sources redshift by the Hubble expansion into the local Lya line at a lower redshift. 
These photons mix the spin states via the Wouthuysen-Field process whereby an atom initially in 
the n = 1 state absorbs a Lya photon, and the spontaneous decay which returns it from n = 2 to 
n = 1 can result in a final spin state which is different from the initial one. Since the neutral IGM is 
highly opaque to resonant scattering, the shape of the radiation spectrum near Lya is determined 
by Tk (Field 1959), and the spin temperature is then a weighed mean of Tk and Tomb: 

= l + y a ' (120) 
where (if T5 » T*) the Lya pumping efficiency is 

P10T* ,,„,, 
Va = ■ (121) 

Here P10 is the indirect de-excitation rate of the triplet n = 1 state via the Wouthuysen-Field 
process, related to the total scattering rate P a of Lya photons by P±o = 4P a /27 (Field 1958). Thus 
the critical value of P a is given by the thermalization rate (Madau et al. 1997) 

Ah - 27A g° MB - 7.6 x ID" 12 (1±£) s- . (122) 

A patch of neutral hydrogen at the mean density and with a uniform T5 produces an optical 
depth at 21(1 + z) cm of 

assuming a high redshift z. Since the brightness temperature through the IGM is T& = TcMBe _r + 
Ts(l — e~ T ), the observed differential antenna temperature of this region relative to the CMB is 
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(Madau et al. 1997, with the O m dependence added) 

m - a + -r"c* - - o « — (M) - (l±£) - (^) , 

(124) 

where r < 1 is assumed and 5Ti, has been redshifted to redshift zero. In overdense regions, the 
observed dT^ is proportional to the overdensity, and in partially ionized regions STj, is proportional 
to the neutral fraction. Thus, if T$ S> Tomb then the IGM is observed in emission at a level that 
is independent of T$. On the other hand, if T$ <C Tqmb then the IGM is observed in absorption at 
a level that is a factor ~ Tqmb/Ts larger than in emission. As a result, a number of cosmic events 
are expected to leave observable signatures in the redshifted 21 cm line. 

Since the CMB temperature is only 2.73(l+z) K, even relatively inefficient heating mechanisms 
are expected to heat the IGM above Icmb well before reionization. Possible preheating sources 
include soft X-rays from early quasars or star-forming regions, as well as thermal bremsstrahlung 
from ionized gas in collapsing halos. However, even the radiation from the first stars may suffice for 
an early preheating. Only ~ 10% of the present-day global star formation rate is required (Madau 
et al. 1997) for a sufficiently strong Lya background which produces a scattering rate above the 
thermalization rate Pth- Such a background drives Ts to the kinetic gas temperature, which is 
initially lower than Tqmb because of adiabatic expansion. Thus, the entire IGM can be seen in 
absorption, but the IGM is then heated above Tcmb m ~ 10 s yr (Madau et al. 1997) by the atomic 
recoil in the repeated resonant Lya scattering. According to § |8.1| (also compare Gnedin 2000a), 
the required level of star formation is expected to be reached already at z ~ 20, with the entire 
IGM heated well above the CMB by the time overlap begins at z ~ 10. Thus, although the initial 
absorption signal is in principle detectable with the SKA (Tozzi et al. 2000), it likely occurs at 
<^ 100 MHz where Earth-based radio interference is highly problematic. 

As individual ionizing sources turn on in the pre-overlap stage of reionization, the resulting 
H II bubbles may be individually detectable if they are produced by rare and luminous sources such 
as quasars. If the H II region expands into an otherwise unperturbed IGM, then the expanding 
shell can be mapped as follows (Tozzi et al. 2000). The H II region itself, of course, shows neither 
emission nor absorption. Outside the ionized bubble, a thin shell of neutral gas is heated above 
the CMB temperature and shows up in emission. A much thicker outer shell is cooler than the 
CMB due to adiabatic expansion, but satisfies Ts = Tk and produces absorption. Finally, at large 
distances from the quasar, Ts approaches Icmb as the quasar radiation weakens. For a quasar 
with an ionizing intensity of 10 57 photons s -1 observed after ~ 10 7 yr with 2' resolution and 1 
MHz bandwidth, the signal ranges from -3 to 3 /xJy per beam (Tozzi et al. 2000). Mapping such 
regions would convey information on the quasar number density, ionizing intensity, opening angle, 
and on the density distribution in the surrounding IGM. Note, however, that an H II region which 
forms at a redshift approaching overlap expands into a preheated IGM. In this case, the H II region 
itself still appears as a hole in an otherwise emitting medium, but the quasar-induced heating is 
not probed, and there is no surrounding region of absorption to supply an enhanced contrast. 
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At redshifts approaching overlap, the IGM should be almost entirely neutral but with T$ S> 
Tqmb- In this redshift range there should still be an interesting signal due to density fluctuations. 
The same cosmic network of sheets and filaments that gives rise to the Lya forest observed at z <J 5 
should lead to fluctuations in the 21 cm brightness temperature at higher redshifts. At 150 MHz 
(z = 8.5), for observations with a bandwidth of 1 MHz, the root mean square fluctuation should 
be ~ 10 mK at 1', decreasing with scale (Tozzi et al. 2000). 

A further signature, observable over the entire sky, should mark the overlap stage of reioniza- 
tion. During overlap, the IGM is transformed from being a neutral, preheated and thus emitting 
gas, to being almost completely ionized. This disappearance of the emission over a relatively narrow 
redshift range can be observed as a drop in the brightness temperature at the frequencies corre- 
sponding to the latter stages of overlap (Shaver et al. 1999). This exciting possibility, along with 
those mentioned above, face serious challenges in terms of signal contamination and calibration. 
The noise sources include galactic and extragalactic emission sources, as well as terrestrial inter- 
ference, and all of these foregrounds must be modeled and accurately removed in order to observe 
the fainter cosmological signal (see Shaver et al. 1999 for a detailed discussion). For the overlap 
stage in particular, the sharpness of the spectral feature is the key to its detectability, but it may 
be significantly smoothed by inhomogeneities in the IGM. 



9.2. Effect of Reionization on CMB Anisotropies 

In standard cosmological models, the early universe was hot and permeated by a nearly uniform 
radiation bath. At z ~ 1200 the free protons and electrons recombined to form hydrogen atoms, 
and most of the photons last scattered as the scattering cross-section dropped precipitously. These 
photons, observed today as the Cosmic Microwave Background (CMB), thus yield a snapshot of the 
state of the universe at that early time. Small fluctuations in the density, velocity, and gravitational 
potential lead to observed anisotropies (e.g., Sachs & Wolfe 1967; Bennett et al. 1996) that can be 
analyzed to yield a great wealth of information on the matter content of the universe and on the 
values of the cosmological parameters (e.g., Hu 1995; Jungman et al. 1996). 

Reionization can alter the anisotropy spectrum, by erasing some of the primary anisotropy 
imprinted at recombination, and by generating additional secondary fluctuations that could be 
used to probe the era of reionization itself (see Haiman & Knox 1999 for a review). The primary 
anisotropy is damped since the rescattering leads to a blending of photons from initially different 
lines of sight. Furthermore, not all the photons scatter at the same time, rather the last scattering 
surface has a finite thickness. Perturbations on scales smaller than this thickness are damped since 
photons scattering across many wavelengths give canceling redshifts and blueshifts. If reionization 
occurs very early, the high electron density produces efficient scattering, and perturbations are 
damped on all angular scales except for the very largest. 

The optical depth to scattering over a proper length dl is dr = aTn e dl, where ot is the 
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Thomson cross-section and n e the density of free electrons. If reionization occurs instantaneously 
at redshift z, then the total scattering optical depth in ACDM is given by (e.g., §7.1.1 of Hu 1995) 



reionization of z = 5.8 (Fan et al. 2000), with r = 0.10 if z = 11.6 and r = 0.15 if z = 15.3. 
Recent observations of small-scale anisotropics (Lange et al. 2000; Balbi et al. 2000) revealed a 
peak in the power spectrum on a ~ 1° scale, as expected from the primary anisotropies in standard 
cosmological models. This indicates that the reionization damping, if present, is not very large, 
and the observations set a limit of r < 0.33 at 95% confidence (Tegmark & Zaldarriaga 2000) and, 
therefore, imply that reionization must have occurred at z ^ 30. 

However, measuring a small r from the temperature anisotropies alone is expected to be very 
difficult since the anisotropy spectrum depends on a large number of other parameters, creating a 
near-degeneracy which limits our ability to measure each parameter separately; the degeneracy of 
r with other cosmological parameters is due primarily to a degeneracy with the gravitational-wave 
background. However, Thomson scattering also creates net polarization for incident radiation which 
has a quadrupole anisotropy. This anisotropy was significant at reionization due to large-scale struc- 
ture which had already affected the gas distribution. The result is a peak in the polarization power 
spectrum on large angular scales of order tens of degrees (Zaldarriaga 1997). Although experiments 
must overcome systematic errors from the detector itself and from polarized foregrounds (such as 
galactic dust emission and synchrotron radiation), parameter estimation models (Eisenstein, Hu, & 
Tegmark 1999; Zaldarriaga, Spergel, & Seljak 1997) suggest that the peak can be used to measure 
even very small values of r: 2% for the upcoming MAP satellite, and 0.5% for the Planck satellite 
which will reach smaller angular scales with higher accuracy. 

Reionization should also produce additional temperature anisotropies on small scales. These 
result from the Doppler effect. By the time of reionization, the baryons have begun to follow dark 
matter potentials and have acquired a bulk velocity. Since the electrons move with respect to 
the radiation background, photons are given a Doppler kick when they scatter off the electrons. 
Sunyaev (1978) and Kaiser (1984) showed, however, that a severe cancellation occurs if the electron 
density is homogeneous. Opposite Doppler shifts on crests and troughs of a velocity perturbation 
combine to suppress the anisotropy induced by small-scale velocity perturbations. The cancellation 
is made more severe by the irrotational nature of gravitationally-induced flows. However, if the 
electron density varies spatially, then the scattering probability is not equal on the crest and on the 
trough, and the two do not completely cancel. Since a non-zero effect requires variation in both 
electron density and velocity, it is referred to as a second-order anisotropy. 

The electron density can vary due to a spatial variation in either the baryon density or the 
ionized fraction. The former is referred to as the Ostriker-Vishniac effect (Ostriker &: Vishniac 1986; 
Vishniac 1987). The latter depends on the inhomogeneous topology of reionization, in particular on 
the size of H II regions due to individual sources (§ |6.2[ ) and on spatial correlations among different 
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With our standard parameters (end of §[]]) this implies r = 0.037 at the current lower limit on 
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regions. Simple models have been used to investigate the character of anisotropics generated during 
reionization (Gruzinov & Hu 1998; Knox et al. 1998; Aghanim et al. 1996). The Ostriker-Vishniac 
effect is expected to dominate all anisotropics at small angular scales (e.g., Jaffe &: Kamionkowski 
1998), below a tenth of a degree, because the primary anisotropics are damped on such small scales 
by diffusion (Silk damping) and by the finite thickness of the last scattering surface. Anisotropics 
generated by inhomogeneous reionization may be comparable to the Ostriker-Vishniac effect, and 
could be detected by MAP and Planck, if reionization is caused by bright quasars with 10 Mpc-size 
ionized bubbles. However, the smaller bubbles expected for mini-quasars or for star-forming dwarf 
galaxies would produce an anisotropy signal which is weaker and at smaller angular scales, likely 
outside the range of the upcoming satellites (see, e.g., Haiman & Knox 1999 for discussion). Gnedin 
& Jaffe (2000) used a numerical simulation to show that, in the case of stellar reionization, the 
effect on the CMB of patchy reionization is indeed sub-dominant compared to the contribution of 
non-linear density and velocity fluctuations. Nevertheless, a signature of reionization could still be 
detected in future measurements of CMB angular fluctuations on the scale of a few arcseconds (see 
also Bruscoli et al. 2000, who find a somewhat higher power spectrum due to patchy reionization). 

9.3. Remnants of High-Redshift Systems in the Local Universe 

At the end of the reionization epoch, the heating of the IGM resulted in the photo-evaporation 
of gas out of halos of circular velocity V c above ~ 10-15 km s~ 1 (§ |6.4j ). The pressure of the hot gas 
subsequently shut off gas infall into even more massive halos, those with V c ~ 30 km s _1 (§|6.5|). 
Thus, the gas reservoir of photo-evaporating halos could not be immediately replenished. Some 
dwarf galaxies which were prevented from forming after reionization could have eventually collected 
gas at z = 1-2, when the UV background flux declined sufficiently (Babul &: Rees 1992; Kepner, 
Babul, & Spergel 1997). However, Kepner et al. (1997) found that even if the ionizing intensity 
J21 declines as (1 + z) A below z = 3, only halos with V c ^> 20 km s -1 can form atomic hydrogen 
by z = 1, and V c ;> 25 km s _1 is required to form molecular hydrogen. While the fact that the 
IGM was reionized has almost certainly influenced the abundance and properties of dwarf galaxies 
observed today, the exact manifestations of this influence and ways to prove that they occurred 
have not been well determined. In this section we summarize recent work on this topic, which 
should remain an active research area. 

The suppression of gas infall mentioned above suggests that the abundance of luminous halos 
as a function of circular velocity should show a break, with a significant drop in the abundance 
below V c = 30 km s . Such a drop may in fact be required in order to reconcile the ACDM 
model with observations. Klypin et al. (1999) and Moore et al. (1999) found that the abundance 
of halos with V c ~ 10-30 km s _1 in numerical simulations of the Local Group environment is 
higher by an order of magnitude than the observed dwarf galaxy abundance. The predicted and 
observed abundances matched well at V c > 50 km s . Bullock et al. (2000a) considered whether 
photoionization can explain the discrepancy at the low-mass end by preventing dark matter halos 
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from forming stars. They assumed that a sub-halo in the Local Group can host an observable 
galaxy only if already at reionization its main progenitor contained a fraction / of the final sub- 
halo mass. Using semi-analytic modeling, they found a close match to the observed circular velocity 
distribution for z ve [ on = 8 and / = 0.3. 

These results neglect several complications. As mentioned above, halos with V c ;> 20 km s _1 
may be able to accrete gas and form stars once again at z <^ 1. Any accreted gas at low redshift 
could have been previously enriched with metals and molecules, thus enabling more efficient cooling. 
On the other hand, if the progenitors had a V c <^ 17 km s -1 at 2: rc i on then they were not able to cool 
and form stars, unless molecular hydrogen had not been dissociated (§ |3.3[) . In order to reconcile 
the photoionization scenario with the recent episodes of star formation deduced to have occurred 
in most dwarf galaxies (e.g., Mateo 1998; Grebel 1998), a continuous recycling of gas over many 
generations of stars must be assumed. This would mean that supernova feedback (§fO]) was unable 
to shut off star formation even in these smallest known galaxies. In addition, the existence of a 
large abundance of sub-halos may be problematic even if the sub-halos have no gas, since they 
would interact with the disk dynamically and tend to thicken it (Toth & Ostriker 1992; Moore et 
al. 1999; Velazquez & White 1999). 

However, the photoionization scenario is useful because it may be testable through other im- 
plications. For example, Bullock et al. (2001) used semi-analytic modeling to show that many 
subhalos which did form stars before reionization were tidally disrupted in the Milky Way's grav- 
itational field, and the resulting stellar streams may be observable. The formation of the Milky 
Way's stellar component has also been investigated by White & Springel (1999). They combined a 
scaled-down dark matter cluster simulation with semi-analytic prescriptions for star formation in 
halos, and showed that the oldest stars in the Milky Way should be located mostly in the inner 
halo or bulge, but they cannot be easily identified because the populations of old stars and of low 
metallicity stars are only weakly correlated. 

Gnedin (2000c) pointed out several observed features of dwarf galaxies that may be related 
to their high-redshift histories. First, almost all Local Group dwarf galaxies with measured star 
formation histories show a sharp decline in the star formation rate around 10 Gyr ago. Gnedin 
noted that this drop could correspond to the suppression of star formation due to reionization 
(§ |6.5[) if the measured ages of the old stellar populations are somewhat underestimated; or that it 
could instead correspond to the additional suppression caused by helium reionization (§ |6.3.2j ) at 
z ^ 3. He showed that if only the old stellar population is considered, then the Schmidt Law (§ |5.2|) 
implies that the luminosity of each dwarf galaxy, divided by a characteristic volume containing a 
fixed fraction of all the old stars, should be proportional to some power of the central luminosity 
density. This assumes that the present central luminosity density of old stars is a good measure 
of the original gas density, i.e., that in the core almost all the gas was transformed into stars, and 
feedback did not play a role. It also assumes that the gas density distribution was self-similar in all 
the dwarf galaxies during the period when they formed stars, and that the length of this period was 
also the same in all galaxies. These assumptions are required in order for the total stellar content 
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of each galaxy to be simply related to the central density via the Schmidt Law evaluated at the 
center. Taking ten galaxies which have well-measured star formation histories, and which formed 
most of their stars more than 10 Gyr ago, Gnedin found a correlation with a power law of 3/2, 
as expected from the Schmidt Law (§ |5,2| ). Clearly, the theoretical derivation of this correlation 
combines many simplistic assumptions. However, explaining the observed correlation is a challenge 
for any competing models, e.g., models where feedback plays a dominant role in regulating star 
formation. 

Barkana & Loeb (1999) noted that a particularly acute puzzle is presented by the very smallest 
galaxies, the nine dwarf spheroidals in the Local Group with central velocity dispersions a <^ 
10 km s -1 , including five below 7 km s _1 (see recent reviews by Mateo 1998 & van den Bergh 
2000). These galaxies contain old stars that must have formed at z ^ 2, before the ionizing 
background dropped sufficiently to allow them to form. There are several possible solutions to the 
puzzle of how these stars formed in such small halos or in their progenitors which likely had even 
smaller velocity dispersions. The solutions are that (i) molecular hydrogen allowed these stars to 
form at z > z re i on , as noted above, (ii) the measured stellar velocity dispersions of the dwarf galaxies 
are well below the velocity dispersions of their dark matter halos, or (iii) the dwarf galaxies did not 
form via the usual hierarchical scenario. 

One major uncertainty in comparing observations to hierarchical models is the possibility 
that the measured velocity dispersion of stars in the dwarf spheroidals underestimates the velocity 
dispersion of their dark halos. Assuming that the stars are in equilibrium, their velocity dispersion 
could be lower than that of the halo if the mass profile is shallower than isothermal beyond the 
stellar core. The velocity dispersion and mass-to-light ratio of a dwarf spheroidal could also appear 
high if the galaxy is non-spherical or the stellar orbits are anisotropic. The observed properties of 
dwarf spheroidals require a central mass density of order O.IMq pc -3 (e.g., Mateo 1998), which is 
~ 7 x 10 5 times the present critical density. Thus, only the very inner parts of the halos are sampled 
by the central velocity dispersion. Detailed observations of the velocity dispersion profiles of these 
galaxies could be used to determine the circular velocity of the underlying halo more reliably. 

A cosmological scenario for the formation of dwarf spheroidal galaxies is favored by the fact 
that they are observed to be dark matter dominated, but this may not rule out all the alternatives. 
The dwarf dark halos may have formed at low redshift by the breakup of a much larger galaxy. 
Under this scenario, gas forming stars inside the large parent galaxy would have been unaffected 
by the photoionizing background. At low redshift, this galaxy may have collided with the Milky 
Way or come close enough to be torn apart, forming at least some of the dwarf spheroidal systems. 
Simulations of galaxy encounters (Barnes &; Hernquist 1992; Elmegreen, Kaufman, & Thomasson 
1993) have found that shocks in the tidal tails trigger star formation and lead to the formation of 
dwarf galaxies, but these galaxies contain only small amounts of dark matter. However, the initial 
conditions of these simulations assumed parent galaxies with a smooth dark matter distribution 
rather than clumpy halos with dense sub-halos inside them. As noted above, simulations (Klypin 
et al. 1999; Moore et al. 1999) suggest that galaxy halos may have large numbers of dark matter 
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satellites, and further simulations are needed to test whether these subhalos can capture stars which 
form naturally in tidal tails. 

A common origin for the Milky Way's dwarf satellites (and a number of halo globular clusters), 
as remnants of larger galaxies accreted by the Milky Way galaxy, has been suggested on independent 
grounds. These satellites appear to lie along two (e.g., Majewski 1994) or more (Lynden-Bell & 
Lynden-Bell 1995, Fusi-Pecci et al. 1995) polar great circles. The star formation history of the dwarf 
galaxies (e.g., Grebel 1998) constrains their merger history, and implies that the fragmentation 
responsible for their appearance must have occurred early in order to be consistent with the variation 
in stellar populations among the supposed fragments (Unavane, Wyse, & Gilmore 1996; Olszewski 
1998). Observations of interacting galaxies (outside the Local Group) also demonstrate that "tidal 
dwarf galaxies" do indeed form (e.g., Due & Mirabel 1997; Hunsberger, Charlton, & Zaritsky 1996). 

10. Challenges for the Future 

When and how did the first stars and black holes form and when and how did they ionize most 
of the gas in the universe? In this review we have sketched the first attempts to answer these 
questions and the basic physical principles that underlie these attempts. The coming decade will 
likely be marked by major advances in our ability to make theoretical predictions in an attempt to 
answer these questions, and will culminate with the launch of NGST, a telescope which is ideally 
suited for testing these predictions. At about the same time, the Planck satellite (and perhaps 
MAP before it) is expected to directly infer the reionization redshift from measurements of the 
CMB polarization power spectrum on large angular scales. Also in about a decade, next-generation 
arrays of radio telescopes may detect the 21 cm emission from the pre-reionization, neutral warm 
IGM. The difficult questions just mentioned will receive their ultimate answers from observations, 
but it surely is fun to try to find the answers theoretically in advance, before we can deduce them 
by looking through our most technologically-advanced telescopes. 
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